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Introduction 

Geometric algebra is the Clifford algebra of a finite dimensional vector space over real scalars 
cast in a form most appropriate for physics and engineering. This was done by David Hestenes 
(Arizona State University) in the 1960’s. From this start he developed the geometric calculus 
whose fundamental theorem includes the generalized Stokes theorem, the residue theorem, and 
new integral theorems not realized before. Hestenes likes to say he was motivated by the fact 
that physicists and engineers did not know how to multiply vectors. 


Researchers at Arizona State and Cambridge have applied these developments to classical me- 
chanics,quantum mechanics, general relativity (gauge theory of gravity), projective geometry, 
conformal geometry, etc. 
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Chapter 1 


Basic Geometric Algebra 


1.1 Axioms of Geometric Algebra 


Let V (p,q) be a finite dimensional vector space of signature (p,q)! over #. Then Va,b,c € V 
there exists a geometric product with the properties - 


(ab)c = a(bc) 
a(b+c)=ab+ac 
(a + b)c = ac+ be 

aaEe RK 


if 
2 = gee 
If a“ £0 then a = at 


1.2 Why Learn This Stuff? 


The geometric product of two (or more) vectors produces something “new” like the /—1 with 
respect to real numbers or vectors with respect to scalars. It must be studied in terms of its 
effect on vectors and in terms of its symmetries. It is worth the effort. Anything that makes 
understanding rotations in a N dimensional space simple is worth the effort! Also, if one proceeds 


'To be completely general we would have to consider V (p,q,r) where the dimension of the vector space is 
n=p+t+q+r and p, g, and r are the number of basis vectors respectively with positive, negative and zero squares. 


i 
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on to geometric calculus many diverse areas in mathematics are unified and many areas of physics 
and engineering are greatly simplified. 


1.3 Inner, -, and outer, A, product of two vectors and their basic 
properties 


The inner (dot) and outer (wedge) products of two vectors are defined by 


a-b= 5 (ab + ba) (1.1) 
1 
a\b= 5 (ab — ba) (1.2) 
ab=a-b+aAb (13) 
aN\b=—bAa (1.4) 
c=atb 
c= (a+b) 
C=a?+ab+bat+b? (1.5) 
2a-b=C—-a—0? 
a-bER 
a-b= |a| |b| cos (0) if a?,b? > 0 (1.6) 


Orthogonal vectors are defined by a-b = 0. For orthogonal vectors a A\ b = ab. Now compute 
(a Ab)? 


Gh —(aAb)(bAa) (1.7) 
—(ab—a-b) (ba—a-b) (1.8) 

— (abba — (a- b) (ab + ba) + (a- b)”) (1.9) 

— (ab? — (a-b) ") (1.10) 

= —a’b’ (1 — cos” (6)) (1.11) 

= —a’b’ sin? (0) (1.12) 


Thus in a Euclidean space, a?,b? > 0, (a. b)” <0 and aA b is proportional to sin (0). If ey and 


€, are any two orthonormal unit vectors in a Euclidean space then (eye fe = —1. Who needs 
the //—1? 


1.4. OUTER, A, PRODUCT FOR R VECTORS IN TERMS OF THE GEOMETRIC PRODUCT9 
1.4 Outer, A, product for r Vectors in terms of the geometric product 


U1 ..bp 


Define the outer product of r vectors to be (€{/"7"" is the mixed permutation symbol) 


aA...Aa= “ pS Se Ree Fa (1.13) 
Hise 
Thus 
i No (GRA Oct tip 
DT CIT ie I ig De I SOG Rows IV Op (1.14) 
and 


ayA...Na; NGj41 NA... N Ap = 
—aA...Naj41 AN aj... ANG (1.15) 


The outer product of r vectors is called a blade of grade r. 


1.5 Alternate Definition of Outer, A, product for r Vectors 


Let €1,€2,...,e, be an orthogonal basis for the set of linearly independent vectors a1, d2,..., Gy 
so that we can write 
= S> aue; (1.16) 
J 
Then 


a102...a, = ( ) ones] ( ) ones] nar’ ( ) on) 
J J2 Jr 


) O17, A2jq + + + Arj, Cj, O72 ++ + Ej, (1.17) 
ji a2"), Jr 


Now define a blade of grade n as the geometric product of n orthogonal vectors. Thus the product 
€5,€jy ---€;, IM equation 1.17 could be a blade of grade r, r — 2, r — 4, etc. depending upon the 
number of repeated factors. 
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If there are no repeated factors in the product we have that 


— pdledr 
ee oF eS ee (1.18) 


Tr 


Due to the fact that interchanging two adjacent orthogonal vectors in the geometric product will 
reverse the sign of the product and we can define the outer product of r vectors as 


ay N\...\Gp = » Caer ts 10 pj, €1 . + Ep (1.19) 
Fly Jr 
= det (a) €1...e, (1.20) 


Thus the outer product of r independent vectors is the part of the geometric product of the r 
vectors that is of grade r. Equation 1.19 is equivalent to equation 1.13. This can be proved by 
substituting equation 1.17 into equation 1.13 to get 


1 ee 
ay/A...A\Q,r = A S- S- Bie Ona ee Min jp EF ++ + CF, (1.21) 


BM ygtr Jly rode 


1 fine sah <b) 
dy-br JL Ir 
wal ) ) ey Ey op Wiggy ++ CipjpE1- ++ er (1.22) 


UW yeeybr Jl yey dr 


i bascit 42, 2 

=s ys ef Pret” det (a) €1 ... ep (1.23) 
es Jl yedr 

= det (a) e1...e, (1.24) 


We go from equation 1.22 to equation 1.23 by noting that > er ay, i, «. + i,j, 18 just det (a) 
i ate 


j\-Jr gives the correct sign for the determi- 


with the columns permuted. Multiplying det (a) by e417 
nant with the columns permuted. 


If e,,...,@, is an orthonormal basis for vector space the unit psuedoscalar is defined as 
Pere, (1.25) 
In equation 1.24 let r = n and the aj,...,a, be another orthonormal basis for the vector space. 


Then we may write 
Q1... Qn = det (a) €,... en (1.26) 


Since both the a’s and the e’s form orthonormal bases the matrix a is orthogonal and det (a) = 
+1. All psuedoscalars for the vector space are identical to within a scale factor of +1.? Likewise 
a, \.../ Qn is equal to I times a scale factor. 


It depends only upon the ordering of the basis vectors. 


1.6. USEFUL RELATION’S vt 


1.6 Useful Relation’s 


1. For a set of r orthogonal vectors, e€1,..., ey 
CL az INC eo = Cys he (1.27) 
2. For a set of r linearly independent vectors, a,,...,a,, there exists a set of r orthogonal 
vectors, €1,.-.,@,, such that 
Oty Os Chae (1.28) 
If the vectors, a1,...,a,, are not linearly independent then 
Gr hea 0 (1.29) 


The product a; A... Aa, is call a “blade” of grade r. The dimension of the vector space is the 
highest grade any blade can have. 


1.7 Projection Operator 


A multivector, the basic element of the geometric algebra, is made of of a sum of scalars, vectors, 
blades. A multivector is homogeneous (pure) if all the blades in it are of the same grade. The 
grade of a scalar is 0 and the grade of a vector is 1. The general multivector A is decomposed 
with the grade projection operator (A),. as (N is dimension of the vector space): 


A= S*(A) 


r=0 


(1.30) 


Tr 


As an example consider ab, the product of two vectors. Then 
ab = (ab), + (ab). (1.31) 
We define (A) = (A), for any multivector A 


1.8 Basis Blades 


The geometric algebra of a vector space, V (p,q), is denoted G (p,q) or G(V) where (p,q) is the 
signature of the vector space (first p unit vectors square to +1 and next g unit vectors square to 
—1, dimension of the space is p+ q). Examples are: 
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p q_ Type of Space 
3 0 3D Euclidean 
1 3. Relativistic Space Time 
4 1 3D Conformal Geometry 
If the orthonormal basis set of the vector space is €1,...,ey, the basis of the geometric algebra 


(multivector space) is formed from the geometric products (since we have chosen an orthonormal 
basis, ei? = +1) of the basis vectors. For grade r multivectors the basis blades are all the 
combinations of basis vectors products taken r at a time from the set of N vectors. Thus the 
number basis blades of rank r are oe My the binomial expansion coefficient and the total dimension 


N 
- 


of the multivector space is the sum of ( ) over r which is 2%. 


1.8.1 G(3,0) Geometric Algebra (Euclidian Space) 


The basis blades for G (3,0) are: 


Grade 
0 1 2 3 
1 €y E€1€2 E€1E2€3 
€2  €1€3 
€3  €2€3 


The multiplication table for the G (3,0) basis blades is 


1 Ee, €2 €3 e€1,€2 €1€3 €2€3 €1€2€3 

1 1 ial €2 €3 e€1€2 €1€3 €2€3 €1€2€3 
El ial 1 €1€2 €1€3 e€2 €3 €1€2€3 €2€3 

€2 €2 —e€1€2 1 e€2€3 —e] —-€1€9€3 €3 —€1€3 
€3 €3 —€1€3 —E9E3 1 €1€2€3 —€1 —eé2 €1€2 
e€1,€2 e€1,€2 —e2 ey, €1€2€3 —1 —€9€3 €1€3 ES 
€1€3 E€1€3 —€3 —€1€2€3 et €2€3 —1 —€1€2 e€2 
€2€3 €2€3 ~=€1€2€3 —€3 e2 —€1€3 e€1e2 —1 —€] 
€1€2€3 | €1€2€3 e€2€3 —€1€3 €1€2 —€3 e€2 el —1 


Note that the squares of all the grade 2 and 3 basis blades are —1. The highest rank basis blade 
(in this case €,€2e3) is usually denoted by J and is called the pseudoscalar. 


1.8. BASIS BLADES 13 
1.8.2 G(1,3) Geometric Algebra (Spacetime) 
The multiplication table for the G (1,3) basis blades is 
1 Yo a SP. 3 yor yoY2 172 
i 1 Yo onl a 13 Yov1 YoY2 172 
Yo Yo 1 yoV1 72 Yo73 V1 Y2 yov172 
pa onl —771 = y172 13 Yo oye =o 
Y2 273 —07Y2 — S72 = 7273 Yo" 72 Yo V1 
3 3 —V073 =VAY3 —Y273 mk: YOYLY3 YOY2Y3 Y1273 
YOY1 YOY1 —~Y1 —Y0 YoY1Y2 YOY1LY3 1 ~V1Y2 —Yo7Y2 
YoY2 YoY2 ~Y2 —~YOY1 72 —Y0 YoY2Y3 V172 1 YoY1 
Y1Y2 V172 YOY1Y2 2 ~Y1 Y1Y273 YoY2 —~YOY1 =I, 
YoY3 YoY3 —~Y3 —~YOY1 73 —VoY273 —~Y0 Y13 273 YOV1Y273 
173 V173 YOY1LY3 3 —~YV1 7273 = YoY3 —VV1 7273 —Y273 
7273 7273 YoY2 3 Y1Y273 3 a YOY1 7273 YOY3 YAY3 
Yov172 YoY172 Y1Y2 YyoY2 — 71 YoY1 7273 Y2 =I =70 
YOY1LY3 YOY1Y3 Y173 YoY3 —VO71 7273 —~Y7Y1 3 ~V1Y273 —VOV273 
Yo7273 yo Y273 7273 YoY1 7273 YoY3 =Vory2 17273 13 YoV1 73 
Y1Y293 V17293 —Yo71 7273 —Y273 YAY3 = YLY2 YOY2Y3 —YVOV1 3 —~Y3 
YOY1Y273 | YoYV1Y2Y3 —V17273 —Yo7273 YOYV1Y3 —V71 72 2°73 —~V1Y3 —Y073 
YoY3 W173 7273 Yov1y2 Yor 73 yov273 Y17273—_- YoV1 7273 
i YoY3 Y1Y3 273 YOY1Y2 YOY1LY3 Yo7273 VAV273 = YOV1Y2Y3 
Yo 3 YOY1Y3 YOY2Y3 VW72 173 1273 YOVIV2Y3. «V1 7273 
71 —YY1Y3 —%73 V1 7273 YoY2 YoY3 —YV17273 =—J273 yo v273 
2 —~YY2Y3 —V172Y3 —~Y3 —~V71 YoY1 723 YOY3 Y1VY3 —YVOYV1Y3 
13 Yo V1 Y2 —YV172Y3 =—YOY1 — 072 — 172 Yov172 
YoY1 —~YV13 —~YoY3 YoY1Y2Y3 2 3 —V17273 —~YoY2Y3 Y23 
YoY2 —2Y3 —Yo71 7273 —V073 ~YV1 Y1Y293 3 YOYLY3 —V173 
V1 72 YoY1 723 7273 —V13 me) YoY273 —VOV1 3 —~Y3 —Y073 
YoY3 1 YOY1 YoY2 V1 7273 meg (ik —~Y2 —Vo71 72 Y1Y2 
Y1Y3 —~YO7Y1 =k 172 —~Vo7273 —Y0 YoY1Y2 2 YoY2 
7273 —Vo72 —V1 72 = YOY1Y3 —~YOY1Y2 0 a 2 —~YOY1 
YOoY1Y2 Y17293 YOY273 —~YOY1 73 = 273 —~YV13 —V073 —Y3 
YoV173 =i =o yor172 Aaa = N72 Yo"2 72 
Yo1273 =2 —V172 —Y0 VIS —1172 =I — V1 =i 
Vi9273 yorv172 “YD =i Ys —Y072 Yo71 1 —70 
YoY1 9273 172 YoY2 —Y71 3 = 2 ga Yo =s 
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1.9 Reflections 


We wish to show that a,v € V > ava € VY and v is reflected about a if a? = 1. 


Figure 1.1: Reflection of Vector 


1. Decompose v = vj +v_ where vy is the part of v parallel to a and v_ is the part perpendicular 
to a. 


2. av = avy + av_ = va — v1a since a and v_ are orthogonal. 
3. ava = a? (v = v1) is a vector since a? is a scalar. 
4. ava is the reflection of v about the direction of a if a? = 1. 


5. Thus a,...a,va;,...a, € Y and produces a composition of reflections of v if a? = +--+ = 
2 
Goal, 


1.10. ROTATIONS 15 


1.10 Rotations 


1.10.1 Definitions 


First define the reverse of a product of vectors. If R = a,...as then the reverse is Rt = 
(Giese as)! =a,...a,, the order of multiplication is reversed. Then let R = ab so that 


RR' = (ab)(ba) = ab?a = a*b? = R'R (132) 

Let RR! =1 and calculate (RuRt)’, where v is an arbitrary vector. 
(RuR')* = RoR'RuRt = Rv? Ri =v? RR = 0? (1.33) 
Thus RvR! leaves the length of v unchanged. Now we must also prove Rv, R'- RvzRt = v1 - vo. 


Since Rv, Rt and RvR! are both vectors we can use the definition of the dot product for two 
vectors 


RvR RvR! = = (Ru, Ri Rup R! + Rv2R' Rv, RB") 


(Rviv2R' + Rv, R') 


NO eed NON Maes 


= ak (vyv9 + V2U1) Ri 
=R (v1 : V2) Ri 
=Ui,° voRR' 


= U1° V2 


Thus the transformation RvR! preserves both length and angle and must be a rotation. The 
normal designation for R is a rotor. If we have a series of successive rotations R,, Ro,..., Ry to 
be applied to a vector v then the result of the k rotations will be 


Ry, Rp ox. RivRLRS..: BR} 


Since each individual rotation can be written as the geometric product of two vectors, the com- 
position of k rotations can be written as the geometric product of 2k vectors. The multivector 
that results from the geometric product of r vectors is called a versor of order r. A composition 
of rotations is always a versor of even order. 
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1.10.2 General Rotation 


The general rotation can be represented by R = e2" where u is a unit bivector in the plane of 
the rotation and @ is the rotation angle in the plane.? The two possible non-degenerate cases are 
w=+1 ; ; 
4 7 An ae 0 : @ 
Si. (Euclidean plane) u . 1: cos (3) + usin (5) ; (1.34) 
(Minkowski plane) wu? =1: cosh (%) + wsinh (§) 


Decompose v = vy + (v — vI)) where vj is the projection of v into the plane defined by u. Note 
that v — vy is orthogonal to all vectors in the u plane. Now let u = e,ey where ey is parallel to 
vy and of course e, is in the plane u and orthogonal to ey. v — vy anticommutes with e) and e_ 
and v) anticommutes with e, (it is left to the reader to show RR! = 1). 


1.10.3. Euclidean Case 


For the case of u? = —1 


vy = |vyl er 


Figure 1.2: Rotation of Vector 


RuRt = (cos (5) + ee) sin (5)) (vy + (v — v)) (cos (5) + ee, sin @) 


Since v — vy anticommutes with e; and e, it commutes with R and 


RuR' = Rv Ri a (v = vI)) (1.35) 


lo e) 
3eA is defined as the Taylor series expansion e4 = y =n where A is any multivector. 
j=0 7" 
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So that we only have to evaluate 


RvR! = (cos (5) + e,e sin (5) oT (cos (5) + eye, sin (5)) (1.36) 


Since vj = |vj| e\ 
Rv Rt = |vj| (cos (9) ey + sin (4) e,) (1.37) 


and the component of v in the u plane is rotated correctly. 


1.10.4 Minkowski Case 


For the case of u? = 1 there are two possibilities, on > 0 or on <0. In the first case ej = 1 and 
e{ = —1. In the second case ej = —1 and ej = 1. Again v — uv, is not affected by the rotation 
so that we need only evaluate 


0 : 0 0 . 0 
RvR = (cost (5) + e,e) sinh (5)) vy (cost (5) + eye, sinh (5) 


Note that in this case || — on and 
2 i 
tf ope U: |v | (cosh (4) e) + sinh (9) e,) 
a { on <0: [uy] (cosh (4) ey — sinh (4) e,) ie) 


1.11 Expansion of geometric product and generalization of - and / 


If A, and B, are respectively grade r and s pure grade multivectors then 
A,B, = (Ar Bs) ips at (Ar Bs) ip si42 Bees ae (Ar Bs) min(r+s,2N—(r+s)) (1.39) 
A,: B, = (ABs) is (1.40) 
A, \ B, = (A,Bs) (1.41) 


Thus if r+s > N then A, A B, = 0, also note that these formulas are the most efficient way of 
calculating A,-B, and A, A B,. Using equations 1.28 and 1.39 we can prove that for a vector a 
and a grade r multivector B, 


r+s 


BS 5 (aB, 22) Re (1.42) 
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1 
aN By = 5 (aB, + (-1)" Bya) (1.43) 


If equations 1.42 and 1.43 are true for a grade r blade they are also true for a grade r multivector 
(superposition of grade r blades). By equation 1.28 let B, = e,...e, where the e’s are orthogonal 
and expand a 


a=ay+ Saye; (1.44) 
j=l 


where a, is orthogonal to all the e’s. Then* 


Tr 
OB. = y (—1)? *ajete; --+2)---ep + aser... ey 
j=l 


=a-B,+aAB, (1.45) 


Now calculate 


r 


Ba So (-1)" FajeFer Eye Op — (Sy dit ces 


i 
=(-1)"" (Sw aige ++ Gyr ++ €p — A141... “| 

jal 
= (-1)"""(a- B, —aAB,) (1.46) 


Adding and subtracting equations 1.45 and 1.46 gives equations 1.42 and 1.43. 


1.12 Duality and the Pseudoscalar 


If e1,...,@, is an orthonormal basis for the vector space the the pseudoscalar J is defined by 
meee (1.47) 


Since one can transform one orthonormal basis to another by an orthogonal transformation the 
I’s for all orthonormal bases are equal to within a +1 scale factor with depends on the ordering 
of the basis vectors. If A, is a pure r grade multivector (A, = (A,),.) then 


A,I = (A,1) (1.48) 


n—-r 


4 3 as 
€1.--€j7-1€j€j41---€p = €1--- €j7-1€j41--- Er 
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or A,/ is a pure n — rr grade multivector. Further by the symmetry properties of J we have 
TA, = (-1)°""" A, I (1.49) 


I can also be used to exchange the - and A products as follows using equations 1.42 and 1.43 


ORS : GAT (Ayr aa) (1.50) 
= 5 (aA-F — (1) (-1)""* A,al) ve) 
zs 5 (aA, elt Aca T (1.52) 
Bee (1.53) 


More generally if A, and B, are pure grade multivectors with r+s < n we have using equation 1.40 
and 1.48 


Ar + (Bsl) = (Ar Bs1) p_(n—s)| (1.54) 
UA, B.), 1 (1.56) 
= (Agi BB) 1 (1.57) 
Finally we can relate I to I’ by 
n(n—1) 

PSs F (1.58) 

1.13 Reciprocal Frames 
Let e€1,...,€, be a set of linearly independent vectors that span the vector space that are not 


necessarily orthogonal. These vectors define the frame (frame vectors are shown in bold face since 
they are almost always associated with a particular coordinate system) with volume element 


En, =e, A... \ €n (1.59) 
So that E, «x I. The reciprocal frame is the set of vectors e',...,e” that satisfy the relation 


eee 67, Vi = Leis tt (1.60) 
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The e’ are constructed as follows 
ei =(-1) Ve, Aen A... NGA... NenE,! (1.61) 


So that the dot product is (using equation 1.53 since E,' x I) 


e;-e) = (-1)"" e;- (ey Neg A...A€;A...N€nE,, ) (1.62) 
=(-1)"(eAerAesA...A€j;A... Aen) ED! (1.63) 
=0, Wit) (1.64) 

and 
e,-e =e: (e2 Pe aelicd e,E;,') (1.65) 
=(e,;Ae.A...A€,) E,' (1.66) 
= (1.67) 


1.14 Coordinates 


The reciprocal frame can be used to develop a coordinate representation for multivectors in 
an arbitrary frame e1,...,e@, with reciprocal frame e',...,e”. Since both the frame and it’s 
reciprocal span the base vector space we can write any vector a in the vector space as 


a=a'e;=a,e' (1.68) 


where if an index such as 7 is repeated it is assumes that the terms with the repeated index will 
be summed from 1 to n. Using that e; - e? = 67 we have 


aj=a:e; (1.69) 
a'=a-e' (1.70) 


In tensor notation a; would be the covariant representation and a‘ the contravariant representa- 
tion of the vector a. Now consider the case of grade 2 and grade 3 blades: 


e'-(aAb)=a-e'b—b-e'a 
e; (a: e'b—b- e'a) =ab—ba=2aAb 
e'-(aAbAc)=a-ebAc—b-eadct+c:-e'anb 
e; (a-e’bAc—b-e'adct+c-eaNb) =abAc—badc+canb=3aNbAc 
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for an r-blade A, we have (the proof is left to the reader) 
ee. A, =rA, (1.71) 
Since e;e’ = n we have 
ee’ \ A, =e; (e'A, —e'- A,) =(n—r)A, (1.72) 
Flipping e’ and A, in equations 1.71 and 1.72 and subtracting equation 1.71 from 1.72 gives 
e;A,e' = (—1)" (n — 2r) A, (i273) 
In Hestenes and Sobezyk (3.14) it is proved that 
(e* Iosad§ e"!) Nex ecKer) = OL be ae oer (1.74) 


so that the general multivector A can be expanded in terms of the blades of the frame and 
reciprocal frame as 


A= Ss" Aj;.ne’ Ne? Avs e8 (1.75) 
i<j<in<k 
where 
Bee = (ex, Ars A ej /\ e;) -A (1.76) 


The components Ajj... are totally antisymmetric on all indices and are usually referred to as the 
components of an antisymmetric tensor. 


1.15 Linear Transformations 


1.15.1 Definitions 


Let f be a linear transformation on a vector space f : ¥V > V with f (aa + 6b) =af (a)+ Gf (b) 
Va,b€ V anda, eR. Then define the action of f on a blade of the geometric algebra by 


f (ay Aw A Gy) = Ff (aq) Aes Af (ay) (1.77) 
and the action of f on any two A, B € G(V) by 


f (aA + BB) = af (A) + Bf (B) (1.78) 
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Since any multivector A can be expanded as a sum of blades f (A) is defined. This has many 
consequences. Consider the following definition for the determinant of f, det (f). 


f (D) = det (f)1 (1.79) 


First show that this definition is equivalent to the standard definition of the determinant (again 
€1,.--,@y is an orthonormal basis for V). 


f (er) = Ges (1.80) 


Then 
N N 
1203) — (> ont A... A (> cant 
sy=l1 sn=l1 
= S- Cites degC e wees (1.81) 
$15--53N 
But 
Cx oc Hogs See Re Cis SON (1.82) 
so that 
FISD en tae, acting h (1.83) 
ST yeeey SN 
or 
det (f) = S- eyo ay Gisys 2c ONes (1.84) 
Sy yeeey SN 


which is the standard definition. Now compute the determinant of the product of the linear 
transformations f and g 


det (fg) l= fg) 


= f(g) 

= f (det (g) F) 

= det (9) f 2) 

= det (g) det (f) I (1.85) 
det (fg) = det (f) det (g) (1.86) 


Do you have any idea of how miserable that is to prove from the standard definition of determi- 
nant? 
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1.15.2 Adjoint 


If F is linear transformation and a and b are two arbitrary vectors the adjoint function, F, is 
defined by 


a: F (b) =b- F(a) (1.87) 
From the definition the adjoint is also a linear transformation. For an arbitrary frame e),...,€n 
we have 


So that we can explicitly construct the adjoint as 


F(a)=e (e; . F (a) 
= e'(a- F (e,)) 
=e' (F (e;)-e’) a; (1.89) 
so that Fi; = F (e;) - e? is the matrix representation of F for the e;,...,€, frame. However 
F (a) =e' (F (e’) - e;) aj (1.90) 
so that the matrix representation of F is Fj; = F (e’)-e;. If the e1,...,e€, are orthonormal then 


e; =e? for all j and Fj; = Fj; exactly the same as the adjoint in matrices. 


Other basic properties of the adjoint are: 
F (a) = e'a- F (e;) = e'e;- F(a) = F (a) (1.91) 


and 


= G (F(a) (1.92) 


so that F = F and FG =GF. A symmetric function is one where F = F. As an example 
consider F'F 


FR=FF=FF (1.93) 
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1.15.3 Inverse 


Another linear algebraic relation in geometric algebra is 
2 If IA) 
1 A) = 


where f is the adjoint transformation defined by a- f(b) = b- f (a) Va,b € V and you have an 
explicit formula for the inverse of a linear transformation! 


VA EG(V) (1.94) 


1.16 Commutator Product 


The commutator product of two multivectors A and B is defined as 
1 
AxB= m (AB — BA) (1.95) 


An important theorem for the commutator product is that for a grade 2 multivector, Ay = (A)a, 
and a grade r multivector B, = (B),. we have 


AoB,. = Apo ix B,. + Ax B, + Apo : B, (1.96) 


From the geometric product grade expansion for multivectors we have 


A2B, = (A2 Br) 49 + (A2Br), + (A2Br)),-91 (1.97) 
Thus we must show that 
(ApB,), = AoxB; (1.98) 
Let e€1,...,€n be an orthogonal set for the vector space where B, = e,...e, and Ap = s Qimelem 
So we can write ‘ — 
Ag Be = ( S- once] x (€1...€,) (1.99) 
l<m=2 


Now consider the following three cases 
1. land m >r where ejeme,... Cr = €1.-- Cr€1em 


2. l<randm>~r where €je,€1...€ = —€1... Cr€]€m, 
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3. land m <r where ejé»,€1...€, = €1 ... €pCEm 


For case 1 and 3 €jém commute with B, and the contribution to the commutator product is 
zero. In case 3 €j€,, anticommutes with B, and thus are the only terms that contribute to 
the commutator. All these terms are of grade r and the theorem is proved. Additionally, the 
commutator product obeys the Jacobi identity 


Ax (BxC) = (AxB) xC' + Bx (AxC) (1.100) 


This is important for the geometric algebra treatment of Lie groups and algebras. 


Chapter 2 


Examples of Geometric Algebra 


2.1 Quaternions 


Any multivector A € G (3,0) may be written as 
A=a+a+B+I (2.1) 


where a, 8 € R, a € V(3,0), B is a bivector, and J is the unit pseudoscalar. The quaternions 
are the multivectors of even grades 


A=a+B (22) 


B can be represented as 
B=oi+ Bj+7k (2.3) 


where i = €9€3, j = e1e3, and k = e,é2, and 


Pep Hie Sika. (2.4) 


The quaternions form a subalgebra of G (3,0) since the geometric product of any two quaternions 
is also a quaternion since the geometric product of two even grade multivector components is a 
even grade multivector. For example the product of two grade 2 multivectors can only consist 
of grades 0, 2, and 4, but in G (3,0) we can only have grades 0 and 2 since the highest possible 
grade is 3. 


yas 
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2.2 Spinors 


The general definition of a spinor is a multivector, = € G (p,q), such that wow! € V(p,q) Vu € 
V (p,q). Practically speaking a spinor is the composition of a rotation and a dilation (stretching 
or shrinking) of a vector. Thus we can write 


wow! = pRuR' (2.5) 
where R is a rotor (RR = Ly; Letting U = R'w we must solve 
UvU' = pv (2.6) 


U must generate a pure dilation. The most general form for U based on the fact that the |.h.s 
of equation 2.6 must be a vector is 


U=a+ fl (2.7) 
so that 
UvUt = a*v+ af (Iv + vl") + 7 Ivl' = pu (2.8) 


(n—1)(n—2) 
2 


Using vit = (—1) Iv, vIt = (-1)""* Itv, and It = (—1)! we get 


(n=1)(n—2) 


a’v + af (1 - (yy) Iv + (—1)"t** Bu = pv (2.9) 


ip (aL (n= 2) 
2 


is even @ = 0 and a ¥ 0, otherwise a, 3 4 0. For the odd case 


vb = R(a+ Bl) (2.10) 


where p = a2 +(—1)"**"' 8? In the case of G (1,3) (relativistic space time) we have p = a? + 62, 
ped 


2.3 Geometric Algebra of the Minkowski Plane 


Because of Relativity and QM the Geometric Algebra of the Minkowski Plane is very important 
for physical applications of Geometric Algebra so we will treat it in detail. 
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Let the orthonormal basis vectors for the plane be yo and y, where yj = —y7 = 1.' Then the 
geometric product of two vectors a = apyo + 171 and b = boy + 6191 is 
ab = (aoyo + 4191) (boyo + 6141) (2.11) 
= agboy§ + a1b177 + (aob1 — a1b0) Yo" (2.12) 
_ agbo = a,b, + (ab = abo) I (2.13) 
so that 
a:-b= aobo = ay,o4 (2.14) 
and 
EKd= (aby = abo) I (2.15) 
and 
P= 71071 = 1671 4 (2.16) 
Thus 
eas oo a’ I? 
ee Dear (2.17) 
i=0 
OF adi Oe eutLg 
= ; 2.18 
Dd aii t 20 Gia os) 
= cosh (a) + sinh (a) I (2.19) 


since [7 = 1. 


In the Minkowski plane all vectors of the form a4 = a (Yo +71) are null (a = 0). One question 
to answer are there any vectors b+ such that az -b, = 0 that are not parallel to ax. 


by + by = 0 
at = +b 


Thus bs must be proportional to a4 and the are no vectors in the space that can be constructed 
that are normal to az. Thus there are no non-zero bivectors, a A b, such that (aA b)? = 0; 
Conversely, if aA b 4 0 then (aA b)* > 0. 


Finally for the condition that there always exist two orthogonal vectors e; and é2 such that 
aNb= ee (2.20) 


we can state that neither e; nor é2 can be null. 


l= 
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2.4 Lorentz Transformation 


We now have all the tools needed to derive the Lorentz transformation with Geometric Algebra. 
Consider a two dimensional time-like plane with with coordinates t? and x; and basis vectors Yo 
and 7;. Then a general space-time vector in the plane is given by 


r=tywtumn=tytri%i (2.21) 


where the basis vectors of the two coordinate systems are related by 


7, = Ry, R! w=0,1 (292) 
and R is a Minkowski plane rotor 
R=sinh (=) + cosh (=) 10 (2.23) 
so that 
RyR' = cosh (a) yo + sinh (a) 7 (2.24) 
and 
Ry,R' = cosh (a) 7, + sinh (a) yo (2.25) 


Now consider the special case that the primed coordinate system is moving with velocity § in 
the direction of y, and the two coordinate systems were coincident at time t = 0. Then x, = (t 
and x', = 0 so we may write 


tyo + Bty, = t'RyoR' (2.26) 
} (90 + 871) = cosh (a) yo + sinh (a) Y1 (2.27) 
Equating components gives 
cosh (a) = , (2.28) 
shia) = “8 (2.29) 


t 
Solving for a and — in equations 2.28 and 2.29 gives 
t! 


tanh (a) = 6 (2.30) 


?We let the speed of light c = 1. 


2.4. LORENTZ TRANSFORMATION 


t 1 
~~ FR 


Now consider the general case of x,t and 2’, t’ giving 


ty tay, =U RyR'4+ 2'Ry,R! 
= t'y (yo + Bu) + 2'y (11 + 80) 


Equating basis vector coefficients recovers the Lorentz transformation 


(2.34) 


Chapter 3 


Geometric Calculus - The Derivative 


3.1 Definitions 


If F (a) if a multivector valued function of the vector a, and a and b are any vectors in the space 
then the derivative of F’ is defined by 


F (a+ eb) — F(a) 


(3.1) 


b- VF =lim 
«>0 


then letting b = e, be the components of a coordinate frame with x = x*e, (we are using the 


summation convention that the same upper and lower indices are summed over 1 to N) we have 


F (xe; + ex) — F (xe;) 


e,:- VF = lim : (3.2) 
Using what we know about coordinates gives 
OF ; 
Bagh = eer. 
VF=e aoe O;F (3.3) 
or looking at V as a symbolic operator we may write 
V =e!0; (3.4) 


Due to the properties of coordinate frame expansions VF’ is independent of the choice of the e; 
frame. If we consider x to be a position vector then F(x) is in general a multivector field. 


33 
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3.2 Derivatives of Scalar Functions 


If f (x) is scalar valued function of the vector x then the derivative is 
Vi =f (3.5) 


which is the standard definition of the gradient of a scalar function (remember that in an or- 
thonormal coordinate system e, = e*). Using equation 3.5 we can show the following results for 
the gradient of some specific scalar functions 


i RAD Oe, Hee 


Wie: sig: se OF 2) 


3.3 Product Rule 


Let o represent a bilinear product operator such as the geometric, inner, or outer product and 
note that for the multivector fields F and G we have 


On (Fo G) = (OF) oG+ Fo (O,G) (3.7) 
so that 
V (FoG) = e* ((O:F)oG+ F 0 (&G)) 
=e (0,F)°oG+e*F 0 (OG) (3.8) 
However since the geometric product is not communicative, in general 
V(FoG) 4(VF)oG+ Fo (VG) (3.9) 
The notation adopted by Hestenes is 
V(FoG)=VFoG+VFoG (3.10) 
The convention of the overdot notation is 
i. In the absence of brackets, V acts on the object to its immediate right 


it. When the V is followed by brackets, the derivative acts on all the the terms in the brackets. 


it. When the V acts on a multivector to which it is not adjacent, we use overdots to describe 
the scope. 


Note that with the overdot notation the expression AV makes sense! 
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3.4 Interior and Exterior Derivative 


The interior and exterior derivatives of an r-grade multivector field are simply defined as (don’t 
forget the summation convention) 


V - A, = (VAr),_1 = €F- RAP (3.11) 


and 
VA A, = (VA;),., =e A OA, (3:19) 


Note that 


VA(VA A,r) = €'6; (e? A Oj Ar) 
= e’ A e/ ix (0,0; A,) 
= 0 (3.13) 


since e’ A e/ = —e! Ae’, but 0,0; A, = 0,0;A,. 


V-(V-A,) =e! +0; (e! . 0;A,) 
=e'- (e!- (0,0;A,)) 
—+e’. (e ie (0,0; A*I )) 
= te’ - ((e’ A (0,0;A7)) I) 
= +(e’ A (e’ A (0,0;A5))) I 
=a) (3.14) 


Where * indicates the dual of a multivector, A* = AI (J is the pseudoscalar and A = +A*IJ since 
I? = +1), and we use equation 1.53 to exchange the inner and outer products. 


Thus for the general multivector field A (built from sums of A,’s) we have V A (V A A) = 0 and 
V-(V-A) =0. If ¢ is a scalar function we also have 


V A (V6) =e! AG, (e404) 
= e /\ e/0;0; 
= (3.15) 


Another use for the overdot notation would in the case where f (x, qa) is a linear function of its 
second argument (f(z,aa+ 8b) =af(x,a)+Ef(x,b)) and a is a general function of position 
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(a(x) = a'(x) e;). Now calculate 


Defining 
= ee) O 
Vif (a) =a aa (a, e€;) = ef pet (24) 


a=constant 
Then suppressing the explicit « dependence of f we get 


The basic identities for the case of a scalar field a and multivector field F’ are 


V(aF) =(Va)F+aVF 
V -(aF) =(Va)-F+aV-F 
VA(aF) =(Va)AF+aVAF 
if f; and fo are vector fields 
VACA LR) HVA) AR-WVA RA 
and finally if F. is a grade r multivector field 
VA BOS(V ARSE 


where I is the psuedoscalar for the geometric algebra. 


(3.16) 
(3.17) 
(3.18) 


(3.19) 


(3.20) 
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3.5 Derivative of a Multivector Function 


For a vector space of dimension N spanned by the vectors wu; the coordinates of a vector x are the 
vc =2-u' so that c= v'u; (summation convention is from 1 to N). Curvilinear coordinates for 
that space are generated by a one to one invertible differentiable mapping from (ae age ) <> 
(0',...,0%) where the 6’ are called the curvilinear coordinates. If the mapping is given by 
«(0',...,0%) = 2'(6',...,0%) u; then the basis vectors associated with the transformation are 
given by 

Ox = Ox" 
aoe age 
The one critical relationship that is required to express the geometric derivative in curvilinear 
coordinated is 


e,. = 


(3.30) 


OOF 
k a 
i Agi (3.31) 
The proof is 
Ox OOF re 
Ox OO" 
Ox™ OOK 
~ BF Bam ew 
Oo" 


We wish to express the geometric derivative of an R-grade multivector field Fg in terms of the 
curvilinear coordinates. Thus 


VFrR=u (3.36) 


OF _ (00°) OFr _ .0Fn 
dx'  \ dxi) Ok ~~ OOF 


Note that if we start by defining the e;’s the reciprocal frame vectors e* can be calculated 


geometrically (we do not need the inverse partial derivatives). Now define a new blade symbol 
by 
= €j, Ryeseds Cin (3.37) 


F= Putney, ipl (3.38) 


yong 
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Then — 
QPin~in ee 
VE = —pr eel ee ir] 4+ FY Reh apr ak ip] (3.39) 
Define 9 
C {efi, Rasy et = ep eli ae tR] (3.40) 


Where C { eri, 


we can write 


grees 


OF%--tR ‘ fase 
VF —_ por © Eli, Ma ip] + F = RC { efi, ere ip] i} (3.41) 
Note that all the quantities in the equation not dependent upon the F“*® can be directly 
calculated if the e; (0', rey ) is known so further simplification is not needed. 


In general the e;’s we have defined are not normalized so define 


lex| = v/le2| (3.42) 


ex 


é, = — 3.43 
ey (3.43) 
and note that 6? = +1 depending upon the metric. Note also that 
é* = |e,| e* (3.44) 
since 
' . e le, : 
é) . &, = (|e,|e’) - (4) - i eal = § (3.45) 
lex.| ex.| 
so that if F'r is represented in terms of the normalized basis vectors we have 
PR = Bere: Pee iz) (3.46) 
and the geometric derivative is now 
OF. wR e* r P a ‘ 
vi 06k Jez) ne jae aC { ey, al ial} (3.47) 
and : 
Rahs es &, 
C {éfi, me inl} - Je,| oF re ir] (3.48) 
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3.5.1 Spherical Coordinates 


For spherical coordinates the coordinate generating function is: 


so that 


where 


and 


x =r (cos (0) uz + sin (6) (cos (d) uz + sin (?) uy)) 


e, = cos (0) (cos (¢) uz + sin (¢) uy) + sin (6) wu, 
eg = r (— sin (A) (cos (@) uz + sin (¢) wy) + cos (A) uz) 
€4 = 7 cos (8) (— sin (¢) wu, + cos (¢) uy) 


le,|=1 l|eg]=r leg] =rsin (6) 


é,, = cos (8) (cos (¢) uz + sin (¢) uy) + sin (8) uw, 


€y = —sin (0) (cos (@) uz + sin (d) uy) + cos (6) uw, 
é4 = —sin(¢) uz + cos (d) Uy 


the connection mulitvectors for the normalize basis vectors are 


eae 
Ce 
Gi COSI) Tet 9 5 
CN aa Bh poe 
oa at cos(™) . 
PENG CRINCASE Sp 9 \ €¢ 
A iad “p. ualae > HOR CE) il 
Cer eine: r sin (0) igre 
A 7 1. cos(@) . 7 
CLéeA Cg} = ~&¢ — ane \ €9 A Eg 
2 
C {ég A Eg} = —E, A Eg A Eg 
r 


C {é, \ &g A és} = 0 
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(3.49) 


(3.50) 
(3.51) 
(3.52) 


(3.53) 


(3.54) 
(3.55) 
(3.56) 


(3.57) 


(3.58) 
(3.59) 
(3.60) 
(3.61) 


(3.62) 
(3.63) 
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For a vector function A using equation 3.41 and that VA=V-A+VAA 


“A= ne Q b i r @ ie 
MSA) Oe eee ae (3.64) 
1 F 1 
7 ao Cae r sin (0) (Ap (sin (0) A®) + 35A*) (3.65) 
NE) (3.66) 
= Op A? 1 $ > 4 
= ( eng er (3.67) 
i“ ? 
( O6A _ A = 0,°) Ey (3.68) 
rsin(@) or 
A’ 9 OA" , 
(= TUR ge ) és (3.69) 
1 : . 
Ve ~ sin (0) (2 (sin (9) A®) - 3A’) er (3.70) 
1 1 
Lo ro Pn ee 
r \sin Oka Or (rA )) 0 (3.71) 
il 
mal: ra (0, (r.A’) = Op A") Ey (3.72) 


These are the standard formulas for div and curl in spherical coordinates. 


3.6 Analytic Functions 


Starting with G (2,0) and orthonormal basis vectors e, and e, so that J = e,e, and 1? = —1. 
Then we have 
P= Ge, yey (3.73) 
0 0 
Vie, = 3.74 
e Ox 1 ey Oy ( ) 


Map r onto the complex number z via 


z=ax2+Iy=e,r (3.75) 


3.6. ANALYTIC FUNCTIONS Al 


Define the multivector field ~ = u+ Iv where u and v are scalar fields. Then 


Ou Ov Ov Ou 


Thus the statement that ~ is an analytic function is equivalent to 


Vv =0 (3.77) 


This is the fundamental equation that can be generalized to higher dimensions remembering that 
in general that w is a multivector rather than a scalar function! To complete the connection with 
complex analysis we define (zt = x — Ty) 


O- —lf-0 O O Lf oO O 
= —I = + I 3.78 
Ox 2 (zs =| * -O2t- 2 (5 =| S78) 
so that 
Oz Ozh 
dz? Oz 
Bz ; Sct (3.79) 
On? Oat 
An analytic function is one that depends on z alone so that we can write wv (x + Iy) = w (z) and 
Ow (2) 
=0 3.80 
Fat (3.80) 
equivalently 
1/0 O 1 
Ninel fi = —e, Vv =0 3.81 
2 (= a =) " 9° ve ee) 
Now it is simple to show why solutions to Vw = 0 can be written as a power series in z. First 
Vz=V (er) 
- Or fads Or 
x aye y "Oy 


= €;€, Cy + Cy€z ly 
= €y — Cy 


= (3.82) 


so that 
V (z— 2)" =kV (esr — 2) (z — 2)" | =0 (3.83) 


Chapter 4 


Geometric Calculus - Integration 


4.1 Line Integrals 


If F(x) is a multivector field and x (A) is a parametric representation of a vector path (curve) 
then the line Integral of F' along the path « is defined to be 


[Fo a = [ru = dim )) FAs! (4.1) 
i=l 


where 


Az’ = Uj VM-1, Ft = (F (4-1) + F (x5)) (4.2) 


1 
2 
Ces 5 Ag 

if x, = x, the path is closed. Since dz is a vector, that is F (x) D x De (x), a more general 


line integral would be 


dX 


The most general form of line integral would be 


: ie 6G a= / F (x) dx G (a) (4.3) 


/ L(dyx;x) dd = / L (dx) (4.4) 


where L(a) = L(a;x) = is a multivector-valued linear function of a. The position dependence in 
L can often be suppressed to streamline the notation. 
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4.2 Surface Integrals 


The next step is a directed surface integral. Let F' (a2) be a multivector field and let a surface be 
parametrized by two coordinates x (x', x*). Then we can define a directed surface measure by 


dX = or A or dx'dx? = e, A ey drdx? (4.5) 


A directed surface integral takes the form 
fF dX = [Fe A € dx‘ dx? (4.6) 


In order to construct some of the more important proof it is necessary to express the surface 
integral as the limit of a sum. This requires the concept of a triangulated surface as shown Each 


ie 


Figure 4.1: Triangulated Surface 


triangle in the surface is described by a planar simplex as shown The three vertices of the planar 
simplex are Xo, 41, and x2 with the vectors e; and e»2 defined by 


€, = %—-T%o, €2 = %2—TXo (4.7) 


so that the surface measure of the simplex is 


1 1 
AX = 561 Neg = 5 (t1 Aka + 2A to + Xo A 21) (4.8) 
with this definition of AX we have 
= Ak k 
prex = Tim as F°AX (4.9) 


where F* is the average of F over the k'” simplex. 
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L2 


e€2 


Xo e1 Ly 


Figure 4.2: Planar Simplex 


4.3 Directed Integration - n-dimensional Surfaces 


4.3.1 k-Simplex Definition 


In geometry, a simplex or k-simplex is an k-dimensional analogue of a triangle. Specifically, a 
simplex is the convex hull of a set of (k +1) affinely independent points in some Euclidean space 
of dimension k or higher. 


For example, a 0-simplex is a point, a 1-simplex is a line segment, a 2-simplex is a triangle, a 3- 
simplex is a tetrahedron, and a 4-simplex is a pentachoron (in each case including the interior). 
A regular simplex is a simplex that is also a regular polytope. A regular k-simplex may be 
constructed from a regular (& — 1)-simplex by connecting a new vertex to all original vertices by 
the common edge length. 


4.3.2 k-Chain Definition (Algebraic Topology) 


A finite set of k-simplexes embedded in an open subset of R” is called an affine k-chain. The 
simplexes in a chain need not be unique, they may occur with multiplicity. Rather than using 
standard set notation to denote an affine chain, the standard practice to use plus signs to separate 
each member in the set. If some of the simplexes have the opposite orientation, these are prefixed 
by a minus sign. If some of the simplexes occur in the set more than once, these are prefixed with 
an integer count. Thus, an affine chain takes the symbolic form of a sum with integer coefficients. 
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4.3.3. Simplex Notation 


If (vo, 21,..., 2%) is the k-simplex defined by the k + 1 points x9, 71,..., 2%. This is abbreviated 
by 
(X) (4) = (Xo, %1,---, Xk) (4.10) 


The order of the points is important for a simplex, since it specifies the orientation of the 
simplex. If any two adjacent points are swapped the simplex orientation changes sign. The 
boundary operator for the simplex is denoted by 0 and defined by 


k 
A(x) yy = d_ (1) (ao,-- By - Be) G—1) (4.11) 
i=0 
To see that this make sense consider a triangle ()(3) = (%o, 1,2). Then 


O(X) (3) = (21, £2),2) = (Zo, £2) 9) 7 (0, £1) (2) 
= (21, £2),9) = (2, Lo) ,2) Tr (Zo, 21) 2) (4.12) 


each 2-simplex in the boundary 2-chain connects head to tail with the same sign. 
Now consider the boundary of the boundary 


om (X) (3) = O (@1, La) a) + O (La, Lo) + O (Zo, 21) 2) 
= (21) a) — (#2) + 2) ay — @o) ay + (0) ay — (1) a) 


2} (4.13) 
We need to prove is that in general 0? (x) (x) = 9. To do this consider the boundary of the ye 
term on thr r.h.s. of equation 4.11 letting Ay SNP scioe lane » Zk) (1): 
Then ; 
O (xo, Sa Xi, Ror » Zk) (4-1) = 
k 
—1 4(k—-2 
i=0 So (=p ae 
j=l 
i-1 k (4 14) 
O<i<k: So(-1f ah? SO (ryt ae 
j=0 jai41 
k-1 . 
i=k (apr Ae 
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The critical point in equation 4.14 is that the exponent of —1 in the second term on the r.h:s. 
is not 7, but 7 — 1. The reason for this is that when x; was removed from the simplex the 
vertices were not renumbered. We can now express the boundary of the boundary in terms of 
the following matrix elements (Be?) = (-1)'4 Ae) as 


aj 


k k-1 
_ k-2 k 1 k-2 
O (x) ay => (-1)7 7 AN + (- 1) SS (1) A 
j=l j=0 
k-1 i-l k 
4 j k-2 _ k-2 
So (S's 4 (aya } 
i=1 j=0 j=i4+l 
k k-1 
k-2 k-2 
= ayy ay 
j=l j=0 
k-1 i-1 k-1_ k 
+ B= 5 Be a0 (4.15) 
i=1 j=0 i=1 j=i+l 


The consider BY ~) as a matrix (i-row index, j-column index). The matrix is symmetrical and 
in equation 4.15 you are subtracting all the elements above the main diagonal from the elements 
below the main diagonal so that 0? (Z) (x) = 0 and the boundary of a boundary of a simplex is 
Zero. 


Now add geometry to the simplex by defining the vectors 


e€; = @i— Xo, a TE A (4.16) 
and the directed volume element 
1 
AX = park Nee (4.17) 
We now wish to prove that 
, dX = AX (4.18) 
() (x) 


Any point in the simplex can be written in terms of the coordinates ' as 


k 
t= Boob Ss Ne; (4.19) 
i=l 
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with restrictions 


k 
C5 NST and SN <I (4.20) 
i=1 
First we show that 
ax= [ er A+++ Ae, dd! -+-dd*® = AX (4.21) 
(©) (x) (©) (xe) 
or 
1 
/ dx sar (4.22) 
kl 
(®)(%) ; 


define Aj = 1 — S~J_, \' (Note that Ag = 1). From the restrictions on the ’’s we have 


No A Ag-1 
if dy ...dyt = i dn! | ee | dN (4.23) 
(x) 0 0 0 


To prove that the r-h.s of equation 4.23 is 1/k! we form the following sequence and use induction 
to prove that V; is the result of the first 7 partial Integrations of equation 4.23 


(k) 


(Neg) (4.24) 


Then 


=i p—jyiti] Ato 
~ G+))j! (Anat 7a’) i 
I j+1 
— Gai (Ap_j-1) (4.25) 


so that V, = 1/k! and the assertion is proved. Now let there be a multivector field F(x) that 
assumes the values F; = F' (x;) at the vertices of the simplex and define the interpolating function 


k 


f(t)=Fot+ > NA - A) (4.26) 


i=1 


4.3. DIRECTED INTEGRATION - N-DIMENSIONAL SURFACES 49 


We now wish to show that 


k 
1 _ 
dX = —— y Fe) AX = FAX 4.27 
he t k+1 ( ; ( ) 


To prove this we must show that 
: 1 : 
i \' dX = ——AX, VX (4.28) 
(©) (n) k+1 


To do this consider the integral (equation 4.25 with V; replaced by \*~/V;) 


An—j-1 ey Ag—j-a ms eee oe 
i ary NIV; = i dnt 4 (Ay—j—1 — Ne) 
0 0 J: 


~G : 2)! (An—j—1)" 29) 


Note that since the extra \’ factor occurs in exactly one of the subintegrals for each different. X’ 


the final result of the total integral is multiplied by a factor of =o) since the weight of the total 


integral is now and the assertion (equation 4.28 and hence equation 4.27) is proved. 


ol cae 
(K+)! 


Now summing over all the simplices making up the directed volume gives 


FdX = lim YS F*AX' (4.30) 
i Dee > 


The most general statement of equation 4.30 is 


| L(dX) = lim S7L' (AX’) (4.31) 
volume Ree a 


where L(F,,;) is a position dependent linear function of a grade-n multivector F,, and L' is the 
average value of L(dX) over the vertices of each simplex. 


An example of this would be 
Lee) = Ge) Fie (@) (4.32) 


where G (x) and H (x) are multivector functions of x. 
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4.4 Fundamental Theorem of Geometric Calculus 


Now prove that the directed measure of a simplex boundary is zero 


A (a («),)) =iK (a (Binass 2.) ay) =0 (4.33) 
Start with a planar simplex of three points 
0 (Xo, £1, 2) (2) — (21, £2),1) > (Xo, £2) (1) Ez (20,21) (1) (4.34) 


so that 
A (8 (#0, 1,£2),a) = (#2 — #1) — (#2 — #0) + (#1 — 20) = 0 (4.35) 


We shall now prove equation 4.33 via induction. First note that 


1 
08 ——A (Zo) (,_9) A (@k — £1) 
A (8) 1 = an: (4.36) 
and 
; 1 
A (£k) @-1) = (k—1) (21 — Zo) A+++ A (@p-1 — Xo) (4.37) 
so that 
1 ko . 
4 (2a) = Fay A Ban A Ge 20) +6 
where 
1 
C= ko15 (0) (n—2) A (@e — 21) + (-1)"A (Zk) (4-1) (4.38) 
if we let 6 = 2% — X1 we can write 
1 3 ¥ y 
C=7—A (Zo) (na) A (%e — Zo) + 7214 (Zo) q_2) AO + (-1)"A (2x) q_ay (4.39) 
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Then 
1 
A (Xo) 4-2) Ad= (k = DI (xq — 21) Pee AN (Le-1 — £1) AO 
1 
= aoe oO ee ao) oe 
1 
Ene 
= (pay ee Eee 0) 
ae 
= papi (1 ~ 0) A (Be — 0) As A (ea — 20) 
Thus 
=| ; 
Ea) olga AO= 
(-1)" 
= Opa yl (1 20) A (2 — 0) Aor A (ha ~ 0) 
=(=1) A (2k) (4-1) 
However , 
(=1)* A (2) (41) ae cS {4 (Zo) (42) A6 
so that ; 
a ne {4 (Zo) (x2) A (ke — &o) 
and 
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(4.40) 
(4.41) 


(4.42) 
(4.43) 


(4.44) 


(4.45) 


(4.46) 


(4.47) 


(4.48) 


(4.49) 


(4.50) 


We have proved equation 4.33. Note that to reduce equation 4.49 we had to use that for any set 


of vectors 0, Y1,-.., Ye we have 


SAY +O) As AYR +d) =OANA AYE 


(4.51) 
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Think about equation 4.51. It’s easy to prove (6 Ad = 0)! 
Equation 4.33 is sufficient to prove that the directed integral over the surface of simplex is zero 
k 4 
¢ dS = 5 ~(-1)' px =A (a (2)ay) =0 (4.52) 
A(X) (x) 5=0 (#1) (1) 
The characteristics of a general volume are: 
1. A general volume is built up from a chain of simplices. 


2. Simplices in the chain are defined so that at any common boundary the directed areas of 
the bounding faces are equal and opposite. 


3. Surface integrals over two simplices cancel over their common face. 


4. The surface integral over the boundary of the volume can be replaced by the sum of the 
surface integrals over each simplex in the chain. 


If the boundary of the volume is closed we have 


gas = Jim fas =0 (4.53) 


Where § dS" is the surface Integral over the a’” simplex. Implicit in equation 4.53 is that the 
surface is orientated, simply connected, and closed. 


The next lemma to prove that if b is a constant vector on the simplex (X) (x) then 
§ b-cdS=b-A((c)) =b- AX (4.54) 
A(x) (1) 
The starting point of the lemma is equation 4.28. First define 
k 
b=) die’, (4.55) 
i=1 
where the e’’s are the reciprocal frame to e; = 2; — Xo so that 


k 
L— X= » re, (4.56) 
i=l 
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b, =b- e;, (4.57) 
and 
k 
S > Xb; = b+ (w@— a). (4.58) 
i=1 
Substituting into equation 4.28 we get 
k i 
bA dX = J b-(a— dX = —— b-e; AX 4.59 
» [oe [ee A) k a 1 » . 


Rearranging terms gives 


iI k 
[o-2ax=o5 ((d24: ta) + (b+ 1b AX 


©) (4) i=1 
5 k 
=—_.[S°2,| Ax 
a (S-) 
=b-ZAX (4.60) 


Now using the definition of a simplex boundary (equation 4.11) and equation 4.60 we get 


¢ b-xdS = 
A(X) (x) 


1 1 
The coefficient multiplying the r.h.s. of equation 4.61 is i and not fat because both 


(to + +++ +2; +--+ + 2%) and (%;)”,_1) refer to k—1 simplices (the boundary of (x)(,) is the sum 
of all the simplices (Xi) (4) with proper sign assigned). 


S(-'0- (a fee- +H; +---4 m) A (Han) (4.61) 


a] re 


Now to prove equation 4.54 we need to prove one final purely algebraic lemma 


C= (-1)'d- (xo Parga alesse ai tr) A (Bi) = DG (4.63) 
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where C; is defined by 


Ae eae b- (1 +-+-+2%) A (Xo) qa) 
an| 0<i<k: (-1)'d- (xp i epee el aan d xp) A (#8) ny (4.64) 


now define F;, = e; \--- Ae, so that (using equation 1.61 from the section on reciprocal frames) 


(—1)* e'Ey =e A-+-A&j A+ Ale, VO<i<k (4.65) 
and ‘ 
Co<ish = hl ‘(tot +B +-+- +24) Ey (4.66) 
The main problem is in evaluating Cp since 
: 1 
A (%0)(@—1) = k—D! (Zz — £1) A+++ A (&e — 21) (4.67) 
using €; = x; — % reduces equation 4.67 to 
1 
A (£0) (K1) = DI (€2 — €1) A-++ A (ex — e1) (4.68) 


but equation 4.68 can be expanded into equation 4.69. The critical point in doing the expansion 
is that in generating the sum on the r.h.s. of the first line of equation 4.69 all products containing 
x, (of course all terms in the sum contain x; exactly once since we are using the outer product) 
are put in normal order by bringing the x, factor to the front of the product thus requiring the 
factor of (—1)’ in each term in the sum. 


k 
(eg — €1) A+++ A (en — €1) = 2A AN @e— >. (=D er AeA NE Ao A ek 
i=2 


= SU (-1 ter neh AB A Ne 


k 
=S\ ek (4.69) 
or 


k 
3 1 Sie! 
i=1 
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from equation 4.69 we have 


on eae Cee 
ie oe 
= fay 2 64) e’ Ey 


k 
1 a 


1 
csi 
1 


= mop Be (4.71) 


and equation 4.62 is proved which means substituting equation 4.62 into equation 4.61 proves 
equation 4.54 


4.4.1 The Fundamental Theorem At Last! 


The simplicial coordinates, \’, can be expressed in terms of the position vector, x, and the frame 


vectors of the simplex, e; = 7; — 2%. Let the vectors e’ be the reciprocal frame to e; (e;-e? = 6!). 
Then 


XN = e'- (x — x0) (4.72) 


and let f(x) be an affine multivector function of x (equation 4.26) which interpolates, F’, a 
differentiable multivector function of 7 on the simplex. Then 


(x) dS 


k 

A(x) (1) 2 
k 

oy 


(FF) f (a0) aS 


) (4) 
(F, — Fo) (AX) (4.73) 


But 
= F,— Fy (4.74) 
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so that the surface integral of equation 4.73 can be rewritten 


$f (0) aS = YR Feet (AX) 
= a ore (XS IV (AR) (4.75) 


If we now sum equation 4.75 over a chain of simplices realizing that the interpolated function 
f (x) takes on the same value over the common boundary of two adjacent simplices, since f (x) 
is only defined by the values at the common vertices. In forming a sum over a chain, all of the 
internal faces cancel and only the surface integral over the boundary remains. Thus 


¢ f(a) dS = 5° fV- (AX?) (4.76) 


with the sum running over all simplices in the chain. Taking the limit as more points are added 
and each simplex is shrunk in size we obtain the first realization of the fundamental theorem of 


geometric calculus 
¢ FdS= h FV dX (4.77) 
av V 


We can write VdX instead of V-dX since the vector V is totally within the vector space defined 
by dX so that V A dX = 0. The method of proof used can also be applied to the form 


¢ is@=[ vax (4.78) 
OV V 


A more general statement of the theorem is as follows: 


Let L(Ag_1) = L(Ag-1;x) be a linear functional of a multivector A,_; of grade k — 1 and a 
general function of position « which returns a general multivector. The linear interpolation 
(approximation) of L over a simplex is defined by: 


L(A) =L(Aj20) + DW (L (Asai) — L(A; 20)) (4.79) 


Then the integral over a simplex is (Note that since integration is a linear operation, a summation, 
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the integral can be placed inside L(A) since L(A) is linear in A) 


t. L (dS) = L(f asix0) +DoU(f vasin) ~ Sou (f xasia) 


= y (L(e'AX; a;) —L (AX; x0) 


(k) 


=i, (vax) (4.80) 


Taking the limit of a sum of simplicies gives 


f (as) - i L (vax) (4.81) 


4.5 Examples of the Fundamental Theorem 


4.5.1 Divergence and Green’s Theorems 


As a specific example consider L(A) = (JAI~') where J is a vector, I is the unit pseudoscalar 
for a n-dimensional vector space, and A is a multivector of grade n—1. Then equation 4.81 gives 


[ (avaxr) = fv slaxi= g (ast) (4.82) 


we have dX = I |dX| where |dX| is the scalar measure of the volume. The normal to the surface, 
n, is defined by 


n|dS| = dSI~* (4.83) 
where |dS| is the scalar valued measure over the surface. With this definition we get 
[v-siax| -¢ ne FlaS (4.84) 
V av 


Now using the form of the Fundamental Theorem of Geometric Calculus in equation 4.78 and let 
G be the vector J in two-dimensional Euclidean space and noting that since dA is a pseudoscalar 
(for 2-D ds is the boundary measure and dA is the volume measure) it anticommutes with vectors 


in two dimensions we get - 
¢ dsJ = ‘i VdAJ = — / VJdA (4.85) 
aA A A 
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In 2-D Cartesian coordinates dA = Idxdy and ds = nI |ds| so that 


¢ dsJ = -{ VJ Idixdy. 
aA A 


¢ nIJ \ds =- | vatavdy 

aA A 

-$ nJI \ds == | vitaviy 
aA A 


¢ nJ |ds = | vaardy 
aA A 


Letting J = Pe, + Qe, and n= n*e, + n¥e, we get - 
nJ =n-J+(n"Q —nP) ene, 


or 


so that 


¢ n-F\ds|= f V- Saedy 
aA A 


0Q OP 
—Q — ny = pts oe cael 
fe (n*Q — n¥P) |ds| e,ey i (F _ dxrdye,e€, 


but dy = n* |ds| and dx = —n¥ |ds| so that 


¢ Pdi Od¢= i} (3 — ~) dxdy 
OA A \ Ox Oy 


which is Green’s theorem in the plane. 


(4.86) 


(4.87) 


(4.88) 


(4.89) 


(4.90) 


(4.91) 


(4.92) 


(4.93) 


4.5.2 Cauchy’s Integral Formula In Two Dimensions (Complex Plane) 


Consider a two dimensional euclidean space with vectors r = ve, + yey. Then the complex 


number z corresponding to r is 
Z=e,r =x2+yerzey =x2+yl 
zl=re,=x-—ye,e,=xr—yl 


wager ty =r 


(4.94) 
(4.95) 
(4.96) 
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Thus even grade multivectors correspond to complex numbers since J? = —1 and the reverse of 
z, z' corresponds to the conjugate of z. Even grade multivectors commute with dS, since dS is 
proportional to I. 


Let ~(r) be an even multivector function of r, then 
[ voas = fase = or wd) (4.97) 
but the complex number z is given by z = e,r and 
ex f dev = f va: = | e,VvdsS. (4.98) 


Thus if a function ~ is analytic, Vy = 0 and ¢ ydz = 0. Now note that (this will be proved for 
the N-dimensional case in the next example) 


Vea = 276 (r — a) (4.99) 
where a = e,a. Now let 
r—a 
y= Goa) a (e271) (4.100) 


so that 


ci TE) gg (4.101) 
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f acme, [v(Fatrest en)) as 


= e. | (2r0 (r —a)e,f (e.r) + Vf (err) aes) I \ds| 


zi —al 
= 2nIf(a)+ | e:Vf (e2r) ——G! |dS| 
|z —al 
=2rIf(a)+ | ecVf (ez) _ mu |dS| 
O O 1 
= 4.102 
antsia)+ f (5 +15) £@ ——zTles| (4.102) 
If Vf (z) =0 we have - 
free = 2p (a) (4.103) 
z-a 
which is the Cauchy integral formula. If Vf (z) is not zero we can write the more general relation 
SIP aS ¢ fap] CE sigs (4.104) 
z-a Oztz—a 
since 5 eee 5 
Dat 3 (7 15) (4.105) 
os a_1/a a 
= —> Al 
Oz 2 (= 1x] eee) 


4.5.3. Green’s Functions in N-dimensional Euclidean Spaces 
Let w be an even multivector function or let N be even so that ~ commutes with J. The analog 
of an analytic function in N-dimensions is Vw = 0. 


The Green’s function of the V operator is (Sy = 27%/?/I'(N/2) is the hyperarea of the unit 
radius sphere in N dimensions) 


G(2;y) = lim (4.107) 
ea ee (Ie —yl + c) 
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So that 
lim ViG (x;y) =d(x@—y). (4.108) 
ne 
To prove equation 4.108 we need to use 
V2(e-y)=N and Vz|e—y|!" =M(e—y)|a—yl" 
So that 
1 
V-G = 5-4 Ve(le- ul" +e) (e-1) + (le-yl" +6) Vee-w} 
oa — |x —y|% 1 
= = 
Sw (Ic -yl® +.) (ee +e) 
N Sek ie 
=5 : >= Vr G=G-V, =GVe (4.109) 
su (Ie —y\v + c) 
so what must be proved is that 
N 
lim : 5 =6(e-y) (4.110) 
Span (Ix - yl +6) 
First define the volume B, (7 > 0) by 
LEB, Sle) <7 (4.111) 


and let y = 0 and calculate (Note that we use |dV]| since in our notation dV = I |dV| and the 
oriented dV is not needed in this proof. Also r = |z].) 


N ee) N N-1 
: slav|= | dr 
Bs; Sn (I2i" +e) 0 (r TT €) 


“hore 


=| (4.112) 
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Thus the first requirement of a delta function is fulfilled. Now let ¢(2) be a scalar test function 
on the N-dimensional space and S a point set in the space and define the functions 


max (¢, S) = {max (¢(xz)) Vz € S} and min(¢,S) = {min (¢(x)) Va € S} 


and calculate the integral 


= € i 
[rN +e], 
€ 
1- 4.11 
TN +€ (ry) 
and note that 
N € € 
— ay |= 4.114 
Sn Reo Bs (i +e) | TN +e ( ) 
Thus Vr > 0 we have 
N 
tim f —_—*___ |dv| =0 (4.115) 
e0 Sw Bo Br (I2I* +e) 
and 
N 
lim 3 | *___ldv|=1 (4.116) 
E> N a (Ie1* +e) 
Thus 
N 
min ($, Bo — B-) — “__|dv| < 
ON J B..—Br (It +e) 
N €—(x 
cca ( ) 2 |dV | 
N J Boo—Br (Iel* +«) 
N 
< max (¢, B., — B,) — 7 (4.117) 
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and 
N 
min (@, B,) 3 —, |av| < 
N JB, (Il + c) 
N 0) ay 
2 
Sw Jz, (Ie + c) 
N € 
< max (¢, B,) ele ( = y |dV | (4.118) 
t (la|" +€ 
Thus 
lim =| oe |dvV| =0 (4.119) 
«>0 N Bo Br (Ie le c) 
and 
mine) eine COE) ay ena (4.120) 
«0 S N 2 
Nat Be (Ie + c) 
Finally 
lim lim al i: a, ; |dV| = 4(0) (4.121) 
T0€> N JB, (I21 +e) 
and we have proved equation 4.110 since 
him (max(¢, B,) — min(¢, B,)) =0 (4.122) 


Now in the fundamental theorem of geometric calculus let L(A) = GAw so that (remember that 
wv commutes with dV and V,w = 0) 
= (Gv Ea GV.) T|dV| 


= GV WI |dV| 
= V,GyI |dV| (4.123) 
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and 


¢ Gasp = i ie =) W(o) lav] 
OV Vv 


= vy) I 
= Toty) (4.124) 
or - ss a 
Wo)= sh, prove (4.125) 


because 7 is a monogenic function (Vy = 0). 


Chapter 5 


Geometric Calculus on Manifolds 


5.1 Definition of a Vector Manifold 


The definition of a manifold that we will use is - 


A vector manifold (generalized surface) is a set of points labeled by vectors lying in a geometric algebra of arbitre 
dimension and signature. If we consider a path in the surface x (A), the tangent vector is defined by 


» _ Ox (A) 


ae. CIAG RE) = BA) 
Dr = lim 


e>0 € 


8 
II 


Ao 
and the path length 


2 
ssf VJ |x! x'|dr (5 
M1 


5.1.1 The Pseudoscalar of the Manifold 


Now introduce a set of paths in the surface all passing through the same point x. These paths 
define a set of tangent vectors {e1,...,€n}. We assume these paths have been picked so that the 
vectors are independent and form a basis for the tangent space at point x. The outer product of 


65 
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the tangent vectors form the pseudoscalar, J (x), for the tangent space 
Nex KA NRE, 
I(r)= ea 2 
le: Neg A+++ Ae, 


(5.3) 


Thus 

| eel (5.4) 
We require that for any point on the manifold the denominator of equation 5.3 is nonzero. We 
also assume that for all points on the manifold that J (az) is continuous, differentiable, singled 
valued, and has the same grade everywhere. 


5.1.2 The Projection Operator 


Define the projection operator P (A) operating on any multivector A in the embedding multi- 
vector space as 

P (A) =(A-I(x))I7 (2) (5.5) 
We can show that P (A) extracts those components of A that lie in the geometric algebra defined 
by I (a). Since P (A) is linear in A if we show that if P(A,) projects correctly for an r-grade 
multivector it will do so for a general mulitvector. If n is the dimension of the tangent space and 
A, is a pure grade multivector we can write equation 5.5 as 


P (Ar) = (Ard (2))),-n J" (2) (5.6) 


Now consider the blades that make up the components of A,. They will either consist only of 
blades formed from the tangent vectors e; or they will contain at least one basis vector that is 
not a tangent vector. In the first case 


(ArT (®)) pny = Art (x) (5.7) 

and 

P (A,) = Ar (5.8) 
In the second case there is no component of A,J (x) of grade |r — n| and 

P(A, =0 (5.9) 
This is easily seen if one constructs at point x an orthogonal basis (0; +0; = 6;;07) for the 
tangent space {0,,...,0,} and an orthogonal basis for the remainder of the embedding space 
{Ont+1;+-+;0m}. Then any component blade of A, is of the form 


Ave PN, OG Op. 12+ Oj (5.10) 


Tr 
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where 7) <tg <«*+ <4. Ifa; SnV1l<7 <r then 


Ait aa! © 0; Oig +++ O4,° l= Ait ay 0 Op, woe 0;,1 (5.11) 
(Ap al 04, 0ip 21+ O07, ° I) i = Ane ™ "Oe Oy, - +. 0%, (5.12) 
P (AC@ oo ” On Ogg x2 0i,.) Sh OO OE (5.13) 
If any i; > m then A'r9;,0;, ...0;,1 contains no grade |r — n| and 
P (An ae "04, Oig autis 0;,.) =0 (5.14) 
5.1.3. The Exclusion Operator 
For an arbitrary multivector A the exclusion operator P, (A) is defined by 
P, (A) =A-P(A) (5.15) 


5.1.4 The Intrinsic Derivative 


Given a set of tangent vectors {e;} spanning the tangent space and the geometric derivative, V, 
for the embedding space, the derivative intrinsic to the manifold is defined everywhere by 


d=e'e,-V=P(V) (5.16) 


Also note that 
P(0O)=90 (5.17) 


When we write P(V) or P (0) the V or O is not differentiating the J (x) in the P operator 
anymore than V is differentiating dX in the fundamental theorem of Geometric Calculus. 


We also note that if the vector a is in the tangent space that 
a-0=a-V (5.18) 


and that a-O is a scalar operator that gives the directional derivative in the a direction. Also since 
it is scalar it satisfies Leibniz’s rules without using the dot notation (remember the convention 
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that if parenthesis are not present the operator precedence is dot product then wedge 
product then geometric product). 

a:-O(AB) = (a-0A)B+A(a-OB) (5.19) 


An alternative definition for the intrinsic derivative! is to let y(s) be a curve on the manifold. 


d 
Then — is a tangent vector to the mainfold and we can define 


ds 
dy dA (7 (s)) 
—| -O = 5:20 
7 ede ds lave (5.20) 
We can show that equation 5.20 is equivalent to equation 5.16 with the following construction. 
Let (z!,...,2”) be a local coordinate system for the vector manifold x = x(z',...,2”). Then a 
a) 
basis for the tangent space is e; = aA and the intrinsic derivative is 0 = Caaf Now write 
x" x 
A(y(s)) = A(z" (s),...,2" (s)) (6:21) 
so that dA(y(s)) OA dx’ 
y(s oc 
— 0.22 
ds Ox’ ds ( ) 
sa d Ox dx’ dx* 
ry x dx’ x 
= ; =e; 0.23 
ds dsids ‘ds ( ) 
so that 


dy dt’ ,OA dx’ 0A dA(7(s)) 
—! BA Se [i —_ 

ds @ “ds © Oxi ds Oat ds 
and the two definitions are equivalent. 


(5.24) 


5.1.5 The Covariant Derivative 


The O operator is entirely within the tangent space and if the general multivector function A (x) 
is also entirely within the tangent space, it is still possible (even likely) that OA is not entirely 
within the tangent space. We need a covariant derivative D that will result in a multivector 
entirely within the tangent space. This can be done by defining 


a-DA(x) =P (a-0A(zx)) (5.25) 


1Private conversation with Dr. Alan MacDonald. 
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so that 
a-OA=P(a-0A)4+P, (a:0A) =a-DA+P, (a: OA) (5.26) 
Again since a- D is a scalar operator we have 


a D(AB) = P(a:0(AB)) = (a- DA) B+ A(a- DB) (5.27) 


A component expansion of D is given in the usual way by (do not forget the summation conven- 
tion) 


D=e'e,:D (5.28) 

and 
DA, = e' (e;- DA,) = P (OA,) (5.29) 

and 
Ds A, = (DA) a4 (5.30) 
DRA SDA Le (5.31) 


if a (x) is a scalar function on the manifold then 
Oa (x) = Da (x) (5.32) 


because in equation 5.32 no basis vectors are differentiated. To relate O and D if the function 
operated on is not a scalar first construct a normalized basis {e;} for the tangent space at point 
x. Then 

T=e,ANe.A...Ae, and I? = +1 (5.95) 


and (since a-0 and a- D are scalar operators we can move them across the wedge products 
without any problems) 


(acon T= (soe A...A(a- De, + Py (a - 0e;)) Ano Aen) ie (5.34) 
i=1 
=(a-DI)I4*4+5_ (-1)°" Py (a: de) AeA... NGA... A el (5.35) 
1=1 
=(a-DI)I-' +P (a- de) Ae (5.36) 


We go from equation 5.35 to equation 5.36 by using equation 1.61 on page 33 in the section on 
reciprocal frames. 
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Since (a- D)I is a grade n multivector in the tangent space it must be proportional to J and 
thus commute with I so that ((a-0)I)I =I ((a-0) I). Also I~’ = +I so that we have 


(@-DDI =A (a2DI)s 


= 45 ((a-DI)I +I (a- DI) 
= +5 (a-D(P)) 
_9 (5.87) 
Thus 
(a- OI) = P, (a- 0e;) Ne’ I = —S(a)I (5.38) 


Where S'(a) is the shape tensor associated with the manifold. Since S (a) is a bivector we can 
write (A x B = (AB — BA) /2) 

a-0I =I x S(a) (5.39) 
since 


S(a)-l=S(a)AT=0 (5.40) 


and by equation 1.96 page 58. Given that a(x) and b(x) are vector fields on the manifold (both 
are in the tangent space at point x), form the expressions in equation 5.41 (remember that for 
any three vectors u, v, and w we have u- (vA w) = (u-v)w— (u-w)v) 


b- S(a) =b- (e' AP (a- 0e;)) 
= (b-e') P, (a: de;) — (b- P, (a- de;)) e* 
= ((b’e;) -e') Pi (a: dei) 
= WP. (a - 0e;) 


Py (a ve) (5.41) 
but 
P, (a- db) =P, (a- 0 (b'e;)) 
=P, (a Abie; +a- ave) 
=P, (a- dvié;) (5.42) 
and 


b- S(a) =P, (a- 0b) (5.43) 
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Thus 
a: 0b=P(a- 0b) + Py (a: 0b) =a-Db+b- S(a) (5.44) 


and (using the fact that the dot product of a vector and bivector are antisymmetric) 
a-Db=a-0b+S(a)-b (5.45) 


Now consider the expression 


a-D(b,...bp) = 5° by... (a+ 8B; + 5 (a) - Bj)... b, 
1=1 


i=1 


ey ee 
som AS). 
= 0-9 (br..br) +5 (S (a) (bs bj bred Sta 4 
(3 S (a) b; bh b;S (a) nya 
: (by . b,p1S) (a) b,. = 6S (a) bo Perr b,) 
= 0-9 (br -.-Br) +5 (S (a) (br. Br) ~ (br. by) S (a) + 
=a-0(b,...b,) +S (a) x (by...5,) (5.47) 


To get from equation 5.46 to equation 5.47 note that in the sums in parenthesis in equation 5.46 
the i‘” term in the first sum cancels the i” + 1 term in the second sum. 


Since any multivector is a linear superposition of terms containing b,...b, with 1 <r <n and 
a scalar we have 


a-DA=a-0A+S(a)xA (5.48) 
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Where a(x) and b(z) are vector fields on the manifold write 
a-0b=a-OP(b) =a-OP (b) +P (a- 0b) =a- OP (b) +a- Db (5.49) 
Now substitute equation 5.48 into equation 5.49 to get 


a- OP (b) =b- S(a) (5.50) 
5.2 Coordinates and Derivatives 
In a region of the manifold we introduce local coordinates x’ and define the frame vectors as 


= Ox 
Oxi 


From the definition of 0 is follows that e’ = Ox’. The {e;} are referred to as tangent vectors and 
the reciprocal frame {e’} as the cotangent vector (or 1-forms). The covariant derivative along a 
coordinate vector, e;- D, satisfies and defines both D; and Sj. 


The tangent frame vectors satisfy 
Oje; — 0j;e; = (0,0; — 0;0;) x =0 (5.53) 
Using the P operator on equation 5.53 gives 
Dye; — Djex = 0 (5.54) 


while using P, gives 
ei: ay =€;° Si (5.55) 


For arbitrary vectors a and b in the tangent space equation 5.55 becomes 
a-S(b)=b-S(a) (5.56) 


In terms of the coordinate vectors the shape tensor becomes 


S (a) =e* AP, (a- Oex) (5.57) 
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and 
S;=e* NP, (e&;- dex) = e* A PL (ex « Oe;) (5.58) 
Then 
ONE, =e* A Ope; = e* A (P (Apex) + Pi (pei) = DAG + S; (5.59) 
Letting a = a’e; be a constant vector in the tangent space gives the general result 
O\a=DAa+S(a) (5.60) 
Additionally 
O\a=ON(P(a)) 
=OAP(a)+P(0Aa) 
=DAat+dAP(a) (5.61) 
Thus a 
OAP (a) = §(a) (0262) 


Note that if a and b are any two vectors in the embedding space then P (a A b) = P (a) A P (b) 
and if #(x) is a scalar function on the manifold we have 


O\N\0¢=OAP (VO) 
=dAP (Vo) +dAP (Vo) 
=P(V) AP (Vo) +0AP (V9) 


=P(VAV¢)+9AP (V4) (5.63) 
but V A V =0 so 
d\ Ab = S (V6) (5.64) 
Since S'(a) for any vector a lies outside the manifold we have 
DA(Dé¢) =0 (5.65) 


Letting d(x) = 2’ (x), then since x* (x) is a scalar function 
DA (Dz) = DAe' =0 (5.66) 
so that for a general vector a = a; (x) e’ we have 
DAa=DA (aje’) =e' Ae! (d;a;) = oe /\ e? (O;a; — 0;a;) (5.67) 


Equation 5.67 is isomorphic to the definition of the exterior derivative of differential geometry. 


74 CHAPTER 5. GEOMETRIC CALCULUS ON MANIFOLDS 
5.3 Riemannian Geometry 


We shall now relate the shape tensor to the metric tensor and Christoffel connection. The metric 


tensor is defined by 
Gig = Ci C5 (5.68) 


and the Christoffel connection by . . 
I = (Dyeg) e (5.69) 
so that the components of the covariant derivative are given by 
(a- Db) - e' = a! (Dj; (b*ex)) - 
= a! (0;b' + 4,6") (5.70) 


The T% jk can be expressed in terms of the g;; by considering the following relations. First, the 
I, are symmetric in the j and k indices. 


Il, — 4; = (Dex — Deej) -e’ = 0 (5.71) 
Second, the curl of the basis vectors is given by (equation 5.66) 
DA Q= DA (gije’) = (Dgi;) /\ e/ (5.72) 


By equation 5.71 we can write 


. ee 
ee = 3° PDs ep + Dye;) 
Lite 
Sage *. ((e;-D) ex + (ex - D)e;) (5.73) 
Now apply equation B.6 (Appendix A) and equation 5.72 to each term in equation 5.73 to get 


) 
(e;- D) ex + (ex- D) e; =e; -(D A ex) + ex: (D Ae;) + (ej - x) D + (ex: és) D 
=€52(D Ne) +e. (DA e;) + D (ax) 
3+ (Dai) Ae F) + ex ((Dgjt) Ae?) + D (95x) (5.74) 
) 


Now apply equation B.2 en A) to e; - (Don A e') and e, - (Dg, Ae’) giving in the first 
case 


ej (Dam Ae!) = (e;- (Dam) e' — (e;- e') Dara 
= ((e; -D) gia) € ra 5D gu 

= (Dj gn) e! — Dox 

= (0 


gui) od 079k; (5.75) 
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so that equation 5.73 becomes 


; 1. 
Vik = af . ((O; 9x1) e+ (Ongjt) e! — 9x3) 
fe 
= 5° . ((O; 9x0) e+ (Ongjt) e — e'O1gx3) 
i Woe 
= ae (0; 94 + On951 — U9n;) (5.76) 


which is the standard formula for the I,.. 


1 

Now define the commutator bracket [A, B] of the multivectors A and B by (note there is no 5 
factor) 

[A,B] = AB-—BA (5.77) 


Now form equation B.2 (Appendix A) and use the Jacobi identity (equation 1.100) to reduce the 

double commutator products on the r.h.s. of the equation 

[D;, D;| A= 0; (0;A + S; x A) + S; x (O;A+ Sj x A) — 0;(0,A + §; x A) 185 x (O,;A + S; x A) 
= (0,5; = 0; Si) x A+ GS; x 55) xA (5.78) 


However (see equations 5.38 and 5.39, page 153) so that $; = —(0;I)I~', and 
(05; — O5;) = —O; (0,1) I~) + 0; ((A:L) I~") 
= (O;1) (O:L) I~* — (8,1) (O;1) I~ 
= (0,1) I7* (0,1) I7' — (0,1) 17" (0, I 


where we have used that J~' and the partial derivatives of J commute to reduce the second line 
of equation 5.79 


The commutator of the covariant derivatives defines the Riemann tensor 

R (a,b) = P (S (b) x S(b)) (5.81) 
Since R (a, b) is a bilinear antisymmetric function of a and b we may write 

R (aA b) =P (S(b) x S(a)) (5.82) 


or 


RleNe;) =P. (8 (e;) * 5 (e;)) (5.83) 
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Since both S (a) and S'(b) are bivectors we can use equation B.13 (Appendix A) to reduce 
S(b) x S(a) 
S(b) x S(a) = (e* AP (b- Oex)) x (e' AP (a- de;)) 

= (c* -Pi(a: de))) P. (b- Oex) A es (eS . e') P. (b- Oex) \ Px (a: 0ex) 

+ (P, (b- de;,) -e') e* AP, (a: Oex) — (Py (b+ Dex) - Pi (a: Oex))e* Ae’ (5.84) 
In equation 5.84 the first and third terms are zero. The second term is entirely outside the 
tangent space and the fourth term is entirely inside the tangent space. Also note that since the 
second term consists of bivectors that are entirely outside the tangent space that term commutes 


with all multivectors A in the tangent space so that the commutator of the second term with A 
is zero. Thus the Riemann tensor reduces to 


R (ad b) = — (Pi (b- Oeu)- Pi (a: Oey)) e* Ae” (5.85) 


or 


R (e; \e;) = — (Pi (Ojeu) - Pi (Oien)) e* Ae” (5.86) 


To calculate the Riemann tensor in terms of the Christoffel symbols note that 


[D;, Dj] ex = (Sj x Si) x ex (5.87) 
= R(e; Ae;) - ex (5.89) 


Equation 5.87 comes from letting A = e;, in equation 5.80. Equation 5.88 comes from the fact 
that S; x S; is a bivector and that the commutator product of a bivector and a vector is the same 
as the dot product. Finally we have that (5; x 5;)-e, =P (S; x S;)-e, since from equation 5.84 
P. (8; x Si) = — (e™-e') Py (e; - Oem) A P(e - Oe) so that P, (S; x S;)- ex =0. Thus 
R (e; /\ e;) "ee = Fare Dj] Ck 
= D; (T,ea) — D; (THea) 
= (O,T%,,) €a + Pi, Diea - (OP %) €a — Vi, Djea 


sO 


(R (e; A ej) ex) -e = ((O0%,) cs m4) i ((O,0'%,) + re.F%) 5 
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Using equation 5.86 we have 
(R (e; A ej) « ex) -e' = — (Px (Ojeu) - Px (Oiev)) ((e” Ae”) - ex) -e! (5.92) 
Using equation B.10 (Appendix A) to reduce ((e” A e”) - e,) - e! gives 
((e* Ae’) - e,) -e' = gtd? — g™ se (5.93) 
Substituting equation 5.93 into equation 5.92 gives 
(R (e; A e;) + ex) +e! = Px (Ojex) «Pi (Oiev) g”" — Pi (Ojeu) - Ps (Oiex) 


=P, (djex) +P. (dje') — Px (je) - Px (Oiex) (5.94) 
because 
P. (Ojev) g” = Pi (9° Ojeu) 
= P, (8; (g"ev) — (8i9") ev) 
= P, (8 (g"ev)) 
=P, (d;e") (5.95) 
Finally 


Risk = PL. (O;ex) , PL (d,e') mas Py. (d;e") i PL (O;ex) 
= OT, +T5Fy — Ox — THE ie (5.96) 
Which is the standard form of the Riemann tensor in terms of the Christoffel symbols. Note that 
Rijn = Rijn Gul 
= P. (djex) - Pi (Oier) — Ps (Ojer) - Pa (Oiex) (5.97) 
From equation 5.97 and equation 5.53 (0;e; = Oj;e;) we can see that the symmetries of the co- 
variant Riemann tensor are 
Rijn = STG ets Rizk = Tres and Rijn = Ryuj 
To prove the first Bianchi identity form R (e; A e;) - e, and use equation 5.54 (Dje, = D,e;) to 
get 
R (e; /\ e;) "€~ = D,Djex = D Dee 
= DDyé, = D Dre; 
= [Di, Dy ej — [D:, Dy e; + Dy (Die; =; D,e;) 
= R(e;, A ex) -e; — R (ej A ex) - ei (5.98) 
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now defining the function F (a,b,c) by 
F (a,b,c) =a-R(bAc)+ec-R(arAb)+b-Ri(cAa) =0 (5.99) 


However F' (a, b,c) is a linear function of a, b, and c. Also F' (b, a,c) = —F (a,b,c) and F (a,c, b) = 
—F (a,b,c) so since F is antisymmetric in all arguments we may write 


F (a,b,c) = F(aAbAc) (5.100) 


2 
n*(n—1)(n—-2 

a = ( ) ) scalar coefficients. Since the Riemann 

tensor is a bivector valued function of a bivector the degrees of freedom of the tensor is no more 


ae i as 
pide ea) 


Thus equation 5.99 contains n( 


and equation 5.99 reduces the degrees of freedom by n(3) so that the total 


degrees of freedom of the Riemann tensor is 
n(n—1)\? n 1 
ae) = (3) = ri (n° = 1) ; 


5.4 Manifold Mappings 


One way of illuminating the connection between geometric calculus on manifolds and the standard 
presentation of differential geometry is to study the effects of mappings from one manifold to 
another (including mapping of the manifold onto itself). Let f : M — M’ define a mapping 
from the manifold M to M’. For our purposes f is a diffeomorphism. That is f and all of its 
derivatives are continuous and invertible. Thus we can show that the tangent spaces of M and 
M’ are of the same dimension. We shall denote the image of « € M as xv’ so that f(x) = 2’. 
Then if x (A) defines a curve on M, then f (# (A)) = 2’ (A) defines a curve on M’. In summary - 


M: Manifold embedded in vector space V (M Cc V) 
s': Coordinates of manifold M such that x (s',...,s") € M CV, dim(V) >r 
u;: Basis vectors for embedding space V (x = x*u; € M) 

x(A): Trajectory in M such that x (A) = x(s!(A),...,8" (A)) 

M': Manifold embedded in vector space V’ 

Basis vectors for V’ 

f:  Diffeomorphism from V to V’ 

x’: Element of M’ as an image of x € M (2’/= f(r) Ee M'CYV’) 

f’: Component of f in V’ such that f (x) = f’ (x) ui 
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Since is in the tangent space of M and (remember that the dot product comes before the 


co 
TD 
geometric product) 


da! d 
T= SF (@(A)) 


_ (F 0) f (2) 
Z (# ae -o,) (f(x) uw) 


2 (se cay) (f? (x) ui) 


ds' 
= (Grater) (1 a) a) 
ds' Of? , 
= ———u, 101 
dy asi one 
is in the tangent space of M’. Thus if a is in the tangent space of M at point x then 
a =(a+0) f(a) = fase) =F @) (5.102) 


is in the tangent space of M’ at point x’ and the frame (basis) vectors for M at point x map 
from one tangent space to the other via 


6. = (6-0) f(s) = flea) —f(e) (5.103) 


where in the rhs of equations 5.102 and 5.103 we suppress the position dependence, x, in the 
linear functional f. Since f (x) is invertible for all derivatives there is no e; such that f (e;) = 0 
(the dimension of the image tangent space is the same as the original). Note that we actually 
have e; (x) and e} (x). An example is shown in figure 5.1 


The cotangent frame, e’, in M is mapped to e” in M’ via the adjoint of the inverse 
ea rile), (5.104) 
This is simply proved by noting (remember that by definition a- f (b) = b- f (a)) 
ee? =f (e;) - F (e’) 
=o fF (e,) 
=e! -e; 


hs (5.105) 
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M: {ex = cos (s?) (cos (s') uy; +sin (s') uz) + sin (s?) U3: |s?| < ae < =} 


M!' 


M’:{ F (2) = cos (s!) uf + sin (s!) uf + tan (s?) uf [3° < Shae < 2} 


e} (x) = f (e:) = —sin (s') uy + cos (s") uh 


Figure 5.1: Mercator mapping of sphere to cylinder manifold. Top and bottom of sphere excluded. 
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The exterior product of two tangent vectors, e; \ e;, maps to 
e; Ne; +e, Ae, =f (e; Ae) (5.106) 


since f is linear in tangent space arguments. Likewise if J’ is the “unit” pseudoscalar for the 
tangent space of M’ at f (x) (I = +1) and J the “unit” pseudoscalar for the tangent space of 
M at x (I? = +1) then 

f (I) = det (f) I’. (5.107) 


Since f is invertible (the tangent space and its image are isomorphic) we can apply the definition 
of determinant in equation 1.79 to equation 5.107. For cotangent vectors e’ and e/ 


eNeTI SE eV At (8) te Ae) (5.108) 
Since the derivative of a scalar field, ¢, is a cotangent vector 0¢ = e'0;¢ and (x) = ¢' (z’) we 
can write = on, 
=f! (0) =F (e'0,) =e"0,, (5.109) 
Since D is also a cotangent vector we also have 


DST tp): (5.110) 


For the directional derivative of a scalar field 


=(a-d)¢ (5.111) 
or for a vector field 
(a’- 0’) b' = (a- 0) f (b) (5.112) 


The covariant derivative is constructed using the projection operator, P, which contains a con- 
traction with the pseudoscalar J. Thus the covariant derivative depends upon the metric encoded 


by I (2). 


Consider the following operation where a and 0 are tangent vectors (from now on the parenthesis 
around the dot operands are implied) and using equations 5.49 and 5.50 we get 


a:-0b—b-0a=a-Db—b-Da—a-S(b)+b-S (a) 
=a-Db—b-Da (5.113) 
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since a-S'(b) = b- S (a) by equation 5.56. Now define the Lie derivative of a with respect to b by 


L,b=a-0b—b-0a 


We will show that 
Lab Lib! =F (L,b). 


First note that 
(a) (0) ~ (b-0) F(a) =f ((a- 0) (b-0)a) + (a-d) #(®) ~ (0-4) F(a) 
Note that since f (a) is the differential of f (~) we have (Oe; — 0;e; = 0) 


0 = (0,0; — 0;0;) f (x) = Of (e;) — O;F (es) 
0= O,f (e;) = O;f (e;) + f (Oje; = O;e;) 
0= Oie!, = Oje; 


so that Oje, — Oje; = 0 and df (e;) — O;f (e;) = 0. Thus 


(a-d) £0) - (0-8) F@) 


(ake . e'd,) f (b’e;) — (vie . ed}) f (a‘e;) 
a’ df (b’e;) — bO;F (a’e;) 

= a'b’0;f (e;) — b’ a’ dz (e;) 

= a'pi (af (e;) — jt (e:)) = 


and 


(a- 0) f(b) — (b- 0) f (a) =f ((a- 0) b— (b-O)a) 
and by equations 5.112 and 5.119 we have 


Cite 


and the Lie derivative maps simply under f. 


k; 


Since e’* - e/, = df and e’* - e} = 6f we have (using equation 5.117) 


8; (e* - ef) — d; (e* - ej) =0 
e"* . (G,e', — Oje;) + €, - Gel — ef, - je =0 
ef, - Oje"* — ef - Aje"* = 0 


f (e;) - O,f-1 (e*) — f (e;) - 0;f-! (e*) =0 


(5.114) 


(5.115) 


(5.116) 


(5.117) 


(5.118) 


(5.119) 


(5.120) 


(5.121) 
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Equation 5.121 inplies that P’ (F(a) LS (e’) = 0 as can be shown by expanding 
f-1 (8) A f-1 (e*) = FT (e0,) A FH (e*) 
=e", \e™ (5.122) 


Since equation 5.122 is a grade two multivector we can expand it in terms of the blades (compo- 
nents) e”’ Ae’ as follows from equations 1.75, 1.76, B.10, and the fact that 0; is a scalar operator 
that commutes with /. 


P! (F(a) AF (e*)) = (eh, Net) - (€%: Ae) (eA c™) 
= ((ef, Ae): (e" A de) (e" Ae’™) 


e,,, - Ae") (e,-e”)) — ((e, -e*) (e} - Gje’*)) (e” Ae’) 
(el, Oye") — 5, (e,- Oe"*)) (e" Ae’) 

= ((e,,, - He) — (e,- One’”)) (e” Ae”) (5.123) 
But the coefficient of e” A e’” is the same as equation 5.121 so that P’ (FI (0) Af-t (e’)) =¢. 
Note that since we expanded f~! (9) A f=! (e*) in terms of e” A e’” (cotangent vectors in M’) it 
was automatically projected into the tangent space of M’ at 2’ = f (x). 
Since D! = P’ & (0)) we have 

Dike SD KEE) S06. (5.124) 


and using the product rule the exterior derivative of a blade formed from cotangent vectors is 
(just move the term in the product being differentiated to the front using the alternating property 
of the “A” product) 


DENCE Dvn Net) =, (5.125) 
Now expand a grade-r multivector field A(x) on M in terms of the cotangent vectors 
Ae) = Apps Gye Ne? Wise” (5.126) 
Then . . . 
DA A(a) = (DAijria..i, (2)) NC AC? A... Ne” (5.127) 


since the exterior derivate of a scalar field is the same as the covariant derivative and that the 
exterior derivative of a basis blade is zero. Likewise 


D!' A Al (2) = (D' Aain..z, (Z)) NEM NCPR A... ANC (5.128) 
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so that for a general multivector field A(x) on M we have 


DAAWR DAA =F (DAA). (5.129) 


5.5 The Fundmental Theorem of Geometric Calculus on Manifolds 


The fundamental theorem of geometric calculus is simply implemented on manifolds. In figure 5.2 
a simplicial decomposition of a surface defined by a close curve on a spherical manifold is shown. 
The only difference between the derivation of the fundamental theorem of geometical calculus 


Figure 5.2: Simplical decomposition for surface defined by closed curve on spherical manifold. 


on a vector space (which we have proved) and on a manifold is that the pseudo-scalar is not a 
constant but is now a function of position, J (x), on the manifold. If we are not on a manifold 
the pseudo scalars corresponding to the oriented volume of each simplex are proportional to one 
another (they are equal when normalized) so have the same orientation. For a manifold the 
orientation of J (x) can change with the position vector x. In the case of figure 5.2 I (x) is a 
bi-vector defined by the tangent space (tangent plane) for each point on the sphere. 
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Now consider the directed volume element (in the case of figure 5.2 an area element) for each 
1 
simplex in figure 5.2 given by AX = mia A... ex (k = 2 for figure 5.2). As the volume (area) 
of AX > 0, AX « I (2). , 


In equation 4.75 where on the l.h.s. of the equation we are integrating over the boundary of the 
simplex where f is a linear approximation to an arbitrary multivector function. 


¢ f (2) = f¥- (AX) 
A(x) (x) 


consider the operator V- (AX) where we have left the dot on the V to emphasize that it is not 
differentiating the AX. On a given simplex 

<@) 

V=e—— 5.130 

at (5.130) 
where the ’’s are the simplical coordinates (equation 4.19) and the e’’s are the reciprocal vectors 
to the e;’s that define the simplex. In the case of a manifold as the volume of AX — 0 the e;’s 
(since the e’’s define the same subspace as the e,’s) define a pseudoscalar that is proportional to 

oO 

I(x) so that e'~— is the projection of the geometric derivative, V, from the embedding vector 
space of the manifold to the tangent space of the manifold at point x. Thus in the case of a 


manifold as the volume of AX — 0 we have 0 = e'—. 


Or 
Note that 
O(AX) =0- (AX) +0A (AX) (5.131) 
but OA (AX) = Osince 0 is within the subspace (tangent space) defined by AX. The fundamental 
theorem of geometric calculus as applied to a manifold is 


§ Ls)= [ L(v-ax) = | (dax) (5.132) 


One can write dX = I (x) |dX| and note that O or V in equation 5.132 do not differentiate J (x). 


5.5.1 Divergence Theorem on Manifolds 


Let L(A) = (JAI) where J is a vector field in the tangent space of the manifold and substitute 
into equation 5.132 to get 


¢ (JdSI-') = ((JdaxI-") + (JdaxT-1)). (5.133) 
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Now reduce equation 5.133 by noting that ndS = I |dS| where n is the outward normal to the 
surface element dS. We can define an outward normal in n-dimensional manifold since the grade 
of dS is n — 1 and it defines a subspace of the tangent space of dimension n — 1 for which a 


unique normal exists. This gives 


dS = 1 \ds| 
n 
“i. Ids _ 
(JdSI~*) = Sp (snd *) 
ds 
J-n|dS| 
4 


Also 
JOdXI-* = JOI|dX|I- 
= Jo|dx| 
( jadx I“) = (jd |dX|) 
= (Ja) \dX| 
= 9-J|dX|. 
Finally, since both J~' and dX are proportional to I 
( JOI-'dX ) af (Jail) dX | 
1 Lops : 
=+5 (Ja (i7 zs 11)) \dX| 
~ 4h 2 
= +5 (J0(P)) lax, 
=0 


so that the divergence theorem is 
¢ “. J \d8| = / OJ \dX| 
av V 


5.5.2 Stokes Theorem on Manifolds 


(5.134) 


(5.135) 


(5.136) 


(5.137) 


(5.138) 


(5.139) 


Assume that the manifold tangent space dimension is s+ 1 and B, is a grade r multivector 
field. For clarity we denote the grade of volume and surface elements with a subscript on the 
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d’s so that dX = d,,,X and dS = d,S. Now let L(A) = (BrA)), 9 so that the application of 
equation 5.132 gives 


¢ (BrdeS)), 4 = i: (B,Ads41X) 
av V |s—r| 


grade |s — r| grade r+ 1 
Boi = i B.AO +  B,-0O dg41X 
av V bee ora 


grade r — 1 


|s—r| 


= cay f (OA B,) ds41X + [ ( (2, ; 3) ds41X 
uu a 


lowest grade is |s — r| |s—r| lowest grade is |s — r + 2| ee 
ey i (0. By) + dey X (5.140) 
V 
However 
(OA B,) + ds41.X = (DA B,) + ds41X (5.141) 


since ds4,X = I,41 (x) |ds41X| and the dot product of any component of 0 A B, that is not in 
the tangent space defined by /,+, (x) is zero so that 


¢ B, -d,S (3, AD) Oe cm ay f (DAB)> dak (5.142) 
OV Vv Vv 


The divergence theorem is recovered when r = s—1. This is important for constructing conserva- 
tion theorems in curved spaces. Equation 5.142 is the most general application of equation 5.132 
that allows one to replace 0 with D (covariant derivative) on a manifold. 
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5.6 Differential Forms in Geometric Calculus 


5.6.1 Inner Products of Subspaces 


Start by considering products of the form where A, and B, are r-grade blades 
At. B,=(a,A...A ar) (bl A... A br) (5.143) 


Both A, and B, define r-dimensional subspaces of the vector space via the equations x A A, = 0 
and « A B, = 0. We show that if A, and B, do not define the same subspace then At -B, = 0. 
Assume that they do not define the same subspaces, but that the intersection of the subspaces 


has dimension s < r so that one can have an orthogonal basis for A, of {e1,...,€s,€s+1,---,€r} 
and an orthogonal basis for B, of {e1, Se5G Bae pig ats se Then the geometric product of Al 
and B, is 


/ / 
A'B, = Ops se Ogg lg: CP Or rls O oi Sere. 


(CUBE; 25:6) Cera Oe uence (5.144) 


where the quantity in parenthesis is a scalar and the other factors a blade of grade 2(r— s). 
Thus 
Als By =(AlB,) = 0, (5.145) 


If A, and B,. define the same r-dimensional subspace let {e;,...,e,} be an orthogonal basis for 
the subspace and expand A, and B, in terms of the orthogonal basis vectors and reciprocal 
orthogonal basis vectors respectively - 


a; = (a; - e’) e; (5.146) 


and 
b; _ (b; : e;) e) (5.147) 


Let the matrices of the coefficients in equations 5.146 and 5.147 be denoted by [a; - e’] and [}; - e;]. 
Then we can expand A, and B, as 


A, = det ([a; . e']) Ci5.n€p (5.148) 


and 
B,. = det (bj - e;]) et..." (5.149) 
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and (since the determinant of a matrix and determinant of the transpose of matrix are equal) 


Al B,. = det ([a; e! ) det ([b; - ej]) (er ...e1€"...€") 
= det ([a; - e’}) det ([b; - e;]) 
= det ([a; - e’]) det ([b, : el") 
= det ([a; - e’]) det ([b; - e:]) 
= det ([(ai e*) (b; - ex)]) (5.150) 
But 
a; b; (a; e*) e;, * (b; - e1) e 
= (a; -€*) (bj -€1) 5 
= (a; - e*) (bj - ex) (5.151) 
So that 
Al. B, = det ([a;-bj]) . (5.152) 


From our derivations we see that if A, is a general r-grade multivector (not a blade) we can 
always find a r-grade blade such that 


BEA OE icy ID) AG Des 2) @ be Rs Abe e (5.153) 

Now consider the relation between the basis {a1,...,@,} and the reciprocal basis {at,...,a”} 
for an n-dimensional vector space where in equation 5.154 r <n and 1 <ij,jm <n 

(ij: Nava.) * (a A... Aa’) = det ([a;, -a’”]) = det ((62"]) (5.154) 


where 21 < tg < +--+ <i, and jy < jo <-+: < j,. Equation 5.154 is zero unless 7; = 7; for all 
1<l<rso that Sa ; 
GE hGe) OW eRe) =O (5.155) 


11 “te 
since if the index ordering condition is satisfied for the 7;’s and the j,,’s and there is an index, J, 
such that 7; # j; then the two blades do not define the same subspace and the inner product is 
zero. The square of a blade is given by 


(a, A...Aa,)? = A Gr) + (a1... A ay) 5.156 


5.157 


5.158 


A 
Se (ay Ac Nae)! Maher. van) 
) 
) 


(5.156) 
(5.157) 
ae , (5.158) 
(5.159) 


5.159 
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5.6.2 Alternating Forms 


If V is a vector space then an r-rank tensor is a multilinear map 
Pid ves® 1 QVOR (5.160) 


where (@ is the cartesian product of vector spaces) 


®V=Ve---@V (5.161) 
Ww "—~ -_—_ 
= r times 


so that if vu; € V the tuple (v,...,ur) € &V. 


The sum of two r-rank tensors is the r-rank tensor defined by 
(A, + B,) (v1,.--; Ur) = Ap (U1,---5 Ur) + By (U1,.--5 Up) - (5.162) 
The tensor product of rank r and s tensors is a rank r + s tensor defined by 
(Ap OB) (Oise 5 Ups) =] Ay Ginnss (0p) Bs Uriagees tts) a (5.163) 
An r-form (when we say something is a form from now on we mean alternating form) is an r-rank 
tensor with the property (no summation in this case) 


a. bp 


pO CeO Se Oe (Ugo, Vie) (5.164) 


where ¢(!”;"" is the mixed rank permutation symbol. We can always construct an r-form from an 
r-rank ie isor via 


Glin wees Oe) = ee Ge Ae Wasa Oe) (5.165) 


qeesy 


In the geometric algebra a simple representation of an alternating r-form is 
Op (U1,...,Ur) = AL (ur A... A Up) (5.166) 


where A,. is a grade-r multivector. Since the grade of A, and the grade of (v; A... A v,) are the 
same the inner product results in a scalar and also since (v; \.../ v,) is a blade it alternates 
sign upon exchange of adjacent vectors. 
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The basic operations of the “Algebra of Forms” inherit the sum and cartesian product operations 
from tensor since they are tensors. However, in order to construct an algebra of forms we need a 
product of two forms that results in a form. The tensor product of two alternating forms is not 
an alternating form, but since we know how to convert a tensor to a form we define the exterior 
product of an r-form and s-form to be the r+ s form 


a,ABs = € (a, ® Bs) (5.167) 
as an example 
ay (v1) ABy (v2) = a1 (v1) Bi (v2) — a1 (v2) Bi (v1) (5.168) 
so that 
Oty (U1, -.., Up) ABs (Unga, «+ +5 Urs) = (Ap A Bs)! (up A... A Up 4s) (5.169) 
and 
arABs = (—1)" BsAa,y (5.170) 


The interior product of an r and s-form is defined by (r > s) an r — s form (note that we 
are distinguishing the inner, -, and exterior, A, products for forms from the geometric algebra 
products by using boldface symbols) 


Bs + Or = (Bs- A,)'+ (Us41 A... A vp) (5.171) 


A r-form is simple if A, is a blade (A, = a, A... A@,). 


5.6.3. Dual of a Vector Space 


Let V be a n-dimensional vector space. The set of all linear maps f : V — FR is denoted V* 
and is a vector space called the dual space of V. V* is a n-dimensional vector space since for 
f,g:V—7R, x2 €V,a€R and we define 


(f +9) (x) = f(x) +9 (2) (5.172) 
(af) (x) = of (2) (5.173) 
0(x) =0 (5.174) 


Then by the linearity of f (x) and g(x), (f +g) (x), af (x), and 0 (2) are also linear functions 
of x and V* is a vector space. 


Let e; be a basis for V and define a € V* by 
oa (e;) = 5° (5.175) 


J 
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so that if x = x’e; € V then 


gia) =o (ee) =a (5.176) 
so that for any f € V* 
f(z) =f (we) =2'f (e:) (5.177) 
Now assume that 
0 = ajo" (e;) (5.178) 
= a6; (5.179) 
= Qj (5.180) 


j(@)=7 @e) (5.181) 
= x'f (e;) (5.182) 
= f (e;) 0° (x) (5.183) 
= (f (e;) 0") (2) (5.184) 
Thus 
f=fle)a' (5.185) 


and the o’ form a basis for V*. If the e; form an orthonormal basis for V (remember that 
orthonomal implies orthogonal and orthogonal implies that a dot product is defined since the 
dot product is used to define orthogonal) then 
o' (x) =e; 2. (5.186) 
Since for 1-forms, a: V + R, we have a € V*. The most general 1-form can be written 
a(v) = qa" (v). (5.187) 


If a is a 2-form, a (v1, V2), the bases are o*Ao/, and the most general 2-form is written as 


a (v1, v2) = S> ayo" (v1) Ao? (v2) (5.188) 


i<j 


since o'Ao*t = 0 and a’Aai = —o*Ao! from the definition in equation 5.167. 
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5.6.4 Standard Definition of a Manifold 


Let M” be any set? that has a covering of subsets, M" = U UV U... such that 


1. For each subset U there is a one to one mapping dy : U > R” where gy (U) is an open 
subset of R”. 


2. Each @y (U NV) is an open subset of R”. 


3. The overlap maps 
fuv = ¢v°¢y' : bu (UNV) > R” (5.189) 


or equivalently the compound maps 


-1 
du (UNV) 23 Mn 2% RP (5.190) 
are differentiable. 


4. Take a maximal atlas of coordinate patches {(U, dy), (V, dy), ...} and define a topology for 
M”" by defining that a subset W C M” is open if for any p € W there is a (U, dy) such 
that pe U CW. 


If the resulting topology for M” is Hausdorff and has a countable base we say M” is an n- 
dimensional differentiable manifold (look it up since I don’t know what it means’). 


Figure 5.3 (page 94) shows the relationships between the Manifold and the coordinate patch 
mappings. 


If f : M” > R is a real valued function on the manifold M” it is differentiable if fy = fo On 
is differentiable with respect to the coordinates {xj,...,x7,} for any coordinate patch (U, dy). 
The real scalars {xj,,..., x7, } (denoted by the tuple x = (xj,...,xf,)) form a coordinate system 
for gy (U) C R”. In the future we shall simply say that f is differentiable if fy is differentiable. 
Likewise we will shall usually omit the process of replacing f by its composition fo¢,', thinking 
of f as directly expressible as a function f (x) = f (xj,...,x7,) of any local coordinates. 


Now let pe UNV Cc M” and 


te (Gia eth Hor) and ay = (fyi) = Oyo) 


?This section is based upon sections 1.2 and 2.1 in “The Geometry of Physics, An Introduction (Second 
Edition),” by T. Frankel 
3G. Simmons, Topology and Modern Analysis, McGraw-Hill,1963 
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Abstract set M” not necessarily 
embedded in a vector space. 


Figure 5.3: Fundamental Manifold Mappings. 
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Then 
ty = by ° by" (zy) = tu (zy) and zy = gy 0 $y (tu) = zy (2). 


Let Xy = (Xq,..-,X7) and consider the linear approximation to dy o éy' (zu + hXu) where h 
is a small number. Then in the linear approximation 


fa) a 
ov 0 by (tu +hXy) = ty (au) +h Ea Xp = ty + hXy (5.191) 
a 
so that . 
be . OS ke 
De ani, 8 (5.192) 


Ont 

where =] is the Jacobian matrix of the coordinate transformation from one coordinate patch 
bs 

to another and Xf is the transpose (colume vector) of the tuple Xy (row vector). 


5.6.5 Tangent Space 


Now consider the case of a vector manifold in Euclidian space. The definition of the tangent 
vector generalizes the concept of the directional derivative in Rk”. if X, is a vector at point 
peER" and f:R" > Ris a C™ function in the neighborhood of p then define 


Xp (f) = Xp: VEI,- (5.193) 
If f and g are C® and a, 8 € R we have (from the properties of - and V) 


Xp (af + Bg) = aXp (f) + BX» (9) (5.194) 
Xp (f9) = f (p) Xp (9) + 9 (Pp) X> (Ff) (5.195) 


Let C™ (p) denote the set of real functions that are C® on some neighborhood of p € M”. The 
generalization of the tangent vector to an abstract manifold at a point p € M” is defined as a 
real valued function X, : C% (p) + R such that (f,g € C™ (p) and a € R) 


Xp (f +9) = Xp (f) + Xp (9) (5.196) 
Xp (af) = aX; (f) (5.197) 
Xp (f9) = fp) Xp (9) +9 (Pp) Xp (Pf). (5.198) 
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A tangent vector is often called a derivation on C'® (p). For a given coordinate patch (U, éy) we 


0 
define the tangent space basis vectors (sr | by 


( @ ) 7 = Olfe ds’) 


i Fi 
Oxy Oxi; 


(5.199) 


gu (p) 


O 
Since they are partial derivatives the (sr | trivially satisfy equations 5.196, 5.197, and 5.198. 
vy 


P 
The general tangent vector is then given by the scalar linear differential operator 


sof 70 
= ¢ — 2 
Xp = XU (ser), (5.200) 
Thus we require for X, 
of /{ O 
XH=xb (ge) r=xt (SS) (5.201) 
Oty p ¥ Ont, m 
for allp @ UNV. Thus 
3 tO -{ O 
XT = Xi, (| —— 5.202 
i (aap) f= xt (ger) J (5.202) 
= [Ox 4) ; [{ O 
a r | (sr | fare (= f (5.203) 
id (a é Og ‘ : Oxi, ; 
and 
j 
X= xt (3 (5.204) 
P 


is the transformation law for the components, Xj,, of the tangent vectors. 


We have established an isomorphism between the tuple (Xj,..., Xj) that transforms as in 
equation 5.192 and the scalar linear operator X, so that 


. Oo 
ee, ai 2 
(Xu, ---, XG) SOU ane (5.205) 
since 9 5 
i OF x? vi VE: M">R (5.206) 


Uni “Vea”? 
Oxi, Ox\, 
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implies equation 5.192. A vector field on a manifold is then a tangent vector with coefficients, 


X‘ (p) that are functions of the position p on the manifold so that (xy = (1j,.--, 27%) 
X (p) = Xy (av) = Xi (20) es (x; at) ae (5.207) 
U U U U Ax, : U Uo: +9 tu Oxi, : . 


and the coefficients Xj; (ry) transform under a change of basis according to equation 5.204. 


5.6.6 Differential Forms and the Dual Space 


If f : M" + R we define the differential of f at p< M”, df: M> — R where X, € Mj as the 
linear functional 


df (X) = Xp (f) (5.208) 


so that in general case of X being a vector field (we are supressing the the patch index U) 


a) 4d 0 
a (Xp) =f (x (55) = x'(y) 2 (5.200) 
If f (p) = 2" (p) are the coordinate functions we have 
; 4) Ox" ‘ 
dev (4) | ao] 6 (5.210) 


- da’ (X (p)) = da’ [« (p) (2) ) poe (4) | oe — 


Consider what this means in a vector manifold where the point p on the manifold is a vector 
in the embedding vector space. Then the coordinate functions x‘ (p) can be inverted and the 
manifold defined by p(a',...,2"). Knowing the tuple (z',...,2”) lets one calculate the vector p 


on the manifold. Then the basis tangent vectors are e; = = and the general tangent vector is 


Ox 
0 
Kp = xi. Now let f (p) = f (z!,...,2") :R" > R be a function from the manifold 
nt 


to the scalars and denote the directional derivative of f as 


Of 
Ox? 


Pp 


df (Xp) =X- Vfl, =X" (5.212) 
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which is the same as equation 5.207 and the justification for defining the tangent vectors as the 
linear scalar differential operator X, = X' Bail’ What is called the differential operator in the 
et 


Pp 
language of differential forms becomes the directional derivative on a vector manifold. 


The dx’ form a basis for the dual space M;,". The most general element of M7>* can then be 
written as 
a= ode". (6.213) 


The linear functional a € M>* is called a covariant vector, covector, or 1-form. If a; is a function 


of p it is called a covector field. 


qa is an r-form if 
@ = 0;,,4 (p)da A... Adz (5.214) 


where 7, < ig <---<1,. The exterior derivative of a is an r+ 1 form defined by 
da = (dai,..i, (p)) Adx™ A... Adz. (5.215) 


For example if a = a,dx! + agdx? + a3dx* then the exterior derivative of the 1-form a is 


da = Oay _ om da Ada? + OG om da Adx® + OBS 2 2ae dx*Adx® (5.216) 
Ox! = Ox? Ox! = Ox Ox? —s Ox 


and if a = aydx! Adz? + aygdx!Adx? + ag3dx2Adx? then the exterior derivative of the 2-form a 


is 
7 Ja12 Ja13 Ja23 
oe (= Ox? + oa? 


) dx'Adx* Adz’. (5.217) 


5.6.7 Connecting Differential Forms to Geometric Calculus 


To connect differential forms to the geometric calculus we must establish a correspondence be- 
tween the fundamental theorem of geometric calculus and the generalize Stoke’s theorem of differ- 
ential forms. For simplicity, since we do not wish to use the most general formulation of Stoke’s 


theorem we will use the treatment in ”Calculus on Manifolds” by Spivak (http://faculty. 
ksu.edu.sa/fawaz/482/Books/Spivak_Calculus/%200n/20manifolds.pdf). In this treatment 


a manifold is embedded in the Euclidian space ” so that the manifold and its boundry in Spivak 
are submanifolds of #” as are also the manifolds V and OV which are the domains of integration 
in the fundamental theorem of geometric calculus (eq 4.81). The fundamental theorem of geo- 
metric calculus is more general than the generalized Stoke’s theorem since the linear functional 
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in the integrals can be multivector valued fields. The appropriate form of geometric calculus 
theorem is (for explicitness subscripts in parenthesis indicate the grade of a multivector field or 
the dimension of a manifold or it’s boundary, of the rank of a differential form) 


¢ Q(r) AS (r) = cay [ (VAQ@) «drat: (5.218) 
Vir) Vort4) 


Equation 5.218 corresponds to the generalized Stokes theorem in differential forms 


¢ W(7) = / dw(y41)- (5.219) 
OVir) Vor+1) 


The differential forms w(,) and dw(,+1) are alternating tensors (antisymmetric) of ranks r and r+1 
respectively. We now must expand equations (5.218) and (5.219) in terms of their components 
to show equivalence. 


Chapter 6 


Multivector Calculus 


In order to develop multivector Lagrangian and Hamiltonian methods we need to be able to take 
the derivative of a multivector function with respect to another multivector. One example of this 
are Lagrangians that are function of spinors (which are even multivectors) as in quantum field 
theory. This chapter contains a brief description of the mechanics of multivector derivataives. 


6.1 New Multivector Operations 


Define the index i,,.} = (71, %2,...,%,) where r < n the dimension of the vector space and P (icy) 
is the union of 0 and the set of ordered tuples 74,} = (71, 72,..., 7) defined by 


ac ae { (#1705 +2252) such that 21 <%< <<< <4, and 1<4,;< nforl<7<r} |” 
(6.1) 
Essentially P (igry) is an index set that enumerates the r-grade bases of the geometric algebra 
of an n-dimensional vector space where 0 < r <n. Then define the basis blades 


Cie Oy, NON se NCR (6.2) 


where Cig = Co = 1 and 
ett} =e™ ANe™ A... Ae” (6.3) 


where e’(9} = e? = 1 and the e,’s form a basis for the multivector space. With these definitions 


= 6," (6.4) 


Ur} se. 
. CF Ery J{r}" 
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The multivector X can now be written as 
a 2 a AM Cie = De » Xie et. (6.5) 
r=0 try EP (try) r=0 itr} €P(igry) 


We now generalize the Einstein summation convention so we can write X = X REY 8. Now it 
is clear that X is a 2” dimensional vector with components X‘t"}. For example for n = 3 


 < (X)o (X), (X)» (X)3 (6.6) 
= X°+ Xe, + X%eq + X%e34 
Xe; ix €2 + Xe; A €3 + XP 65 ix €3 + Xe ix €2 ix €3. (6.7) 


From the properties of the dot product of blades developed in section 5.6.1 we have for r > 1 
Cie : CF en = (€;, DekscdN €;,) * (en INSdcads en) (6.8) 


= 2 
io Oigryitry City) 


Now define the scalar product of two multivectors A and B by 


Ax B=(AB). (6.10) 
Then 
(A), *(B), =Oifr4As (6.11) 
(A) (Be =A), BAB) ee A) ee (6.12) 
(A)o - (B)o Pr (Aho (B)o = (A) (B) (6.13) 
Ax B=) (A), *B=5- (A), *(B), = (A) (B) + 2 (A), - (B), (6.14) 


by grade counting arguments and the orthogonality properties of basis blades of grade greater 
that 1. Note that since (A),, is a sum of blades each defining a different subspace we have by 
equation 5.145 that the blades that compose (A),. are orthogonal to one another under the * and 
- operations. Also 


Ax B=dAB) =(BAV= BA (6.15) 
Ax(aB+6C)=aA*B+BA*C (6.16) 
AxB=Alx Bl, (6.17) 
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We also have 
bps ee er pee. ee BD 
e Cjgy = EO * ey) = Oy", for r > 0. (6.18) 


We now use the scalar product to define the scalar magnitude |A| of a multivector A by 


|A? = Ate A= S~|(A),)? (6.19) 
r=0 
Now define a super metric tensor for the entire geometric algebra vector space by 
Lfor r= 
= ee ae —a | 
Gig ites = Cigny * Cita, = Ons Haya) ee p (6.20) 


On adey (c.) forr>1 


and 
1 for r =0 
Gltryits} = ettr} x elf} = §"5 gta toh F =u (6.21) 
O'tr}I{r} (c,,) forr >1 


The super metric tensors have only diagonal entries except for G 
block nature of the G tensors (G 


| en igyjay and G'I04, Due to the 
= Gtr}{s} = 0 for r # 8s) we can write 


Ur}I{s} 
= etf{r} : : — pt{r} : : 
Cin =E Ca * Cie) € (ese * ei) (6.22) 
J{r} = e, “ry x elir}) =e. Jer} x ebtr} 
e Cis g (ert? ef") Citas (ere | (6.23) 


An additional relation that we need to prove is (no summation in this case) 
COG es, ll. (6.24) 
First consider the case for r = 1 
eels =e'e;,, =e" -e;, +e" Nei, 
=1+9e;, Ne, =1+ S- g™ (e;, Aex + ei, Ae;,) = 1. (6.25) 
Ji<t1 


Now consider the case r > 1. Since e; a and e’{r} are blades that define the same r-dimesional 


subspaces they can be written as the geometric product of the same r orthogonal vectors 01, ... , Or 
to within a scale factor so that 


ete, = (eA... Ae") (e; Nee IN CR, 
{r} : 


= af (o,...01) (01 ...07) = aBo;...0% = elt} - Ci, = 1 (6.26) 


since by equation 6.26 etheje. = (fre, ) is a scalar. 
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6.2 Derivatives With Respect to Multivectors 


We start by discussing exactly what we mean we say F’ (X) is a funtion of a multivector X. First 
the domain and range of F'(X) are both multivectors in a 2"-dimensional vector space formed 
from the multivector space of the n-dimensional base vector space. Another way of stating this 
is that the geometric algebra of the base vector space is generated by the n-grade pseudoscalar, 
I,, of the base vector space . The domain of F’ is the multivector space of the geometric algebra 
defined by the normalized pseudoscalar [,, where J? = +1. Thus we can consider X to be an 
element of the 2” dimensional vector space defined by J,,. The partial dervatives of F(X) are then 


OF 
the multivectors ——— and there are 2” of them. The definition of the multivector directional 


OX" 
derivative Ox is 
F(X +hA)—-F(X 
opie eee (6.27) 
h—0 h 
Thus 
F(X) + hAttr} —_— — F(X) 
La OX} 
(Ax 0x) F = lim i (6.28) 
oa SOE 
— Ur} 
A ax (6.29) 
a Alive: J{r} 
AXE; *e aXe (6.30) 
— J{r} 
Axe aXI@ (6.31) 
so that in terms of components 
= el{r} 
Ox =e aXe" (6.32) 


Note that we have put parenthesis around (A * 0x) to remind ourselves (Doran and Lasenby do 
not do this) that A * Ox is a scalar operater in exactly the same way that a: V, a-0, and a- D 
are scalar operators. While not explicitly stated in D&L or Hestenes, * must have the same 
precendence as - and higher than any of the other product of geometric algebra. The multivector 
derivative Ox is calculated by letting A take on the values of the bases for the geometric algebra. 


Another notation used for the multivector derivative is 
Fy (A) = F(A) = (A+ dx) F(X). (6.33) 


The form F(A) would be used if it is implicitely understood that the multivector derivative is 
to be taken with respect to X. 
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We can now define the adjoint of F(A) by 


PUB) = OA (ANB) (6.34) 


This make sense when we consider 


(AF (B)) = (Adc (E (C) B) 


AxF(B) =F(A)*B. (6.35) 


When A and B are vectors a and 6 and F'(a) is a vector-valued linear function of a we have 


(a) -b (6.36) 


which recovers the original definition of the adjoint. 


Note: Consider equation 6.35 for the case that A and B are pure grade, but not the 
same grade. For example let B be a bivector, F (B) a vector, and A a vector. Then 
Ax F(B)=A-F(B) is a scalar, but then F(A) must be a bivector for F(A) * B to 
be non-zero. In general if A and F (B) are pure grade they must be the same grade 
then F(A) is the same pure grade as B and we can write A- F(B) = B- F(A). 
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Product rule - 


ee Nani 6 ee ee 
. ie ee h 


, OF OG 
| (r (X) + hating ons) (c (X) +hAie = ac] ~ F(X) G(X) 
= lim 
h—0 h 
OF OG 
a Ulp If{r 
St At aa A AM eae 
i OF OG 
—_ Ufp a ILr | 
= Aie;,,, * 1 ——— (FG + FC) 
7 (A * dx) Ge cm FG) (6.37) 
so that the product rule is 
Ox (FG) = dy (Fc + FG) (6.38) 
Chain rule - 
F XxX A))-—F XxX 
(A Oe) F(G(X)) = fin POG HAD ~ PGW) 
F OG 
= lim 
h—0 h 
_ OGI OF 
F(G(X hA'ir} ——__——_ — F (G(X 
“gg, FG + AO Sees pg 7 FG AD) 
h—0 h 
gy Ce Ok 
OX} OGII} 
,,  OGI pee OF 
— AX xt Chen *x et AGF 
‘ OG 
= ((A* Ox) G) * OGF (6.39) 


If g (X) is a scalar function of a multivector X and f (g) is a scalar function of a scalar then 


(Axx) 9) F = (Avdx) oF (6.40) 
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and 
Oxf (g(X)) = Oxo) F 


One other general relationship of use is the evaluation of 0,4 ((A * Ox) F) 


(B+ 4) ((A* Ox) F) = lim 7 


so that 


For simple derivatives we have the following - 


(A+ AB) * Ox) F) - (Ax Ox) F 


(B * Ox)(A* X) =lim a Ss ea 
h0 h 
=(AB)=BxA 
Ox (A x X) =A 
2. Ax (X * X"): 
((x +hA)(X +hA)t — XxX*) 
(A * Ox) (X * XT) = lim h 
=(XAt+AXT) = Xx ATH Ax Xt 
=2Ax Xt 
Ov XX) SOX! 
Oxt (X « X') =2X 
3. Ox (X * X): 
ge A+ A) (XA + RA) = XX) 
(A * Ox) (X *X) = lim 1 
=(XA+AX)=X*«A+A*xX 
=2A* X 


Ox (X * X) =2X 
Ox; (Xt * XT) =2Xx" 
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(6.41) 


(6.42) 


(6.43) 


(6.44) 


(6.45) 
(6.46) 
(6.47) 
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4. Ox (XT * XT): 


(A * Ox) (Xt XT) =] ((X + ha)! (X + nA)! — XTX1) 


h—0 h 
=(XtAT + AUX) = Xt AT + Ale XT 
=2A * X 
Oe (Xe A) SOX (6.50) 
Oxt (X * X) =2X7 (6.51) 


5. Ox (IxI*): 


ax (Ixt") = ox ((1x?)4) 


k 2 ky 
= 9X |)? ae! 
5 (XP) 
=k|x|*? xt (6.52) 


by the chain rule. 


6.3 Calculus for Linear Functions 


Let f (a) be a linear function mapping the vector space onto itself and the e; the basis for the 
vector space then 


= e' (e;-f (e;)) a! (6.53) 
The coefficient matrix for a linear function f (a) with respect to basis e; is defined as 
fig = es - £ (e;) (6.54) 


and 
f (a) = e'f,;a? (6.55) 
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Now consider the derivatives of the scalar f (b) - c with respect to the f;; 
Op, f (0) i as Of (fu.d*c') 
= 646;,b%c! = bic’. (6.56) 
Multiplying equation 6.56 by (a- e;) e; gives 


(a+ e;) E07, (b) -c = (a- e;) eb’! 
= (a-e;) cb! 
= a,b’c = aje! - exb*c = (a-b)c (6.57) 


Since both f (b)-¢ and (a - 6) c do not depend upon the selection of the basis e;, then (a - e;) e,0y,, 
also cannot depend upon the selection of the basis and we can define 


Of(a) = (a €;) ef, (6.58) 


so that 
Aria (F (b) -¢) = (abe. (6.59) 


From equation 6.58 we see that Ogq) is a vector operator and it obeys the product rule. 


Now consider Of(q) (f (b A c) B) where B = by A by is a bivector. Then 


Ora) (£ (b A c) B) = Ogcay ((£ (6) Af (e)) - (01 A b2)) 
= di.) ((é (b) Af (0)) (by A by) + (f (b) A (0) (by A b2)) 
= dia ((é (b) Af (c)) (i Rb) = (f (c) Af (6)) (by A b2)) . (6.60) 


but by equation B.10 (Appendix B) 
(f(b) A £ (c)) - (b1 A bz) = (£ (0) - ba) (£ (©) - b1) — (£ (0) - br) (£ (©) - ba) (6.61) 


so that (also using equation B.2) 


( 2) 
= (a-b)£f (c) - (b) Aba) = (ab) (£ (c) - B). (6.62) 
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Thus equation 6.60 becomes 


) 
= ((a-6)£ (c) — (ac) f(0))- B 
=f((a-b)c—(a-c)b)-B 
f(a-(bAc))-B. (6.63) 
In general if Az and By are grade 2 multivectors then by linearity 
Os(a) (£ (Az) Bo) = f (a: Aa) - Bo. (6.64) 


The general case can be proved using grade analysis. First consider the following (where the 
subscript indicates the grade of the multivector and C is a general multivector) 


(AyC) = (Ap (Co + +++ + Ch) (6.65) 


then the lowest grade of the general product term A,C, is |p — q| which is only zero (a scalar) if 
p=q. Thus 
(ApC) = (ApCp) = Ap + Cp (6.66) 
Thus we may write by applying equation 6.52 (C' is the multivector (f (a2)... f (a,)) B,) 
(f (ai) A... Af (ar) Br) = ((£ (ai). (a)) Br) 
= f (a,) - (f (a2)... f (a,) B,), (6.67) 


so that 


- (6.68) 
Thus (using equation A.5 in Appendix A) 


r 


Orca) ((£ (a1) A... Af (ar)) Br) = J (-2)* E((a- aia A... Aa A... Ady) Br); 


i=1 


=(6((en ea) OLA A Aa] a.) 


= (I(G (Gi WNicdadWaie) ) Bali: (6.69) 
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Using linearity the general case is 


Orta) (f (A) B) = $7 (f(a (A),) (B),)1 - (6.70) 


Tr 


For a fixed r-grade multivector A, we can write 


‘ : OX tr 
(f (A,.) X;) b= (f (A,.) satin) elitr) 
~ (f (A,) ei) elitr} 
= (F (An) -€¢5) 9 = £4). (6.71 


Applying equation 6.70 to equation 6.71 gives (let f(a-A,) = Bi, eh A... Ae™ bea 
coordinant expansion for f (a- A,)) 


Dewayf (Av) = Orta) (£ (Ar) Xr) Ox, 
= (t(0- An) dx, 
SB gen (OP Nace he) Nes NGM we, (6.72) 
From what we know about subspaces the coefficients of B;,_;,_, are zero unless the e;, \...Ae;, 


contain the same vectors (possibly in a different order) as the e’-! A... Ae” plus one additional 
basis vector, e;. Thus (let E™~' = e*-1 A... Ae" and E,_1 = e;,_, A... A ei) 


(et Aes Net) (eh Neg et A ut We = (Els (Bai Ne) fe Ae") > “(6:78) 


The l.h.s of equation 6.73 can also be put in the form of the r.h.s. since even if the order of 
the e;,’s are scrambled they can alway be put into the same order as the e;,’s via transposition. 
If the order of the e/*’s are made to be the reverse of the e;,’s any minus signs generated will 
cancell out. Also the E,_; and the E’~! can be exchanged since they are scalar multiples of one 
another. The total number of possible e,’s is m — r + 1 where n is the dimension of the vector 
space. Now reduce the scalar coefficient of the blade in equation 6.73 


(E"-)- (E,-1 Ae;)) = (B"* (Er1 Ae;)), 
= (E"" (E,-1e; — Ep-1- e;)) 


=e; — (Era (E"™ -;)), 
Sey (6.74) 
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Now reduce 


ej (7? AE" ") =e;-(P9 AE") +e;A(P AE") 


: y+ (e A E"*) 
= (e;-e’) ate S(- hee) er hy Neth Ae" 
wt oe (6.75) 
Thus 
Ofayf (A,r) = (n-—r+1)f (a: A,). (6.76) 
Now let A, = J then 
f(a) = Op(aydet (f) 1 =f (a-J), (6.77) 


but by equation 1.94 we have 


det (f) £7" (a) I? = f (al") 
det (f) £~' (a) I = f (al) 
( 
( 


=fladi) 
det (f) f~' (a) I =f (a- I) (6.78) 
Orca) det (f) = det (£) £7 (a). (6.79) 


Equation 6.69 can also be used to calculate the functional derivative of the adjoint. The result 
for r > 1 is 


Orta (Ap) = Oe(a) (f (x) A,) Ax, 
a (a: x Ae A,Ox,.. (6.80) 


Equation 6.80 cannot be used when r = 1 since f (a- X,) is not defined. Let A, = 6 then using 
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components we have 
ecayf (b) = ajeiOy,,€* fib! 
= ajeiOy,,€* fucd! 
= a;€;6,5, eb! 


= a,ee!b' = ba. (6.81) 


Chapter 7 


Multilinear Functions (Tensors) 


A multivector multilinear function’ is a multivector function T (Aj,...,A,) that is linear in each 
of its arguments? The tensor could be non-linearly dependent on a set of additional arguments 
such as the position coordinates x’ in the case of a tensor field defined on a manifold. If x denotes 
the coordinate tuple for a manifold we denote the dependence of T on x by T'(Aj,..., Ap; 2). 


T is a tensor of degree r if each variable A; € V, (V,, is an n-dimensional vector space). More 
generally if each A; € G(V,,) (the geometric algebra of V,,), we call T an extensor of degree-r on 


G (Vn). 


If the values of T (a1,...,a,) (a; € Vn V1 <j <r) are s-vectors (pure grade s multivectors in 
G (Vn)) we say that T has grade s and rank r+ s. A tensor of grade zero is called a multilinear 
form. 


In the normal definition of tensors as multilinear functions the tensor is defined as a multilinear 
mapping 


T: Xn > ®, 


i=1 


so that the standard tensor definition is an example of a grade zero degree/rank r tensor in our 
definition. 


‘We are following the treatment of Tensors in section 3-10 of [4]. 
?We assume that the arguments are elements of a vector space or more generally a geometric algebra so that 
the concept of linearity is meaningful. 
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7.1 Algebraic Operations 


The properties of tensors are (a € R, a;,b € V,, T and S are tensors of rank r, and o is any 
multivector multiplicative operation) 


PGi cing aes Op) SO) (isin Ops ey Ox (7.1) 
Oats cig Og Op og Op aE (ages 4 Car ayn Oe) (Oi os es et Oy Qiao Oy) 5 : 
(Pee) (iss: sin) Hl (Cay bien) ES (Gases oye)’: (7.3) 


Now let TJ be of rank r and S of rank s then the product of the two tensors is 


CP O58) (Gig cere Od a EG d cop) OS (Gags ies Gp (7.4) 


O99 


where “o” is any multivector multiplicative operation. 


7.2  Covariant, Contravariant, and Mixed Representations 


The arguments (vectors) of the multilinear fuction can be represented in terms of the basis vectors 
or the reciprocal basis vectors® 


a; =a'J€;,, (7.5) 

=Aj, e : (7.6) 

Equation (7.5) gives a; in terms of the basis vectors and eq (7.6) in terms of the reciprocal basis 
vectors. The index 7 refers to the argument slot and the indices 7; the components of the vector 


in terms of the basis. The Einstein summation convention is used throughout. The covariant 
representation of the tensor is defined by 


de =P Cue ee) (7.7) 
DP (Gay xii) =P (Gps ncaa’, 
=T (Ei, eects €;,) q'! 8S a’ 
=T;, in q'! eee a’. (7.8) 


3When the a; vectors are expanded in terms of a basis we need a notatation that lets one determine which 
vector argument, 7, the scalar components are associated with. Thus when we expand the vector in terms of the 
basis we write a; = a‘i e;,; with the Einstein summation convention applied over the i; indices. In the expansion 
the j in the a’ determines which argument in the tensor function the a” coefficients are associated with. Thus 
it is always the subscript of the component super or subscript that determines the argument the coefficient is 
associated with. 


7.2. COVARIANT, CONTRAVARIANT, AND MIXED REPRESENTATIONS 


Likewise for the contravariant representation 


T (a1, G.) =P (aj,e Ces ,a;,€'") 
=T ( as e’") Gigpa das 
—T tr Ay 
- 


One could also have a mixed representation 


Toe =f =T (a Clrpa neil Csi (OC eraee) 
= (€:,,...,€:,,e°41,... ,e a Oss Ai,. 
=T; i. tet tr yt a eer nears 
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(7.9) 


(7.10) 


(7.11) 


(7.12) 


In the representation of T one could have any combination of covariant (lower) and contravariant 


(upper) indices. 


To convert a covariant index to a contravariant index simply consider 


P(E ipweyO) coe y ie) =f ON ne kc 5€p,,.--,€i,) 
= BE ( Crs Perea ity Cin) 
ee + 3 
cee ws phe Up =g" PD gs pes 


Similarly one could raise a lower index with gj x,. 


(7.13) 
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7.3 Contraction 


The contraction of a tensor between the j*” and k*” variables (slots) is* 
T (ai, +++, Qj-1, Wins Aj+is++- Oe) = Va; (Va, lL (a1, es st) : (7.14) 


This operation reduces the rank of the tensor by two. This definition gives the standard results 
for metric contraction which is proved as follows for a rank r grade zero tensor (the circumflex 
““” indicates that a term is to be deleted from the product). 


DCist sf sith) =q" ie eee (7.15) 
Val Seg... (0.1, a’) get De aie 
eli) a"" age tess hal Lees (7.16) 


Was! : (V.;T) —ekm ‘ els5/a" hk ati om (Agem ai”) i aT, 
agimigpidimalt als a aT 
=gimigh a ats is qin eo aT. 
=giimg® a4... a'm aT; 
= (G8 Thy ig cin ning) GO. a 7.17) 


Equation (7.17) is the correct formula for the metric contraction of a tensor. 


If we have a mixed representation of a tensor, Ty eles and wish to contract between an 
upper and lower index (i; and 7,) first lower the upper index and then use eq (7.17) to contract 
the result. Remember lowering the index does not change the tensor, only the representation of 
the tensor, while contraction results in a new tensor. First lower index 

Fe 45 a Lower Index Gis; T, kj (7.18) 


te eet 


“The notation of the l.h.s. of eq (7.14) is new and is defined by Va, = e!*0,1, and (the assumption of the 
notation is that the 0,., can be factored out of the argument like a simple scalar) 


T (aj, ---,@j-1, Vag, Aj41;-++54r) =T (Oi ated ip e Orgy Oia c5t eis ta Ge) 
=T (Gin Sp OHO PO a Opis sO ies Hae) 
Se ca eT (Gii.5 of Opi Oy 0741) 14 Bigs. 84) 
SPOn OPE (dys sf Gh pies Gris. ples ety Oe) 
=e eT (ag; <6 O51, Cj Ojeiys- 64 Cig +s Or) 
St aks bara Tee Pe ee a wees ate att a’, 
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Now contract between 7; and 7, and use the properties of the metric tensor. 


kj Contract ijtk 5 
Gigkj Li... wedfgesby 7g? Gigkj Ly... wedfgeesdp 
é ks 
oo Stk J 
pt tags (7.19) 


Equation (7.19) is the standard formula for contraction between upper and lower indices of a 
mixed tensor. 


7.4 Differentiation 


If T (a1,...,@,;2) is a tensor field (a function of the position vector, x, for a vector space or the 
coordinate tuple, x, for a manifold) the tensor directional derivative is defined as 


DY Gos j 08 2) = (Orga VE (Wie es (7.20) 
assuming the a’) coefficients are not a function of the coordinates. 


This gives for a grade zero rank r tensor 


(Qr44 : V) T (a1, ar , Oy) =a S03 q'! aes a Se 


=O es Ha Oa sy 


(7.21) 


Lindy? 


7.5 From Vector/Multivector to Tensor 


A rank one tensor corresponds to a vector since it satisfies all the axioms for a vector space, 
but a vector in not necessarily a tensor since not all vectors are multilinear (actually in the case 
of vectors a linear function) functions. However, there is a simple isomorphism between vectors 
and rank one tensors defined by the mapping v (a) : V > ¥® such that if v,aeV 


v(a) =v-a. (7.22) 


So that if v = v’e; = vje' the covariant and contravariant representations of v are (using e'-e; = 
3) — 

O(a) =e" =] way. (7.23) 
The equivalent mapping from a pure r-grade multivector A to a rank-r tensor A (a,,...,a;) is 


A Gia Ge)'= As (ap Aasn ap) (7.24) 
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Note that since the sum of two tensor of different ranks is not defined we cannot represent a 
spinor with tensors. Additionally, even if we allowed for the summation of tensors of different 
ranks we would also have to redefine the tensor product to have the properties of the geometric 
wedge product. Likewise, multivectors can only represent completely antisymmetric tensors of 
rank less than or equal to the dimension of the base vector space. 


7.6 Parallel Transport Definition and Example 


The defintion of parallel transport is that if a and 6 are tangent vectors in the tangent spaced of 
the manifold then 
(a-Ve)b=0 (7.25) 


if b is parallel transported in the direction of a (infinitesimal parallel transport). Since b = b’e; 
and the derivatives of e; are functions of the x’’s then the b’’s are also functions of the x’’s so 
that in order for eq (7.25) to be satisfied we have 


(a- Vz) b =a'd,: (b’e;) 
=a' ((d,:b’) e; + b/A,.€;) 
=a! ((d,ib’) e; + UTE ex) 
=a! ((Bqb!) e) + 0°T.e)) 
=a' ((0,:b7) + b'T4,) ej = 0. (7.26) 


Thus for 6 to be parallel transported (infinitesimal parallel transport in any direction a) we must 
have 
d,ib) = —BFTY,. (7.27) 


The geometric meaning of parallel transport is that for an infinitesimal rotation and dilation of 
the basis vectors (cause by infinitesimal changes in the x’’s) the direction and magnitude of the 
vector b does not change to first order. 


If we apply eq (7.27) along a parametric curve defined by 2x (s) we have 


ab) de Ob! 
ds ds Ox' . 
usp (7.28) 


ds 
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and if we define the initial conditions b/ (0) e;. Then eq (7.28) is a system of first order linear 
differential equations with intial conditions and the solution, b’ (s) e;, is the parallel transport 
of the vector b’ (0) e;. 


An equivalent formulation for the parallel transport equation is to let y(s) be a parametric curve 
in the manifold defined by the tuple y(s) = (z'(s),...,2"(s)). Then the tangent to 7(s) is 
given by 


dy — dx’ 
and if v (a) is a vector field on the manifold then 
dy ASO 5 28 
& v.) ods Oxi c ei) 
_ dx! (dv! - 08; 
~ ds \art 9  ” Oxi 
dx’ (dv! a 
~ ds & a ether) 
dx' Ov4 dx 
~ ds Oat? | d UT ee 
dw dz’ , 
= — — TY? 
( ds ds | = 
=0. (7.30) 


Thus eq (7.30) is equivalent to eq (7.28) and parallel transport of a vector field along a curve is 
equivalent to the directional derivative of the vector field in the direction of the tangent to the 
curve being zero. 


As a specific example of parallel transport consider a spherical manifold with a series of concentric 
circular curves and paralle transport a vector along each curve. Note that the circular curves are 
defined by 


s 
u(s) =Uo + a cos (=—) 


v(s) =p +asin (—) 
ray 27a 
where u and v are the manifold coordinates. The spherical manifold is defined by 
x =cos (wu) cos (v) 
y =cos (u) sin (v) 


z =sin (u). 
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Note that due to the dependence of the metric on the coordinates circles do not necessarily appear 
to be circular in the plots depending on the values of uo and vo (see fig 7.2). For symmetrical 
circles we have fig 7.1 and for asymmetrical circles we have fig 7.2. Note that the appearance of 
the transported (black) vectors is an optical illusion due to the projection. If the sphere were 
rotated we would see that the transported vectors are in the tangent space of the manifold. 


Figure 7.1: Parallel transport for ug = 0 and vp = 0. Red vectors are tangents to circular curves 
and black vectors are the vectors being transported. 


If y(s) = (u(s),vu(s)) defines the transport curve then 


ey (7.31) 
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Figure 7.2: Parallel transport for uo = 7/4 and vp = 7/4. Red vectors are tangents to circular 
curves and black vectors are the vectors being transported. 


and the transport equations are 
dy — (dudf* dvodf" du... 
(Fv) F=(F a + 7s Bo — sin (u) cos (u) Ff e,+ 


dudf®  dvof®  cos(u) fdu,, . dv ., 
(5 Ou eee Ov Y sin (u) ( Mt asf ))e 


-(< — sin (w) cos (w) =f") e,+ 


ds 
(ar Sey (Grr) ‘a 
oe eos (2 (7.33) 


df’ _—_ cos(u) /du ,, ; 
ds sin (u) ( . i ) on 
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If the tensor component representation is contra-variant (superscripts instead of subscripts) we 
must use the covariant component representation of the vector arguements of the tensor, a = a;e’. 
Then the definition of parallel transport gives 


(a- V,,) b =a'd,i (bje’) 
=q! ((O,b;) e! + b;0,.e") : (7.35) 


and we need 
(O,ib;) e7 + bjO,ve? = 0. (7.36) 


To satisfy equation (7.36) consider the following 


0,: (e7 - ex) =0 
(d,:e7) - e, + e? - (O,e,) =0 
(0,:e7) -e, +e? -e1%, =0 
(0,.e7) - Ee, + 6/r', =0 
(ane?) en +14, =0 

(O,0e7) +e, =— ge (7.37) 


Now dot eq (7.36) into e; giving? 
(O,.bj) €7 « ex +b; (Ope?) » ex =0 


(O,ib;) 6% — bE, =0 
(O,ibp) = b,T%,. (7.39) 


7.7 Covariant Derivative of Tensors 


The covariant derivative of a tensor field T (a,,...,a,;x) (x is the coordinate tuple of which T 
can be a non-linear function) in the direction a, is (remember a; = aki ex, and the e;, can be 
functions of x) the directional derivative of T (a,,...,a,;x) where all the a; vector arguments of 
T are parallel transported. 


>These equations also show that 
d,ve? = -Tie*. (7.38) 
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Thus if we have a mixed representation of a tensor 


DA Oiigeus Og a Tes feras et png 2.6 Cap aistedag (7.40) 
the covariant derivative of the tensor is 
(ie = DOLE (Gites Gee) = Tite of NG Caine 
. 3 Bas 41.105 ge qr a" Aiesy a; oor 
= 
a5 ; oe a SG Gale thks, Oye 4,07 
q=st+ 
=e 7 @ ee: a" di, bike Guar 


¢ 
lg isp etgetr al a a +1 
fe y Ee 4, OP ai Oe ae oO aslo: (7.41) 


q=st+l 
From eq (7.41) we obtain the components of the covariant derivative to be 


OT Is1..tr s r 
Tats : ip Ust1.tr : lq istl-vdgq--tr 
Oxrrt > eee a aera e men : (7.42) 


p=1 q=st+l 


To extend the covariant derivative to tensors with multivector values in the tangent space (geo- 
metric algebra of the tangent space) we start with the coordinate free definition of the covariant 
derivative of a conventional tensor using the following notation. Let T (a,,...,a,;2) be a con- 
ventional tensor then the directional covariant derivative is 


(DP Se 0 OVI Ge sese)— >. TGys 5 (0° V) oj yes2): (7.43) 
j=l 


The first term on the r.h.s. of eq (7.44) is the directional derivative of T if we assume that the 
component coefficients of each of the a; does not change if the coordinate tuple changes. The 
remaining terms in eq (7.44) insure that for the totality of eq (7.44) the directional derivative 
(b- V)T is the same as that when all the a; vectors are parallel transported. If in eq (7.44) 
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we let b- V be the directional derivative for a multivector field we have generalized the defin- 
intion of covariant derivative to include the cases where T'(a,,...,@,;x) is a multivector and 
not only a scalar. Basically in eq (7.44) the terms T (e;,,...,¢;,.;x2) are multivector fields and 
(b- V)T (ei,,..., 6:32) is the direction derivative of each of the multivector fields that make up 
the component representation of the multivector tensor. The remaining terms in eq (7.44) take 
into account that for parallel transport of the a,’s the coefficients a‘ are implicit functions of the 
coordinates «*. If we define the symbol V, to only refer to taking the geometric derivative with 
respect to an explicit dependence on the x coordinate tuple we can recast eq (7.44) into 


(b- D)T = (b- Va) T (a1,...,4r50) — ST (a1,...,(b V) aj,...,ar3 2). (7.44) 


7.8 Coefficient Transformation Under Change of Variable 


In the previous sections on tensors a transformation of coordinate tuples Z (x) = (7 (x) ,...,Z" (2)), 
where x = (x',..., 2"), is not mentioned since the definition of a tensor as a multilinear function 
is invariant to the representation of the vectors (coordinate system). From our tensor definitions 


the effect of a coordinate transformation on the tensor components is simply calculated. 


If R(x) = R(z) is the defining vector function for a vector manifold (R is in the embedding 
space of the manifold) then® 


OR Ox! 
aoa Ogi ee) 
_— OR Ox 


econ, <e! bot au 
°For an abstract manifold the equation é; = age can be used as an defining relationship. 
zt 
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Thus we have 


T (e;,, tee €;,) =Tiy.ip 
T (;,, ee se; =T yy .5p 
OT) _ OXI _ 
T (e;,, derins €;,) =T (Fre eras ores 
OT ORI : 
Ona ssias Age (€;,, eens é;,) 
Ont. “Og = 


Equation (7.51) is the standard formula for the transformation of tensor components. 


LOT 


Chapter 8 


Lagrangian and Hamiltonian Methods! 


8.1 Lagrangian Theory for Discrete Systems 


8.1.1 The Euler-Lagrange Equations 


Let a system be described by multivector variables X;, 7 = 1,...,m. The Lagrangian L is a 
scalar valued function of the X;, X; (here the dot refers to the time derivative), and possibly the 
time, t. The action for the system, S, over a time interval is given by the integral 


te ; 
S= [ dtL (eat) (8.1) 


The statement of the principal of least action is that the variation of the action 6S = 0. The 
rigorous definition of 6S = 0 is let 


X}(t) =X, () + % (0) (8.2) 


where Y;(t) is an arbitrary differentiable multivector function of time except that Yj (ti) = 
Y; (t2) = 0. Then 


dS 
S= | 0 (8.3) 


'This chapter follows “A Multivector Derivative Approach to Lagrangian Field Theory,” by A. Lasenby, C. 
Doran, and S. Gull, Feb. 9, 1993 available at http://www.mrao.cam.ac.uk/ cjld1/pages/publications.htm 


129 


130 CHAPTER 8. LAGRANGIAN AND HAMILTONIAN METHODS 


Then 
L (xi, X1,t) 25 (x, rs ee Cee i,t) (8.4) 
se (Xi, %,t) +e) (Vit Ox,L + Yi +dx,L) (8.5) 
coe [ dt [2 (X;, %,,t) es 2 (vi « Ox, L + Yj x 01)) (8.6) 
< ae - ay (vi « Ox, L + Yj x Ox,L) (8.7) 


te m d 
= [ oes . (ax.t a (x,t) (8.8) 


where we use the definition of the multivector derivative to go from equation 8.4 to equation 8.5 
and then use integration by parts with respect to time to go from equation 8.7 to equation 8.8. 
Since in equation 8.8 the Y;’s are arbitrary 6S = 0 implies that the Lagrangian equations of 
motion are 


d 
a 
The multivector derivative insures that there are as many equations as there are grades present 
in the X;, which implies there are the same number of equations as there are degrees of freedom 
in the system. 


dx,L— = (dx,L) =0, Vi=1,...,m. (8.9) 


8.1.2 Symmetries and Conservation Laws 


Consider a scalar parametrised transformation of the dynamical variables 


A= XG @), (8.10) 
where Xj; (X;,0) = X;. Now define 
dX} 
5X, = St (8.11) 
da a=0 


and a transformed Lagrangian 


LU (Xi, X:,t) =L (xi, Xi,t) (8.12) 
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Then 


dL! 
da a=0 


3 (AX; + Oxi! + IX; « Dy’) 


L$: (aoe £ (:-0et) se ¥ (ut) 


- 3 (ax * (a2 - < (2.,2") + < (ax; « a,,L") : (8.13) 


i=1 


If the X;’s satisfy equation 8.9 equation 8.13 can be rewritten as 


ay! ee 
— = — (OX; Ox. Lt) 8.14 
da|, a 2 * ere) 
Noether’s theorem is - 
dL 
— ak => x (6X; * Ox, L) = conserved quantity (8.15) 
da |, a 


From D& L - 


“If the transformation is a symmetry of the Lagrangian, then L’ is independent of a. In this 
case we immediately establish that a conjugate quantity is conserved. That is, symmetries of the 
Lagrangian produce conjugate conserved quantities. This is Noether’s theorem, and it is valuable 
for extracting conserved quantities from dynamical systems. The fact that the derivation of 
equation 8.14 assumed the equations of motion were satisfied means that the quantity is conserved 
‘on-shell’. Some symmetries can also be extended ‘off-shell’, which becomes an important issue 
in quantum and super symmetric systems.” 


A more general treatment of symmetries and conservation is possible if we do not limit ourselves 
to a scalar parametrization. Instead let X/ = X;(X;,/) where M is a multivector parameter. 
Then let 

Poe (xi, x t) (8.16) 


and calculate the multivector derivative of L’ with respect to M using the chain rule (summation 
convention for repeated indices) first noting that 


£ (Aw du) X a xb") = ((A* du) X1) « eel! + (Ax Ou) Xf) *& (OL!) (8.17) 
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and then calculating 
(Ax Oy) L! = ((A« Dy) X1) Ox L! + ((A* Ou) X1) # Og_L! 


= ((A* Oy) X!) * (a0! 7 c (2,,2") ae & (((A* Ou) X1) ¥ xb") (8.18) 


If we assume that the X}’s satisfy the equations of motion we have 


(Aso = c ((A EO) XY» xf) (8.19) 


and differentiating equation 8.19 with respect to A (use equation 6.43) gives 


OyL' = c (a1 ((A * Oy) X!) xb) 
a 


~ at 


Equation 8.20 is a generalization of Noether’s theorem since if 0,,L' = 0 then the parametrization 
M is asymmetry of the Lagrangian and (OyX/) * Ox,L’ is a conserved quantity. 


((Ax1X!) x xb) (8.20) 


8.1.3 Examples of Lagrangian Symmetries 
Time Translation 


Consider the symmetry of time translation 


Xj (t,a) = X; (t+ a) (8.21) 
so that eo 
x,=—*| =X, (8.22) 
da a=0 
and 
dL ters 
— oo X,;* Oy L 8.23 
oye] per 7 *0xL) ( ) 


0= “ (>: (x x ax,b) 2 1) (8.24) 
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The conserved quantity is the Hamiltonian 
H=)> (x; « Ox,1) ae fy (8.25) 
i=1 


In terms of the generalized momenta 
P, =, L, (8.26) 


so that 


He ss (%, * P,) ae (8.27) 


Central Forces 


Let the Lagrangian variables be x; the vector position of the i‘” particle in an ensemble of N 
particles with a Lagrangian of the form 


N i N WN 
L= > matt — >> Vii (loi — 25) (8.28) 
i=1 


i=l j<i 


which represent a classical system with central forces between each pair of particles. 


First consider a translational invariance so that 


v= 2x, +c (8.29) 
where q is a scalar parameter and c is a constant vector. Then 
(1 =c (8.30) 
and 
= EG (8.31) 


so that the conserved quantity is (equation 8.15) 


N N 
i=l a 
c:p=C: So miti 
i=1 


N 
p= >_ mis (8.32) 
1=1 
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since c is an arbitrary vector the vector p is also conserved and is the linear momentum of the 
system. 


Now consider a rotational invariance where 


! ecB/2 


t= were (8.33) 


where B is an arbitrary normalized (B? = —1) bivector in 3-dimensions and a is the scalar angle 
of rotation. Then again L’ = L since rotations leave #? and |x; — «;| unchanged and 


a = ; (BetBl2,e-28/2 = eB /2 eB)? B) 
a 
1 


Remember that since B-2; is a vector and the scalar product (*) of two vectors is the dot product 
we have for a conserved quantity 


: (B- x) - (0;,L) = ym (B- 2,) - 4; (8.35) 
=B-S mi (xi A %i) (8.36) 

=B. i (8.37) 

J= So mi (xi A 4i) (8.38) 


where we go from equation 8.35 to equation 8.36 by using the identity in equation B.12. Then 
since equation 8.35 is conserved for any bivector B, the angular momentum bivector, J, of the 
system is conserved. 


8.2 Lagrangian Theory for Continuous Systems 


For ease of notation we define 
AV = AV. (8.39) 
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This is done for the situation that we are left differentiating a group of symbols. For example 
consider 


(ABC) V = (ABC)V. (8.40) 


The r.h.s. of equation 8.40 could be ambiguous in that could the overdot only apply to the B 
variable. Thus we will use the convention of equation 8.39 to denote differentiation of the group 
immediately to the left of the derivative. Another convention we could use to denote the same 
operation is 


ABCY = (ABC) ¥V (8.41) 


since using the “hat” symbol is unambiguous with respect to what symbols we are applying the 
differentiation operator to since the “hat” can extend over all the relevant symbols. 


8.2.1 The Euler Lagrange Equations 


Let w; (x) be a set of multivector fields and assume the Lagrangian density, £, is a scalar function 
L(wv;, Vvi, x) so that the action, S', of the continuous system is given by 


4 


where V is a compact n-dimensional volume. The equations of motion are given by minimizing 
S using the standard method of the calculus of variations where we define 


5 (@) = vi (x) + Gi (x), (8.43) 


and assume that 7; (x) yields an extrema of S and that ¢; (x) = 0 for all x € OV. Then to get 
an extrema we need to define 


S(6.) = [ |dx”| L (wh; + €b;, Vd; + Vi, 2) , (8.44) 


Os Os 
so that S (0) is an extrema if ae 0. Let us evaluate De (summation convention) 
€ € 


Os 


Oc 


i ldx"| (6: * Op,) £+ (Vb; * Ovy,) £). (8.45) 


e=0 
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Start by reducing the second term in the parenthesis on the r.h.s. of equation 8.45 using RR5 
(Appendix C) 


= (40r0h~ (6 (2) 

=V-(GOeu£), — (4: (Oey,£) ¥) 

=V-(¢iOeu£), — 6: * ((Orn,£) V). (8.46) 
So that equation 8.45 becomes 


as - [ da"| 6; (84,£ — (Ovg.l) ¥) + i lde"|V - (6iOve,L), 


7 dic] i * * (OL — (Ovy,L) ¥) av / Jas" "| n- (6i0ou,L) 
V av 


The hee |dx”| V - (@;0vy,£), term is found to be zero by using the generalized divergence theorem 
(equation 4.84) and the fact that the ¢;’s are zero on OV. If ¢; is a pure r-grade multivector we 
have by the properties of the scalar product the following Lagrangian field equations 


(Oy£ — (Ovy,L) MJ: =0 (8.48) 


7 (23 ae (Aevy'£) ). = (8.49) 


For the more general case of a ¢; being a mixed grade multivector the Lagrangian field equations 
are 


Oy L i (Ovu;L) V =0 (8.50) 


or 


ayL-V (Avuy£) =) (8.51) 


Note that since equation 8.50 is true for ~; being any kind of multivector field we have derived 
the field equations for vectors, tensors (antisymmetric), spinors, or any combination thereof. 
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8.2.2 Symmetries and Conservation Laws 


We proceed as in section 8.1.2 and let wi = wi (w;,. M) where M is a multivector parameter. Then 
L' (hi, Voi) = Li, V4) (8.52) 
and using equation 6.39 and RR5 


(A* Ou) L! = ((A* Ou) Ui) * (OyL’) + ((A * Ou) VU!) & (OoyL’) 
= ((A * Ou) Wj) * (OyL') + (A * Ou) Vedi) (Ovy:L’)) 
((A * yr) U4) * (By.L!) + (V (((A* ur) YD) (AvwiL’)) — ¥ (((A * Om) V1) (GeueL’) J ) 
= ((A* nr) Wt) * (yl!) +V- (((A* Om) W) (Ooy£)), - a ) (Oey£') 7) 
= ((A* du) Ui) * (QyL!) +0 - (((A* Ou) V4) (Oow£)), — (A * Our) Wt) * (Bowl?) ¥ 
= (A Ou) Hi) * ((BuL’) - (00uj0) 8) +9-(( (8.53) 


and if the Euler-Lagrange equations are satisfied we have 
(A * Ou) L' = V-(((A* Om) ¥4) (Ovw£')), (8.54) 


and by equation 6.43 
Ol! = Aa (V-(((A* An) YD) (Oowl’)),) (8.55) 


If OL’ = 0, equation 8.55 is the most general form of Noether’s theorem for the scalar valued 
multivector Lagrangian density. 


If in equation 8.55 M is a scalar a and B is a scalar 3 we have 
ant! = a3 (V- ee ) vt) (Ovyi£’)),) 
=V-((davi) (AvwL£’)), - (8.56) 
Thus if a = 0 in equation 8.56 we have 


Bok |p = V > (Oath) (Ovu£’)), (8.57) 


a=0 


which corresponds to an differential transformation (0,£' = 0 is a global transformation). If 
O.L'|,-9 = 0 the conserved current is 


j= (Gath) (Brw£)),| (8.58) 
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with conservation law 


Vg=t. (8.59) 
If O.L'|,-9 # 0 we have by the chain rule that (g(x) = Og2"|,_9) 
Ope 3 One |e VE = oe We (8.60) 
and consider 
V- (gL) =(V-g)L+g9-VE. (8.61) 
If V-g =0 we can write O,L'|,. as a divergence so that 
V > (((Oatbile=o) (Ovu£)), — (Aat"|gao) £) = 0 (8.62) 
and 
J = ((Oajla—o) (AvuiL)); — (A240) £ (8.63) 


is a conserved current if V- 0,2"|,-9 = 0. 


Note that since 0,2'|,_, is the derivative of a vector with respect to a scalar it itself is a vector. 
Thus the conserved 7 is always a vector that could be a linear function of a vector, bivector, etc. 
depending upon the type of transformation (vector for affine and bivector for rotation). However, 
the conserved quantity of interest may be other than a vector such as the stress-energy tensor or 
the angular momentum bivector. In these cases the conserved vector current must be transformed 
to the conserved quantity of interest via the general adjoint transformation A* 7 (B) = Bj (A). 


8.2.3. Space-Time Transformations and their Conjugate Tensors 


The canonical stress-energy tensor is the current associated with the symmetries of space-time 
translations. As a function of the parameter a we have 


x =xt+an (8.64) 
Wi (x) = Yi (e + an) (8.65) 
Then 
OyL'| 9 = Ok (4 +an)|,-9 =n: VL=V- (nl) (8.66) 
Ooi lgao = 2° Vi (8.67) 


and equation 8.56 becomes 


V (nL) =V-((n- Vedi) (Ovu,£)), (8.68) 


8.2. LAGRANGIAN THEORY FOR CONTINUOUS SYSTEMS 


so that 


V -T(n) =V-((n- Vv) (Ovy,£) — nL), =0. 
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(8.69) 


Thus the conserved current, 7’ (n), is a linear vector function of a vector n, a tensor of rank 2. 
In order put the stress-energy tensor into the standard form we need the adjoint, T (n), of T (n) 


(we are using that the adjoint of the adjoint is the original linear transformation). 


T (n) 


((n- V¥i) (Avy£))1 — n£ 
= (n : v) CP (Dey£)). — nl. 


Using equation 6.34 we get 


or in rectangular coordinates 


1p a) a) Ory 
= His (yee HY 
: (v) : (« x) eae ie) dae 


Thus in standard tensor notation 


OTH’ 
=0 
Ox : 


so that T (n) is a conserved tensor. 


(8.70) 
(8.71) 


(8.72) 


(8.73) 


(8.74) 


(8.75) 


(8.76) 


Now consider rotational transformations and for now assume the w; transform as vectors so that 


, _oB  aB 
xy =e€E 2 ve 2 
Ui (a) =e yj (2!)e"™, 


(8.77) 
(8.78) 
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Note that we are considering active transformations. The transformation of « to x’ maps one 
distinct point in space-time to another distinct point in space-time. We are not considering 
passive transformations which are only a transformation of coordinates and the position vector 
does not move in space-time. Because the transformation is active the sense of rotation of the 
vector field 7 (x) is opposite that of the rotation of the space-time position vector”. Thus 


Br = so (cB —Bzr)e® 
On |p — eB, (8.79) 
dail = Ze (Bui (0") — Wi (2) BE + ef Oadi(a)) oF 
=e (Bx di (a) eF +eF (gat! Vay (2!) eF (8.80) 
Aabllano = BX Yi (#) + Bait'|aao Vi (2) 
=Bxy,+(e¢-B)-Vi; =: B+ (a: B)-Vui. (8.81) 


Thus 


V- (do2"|,-9) = V+ (a: B)=—-(B-4)-V=-B- (#av) =B-(VAz)=0 — (8.82) 
since V A x = 0° the derivative of the transformed Lagrangian at a = 0 is a pure divergence, 


OSV (GBD): (8.83) 


=V-(Bxy-(B-2)- Vi) (Ovu,L£)), (8.84) 


a=0 


V (derbi) (OvwL)), 


?Consider an observer at location x that is rotated to location 2’. If he is rotated in a clockwise sense he 
observes the vector field to be rotating in a counter clockwise sense. 
3Consider 


0 
VAa=e’>— A aye" 


Ox” 


= v n 
= Dgy ine Ae 


= He’ fell 
= Agu Inet e Ae 


Ox" 


= pv ome’ er 


= Hanne” Ae” 


= gnve”’ Ne" =0 


since gp, is symmetric and e” A e” is antisymmetric. 
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V-F(B) = V+ (Bari) (Ovw£)),| _ — AaL'lano 
=V-(((Bx Ui — (B-2)- Vi) (Bvy,£)), ~ (eB) £) (8.85) 
J(B) = ((Bx Yi — (B+) Vi) (Ovy.L)), — (@- B)L. (8.86) 


By Noether’s theorem V - Ei (B) =0 where J (B) is a conserved vector. So that 
Vira = i(v) -B=O0foralB = (Vv) =0 


The adjoint functon J (n) is, therefore a conserved bivector-valued function of position (we can 
use - instead of * since all the grades match up correctly). 


Now consider the coordinate (tensor) representation of J (v) =0 


a=nye” (8.87) 
DAD) SPOT By IVE, (8.88) 
cies Ojery 
J (v) =e, Ne, = 0 (8.89) 
OJHYy 
aa (8.90) 


To compute J(n) note that A is a bivector using equation 6.34 to calculate the adjoint of the 
adjoint and (A « 0g) B = A® and RR7° we get 


Ax J (n) = (Ax Oz) )(I(B )n) 

= (A* dp) (((B x Ui — (B+) Vda) (Oey,£)), 2 — (w+ B) Ln) 

= ((A x Yj — (A> &) - Vivi) (Ovy£n) — (@ - A) Ln) 

= ((A x hi — (A+ (@ AV)) i) (Ovy,£r) — (@- A) Ln). (8.91) 


Using RR5’ (Reduction Rule 5 in Appendix C) the first term on the r.h.s. of equation 8.91 


for) 


((A)» af b = (A)y- (ab) 
9 andi ao Ve ae 
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reduces to 


(A x Wi) (vu £n)) = 5 (AY: — Wi A) (Ovu,£n)) 


(A {i (Avy, Ln) — (OvyL£n) Vit) 


(A (Wi x (Ovy,£7))) 
Ax (wi; x (Ovu,£7)). (8.92) 


i 
D 
il 

9 


using RR7 the second term reduces to 


(A+ (@AV)) i (Ove,£n)) = (A: (©AV)) (dideu£n) 


ee ((« A v) (Uidey£n) ) (8.93) 
using RR5 again the third term reduces to 
le. Aycay= ; (2A — Az) Ln) 
= ; (A (nx — xn) L) 
=(A(nAx)L) 
=A da AL: (8.94) 


Combining equations 8.92, 8.93, and 8.94 reduces equation 8.91 to 
L(n) = (Wi; x (Ovu,£7n)) 9 — (2 A Vv) (Hidey Ln) —(nAx)£L (8.95) 


where J (n) is the angular momentum bivector for the vector field. Using equation 8.93 we can 


reduce equation 8.95 to 
J(n) =T (n) A+ (iy x (Ovy,£n)),- (8.96) 


If ~; is a spinor instead of a vector the transformation law is 
Vi (a) =e? Wi (2’) (8.97) 
so that - 
Bait =e Si (1!) + €F Oa’ (2') (8.98) 


and B 
Ou lao = 3 Vi fee). 2 Nes (8.99) 
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Equations 8.86 and 8.95 become 


FB) = (Fu: eee Vi) (rm£)) _(e-B)L (8.100) 


J(n) = (s (Avv,£n) - (2 A v) (Uidey,Ln) —(nAz)L (8.101) 
2 2 
where J(n) is the angular momentum bivector for the spinor field. 


Since spinors transforms the same as vectors under translations the expression (equation 8.72) 
for the stress-energy tensor of a spinor field is the same as for a vector or scalar field, T(n) = 


Vv (di (Ovy,L) n) —nZ, so that 


Spinor Field: J(n) =T(n)Axu+t+ (3 (Avy,£n) : (8.102) 
2 
compared to 
Vector Field: J(n) =T'(n) Ax+ (Yi x (Ovu,£n)), - (8.103) 


8.2.4 Case 1 - The Electromagnetic Field 


Example of Lagrangian densities are the electromanetic field and the spinor field for an electron. 
For the electromagnetic field we have 
C=5F-F-A-J (8.104) 
where F = V A A and using eq 6.15 and eq 6.17 
F.F=((VAA)(VAA)) 


=(5(vA ~ (VA) aa ; (vA = vay’) 

z i((v4 A)') ) 

= : VAVA-VA(VA)' — (VA) VA+ (VA)! (VA)*) 
2 
1 
2 


( 
((VA — (VA) «(VA)t —(VA)t « (VA) + (V4)t # (v)t) 
((vA 


— (VA) * (vA)") (8.105) 
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Now calculate using eq 6.47 and eq 6.51 


Oat” = savay ((vA) * (VA) — (VA) * (vA)') 


=(VA)I-VA 
=-2VAA (8.106) 


We also have by eq 6.44 and that AT = A 
Oat (A+ J) = Ogi (J * A) 


= J (8.107) 

so that 
Onl£L—V (Ava£) =JAV(VAA) =O (8.108) 
VF=J (8.109) 


8.2.5 Case 2 - The Dirac Field 


For the Dirac field 


L= (Vb lye! — eAdyuht — may") 
= (Ve) « Ly") — (eAdbyet) — mp «yt (8.110) 
= (Velyz) * U' — (eAdyy!) — mp xy! (8.111) 
The only term in eq 8.110 and eq 8.111 that we cannot immediately differentiate is 0,; (eApyyt). 
To perform this operation let w' = X and use the definition of the scalar directional derivative 
(eA (Xt + ABT) % (X + AB) — eAXty,X) 
h 


(B « Ox) (eAX'y,X) = lim 
= (eABly,X + eAX'y,B) 
= Bt « (ey.X A) + Bx (eAX',) 
=Bx (cAX') + Bx (eAX') 
= Bx (2eAX'y) (8.112) 
Ox (eAX'y,X) = 2AX yy (8.113) 
Oy (eAdyy") = 2eAdy: (8.114) 
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The other multivector derivatives are evaluated using the formulas in section 6.2 


cout (VY) « Fred")) = viel! = del = 12 (8.115) 
Oyt (Vely2) * ve") = Vole (8.116) 
Opt (raeh * p!) = 2my. (8.117) 


The Lagrangian field equation are then 
Os =V (Avuyr£) SWityadekin: tail Viiv: 0 


= 2(Vely, — eA, — mp) = 0 
0= Vuly, — eA — mw. (8.118) 


8.2.6 Case 3 - The Coupled Electromagnetic and Dirac Fields 


For coupled electromagnetic and electron (Dirac) fields the coupled Lagrangian is 


is (ver! Alaa = mans. i) 


1 

= (Ve) + (Fret) — (eArpt) — map pt + oF (8.119) 
1 

= (Vily.) +t — 0A + (Wrept) — mb a dt + SPY, (8.120) 


where A is the 4-vector potential dynamical field and w is the spinor dynamical field. Between 
the previously calculated multivector derivatives for the electromagnetic and Dirac fields the only 
new derivative that needs to be calculated for the Lagrangian field equations is 


Oat (eA * (dybt)) =e (dyebt)! = ew. (8.121) 


The Lagrangian equations for the coupled fields are then 


dyl=V (Avyt£) = 2(Vwly, — eAvy: — my) = 0 (8.122) 
behav (Avay£} = —evywt +V(V AA) =0 (8.123) 

or 
Vly, — eAdry = mb (8.124) 


VE = evry. (8.125) 


Chapter 9 


Lie Groups as Spin Groups! 


9.1 Introduction 


A Lie group, G, is a group that is also a differentiable manifold. So that if g € G and x € RN 
then there is a differentiable mapping ¢ : RY — G' so that for each g € G we can define a tangent 
space Ty. 


A example (the one we are most concerned with) of a Lie group is the group of n x n non- 
singular matrices. The coordinates of the manifold are simply the elements of the matrix so that 
the dimension of the group manifold is n? and the matrix is obviously a continuous differentiable 
function of its coordinates. 


A Lie algebra is a vector space g with a bilinear operator [-,-] : g x g  g called the Lie bracket 
which satifies (x,y,z € g and a, 6 € ®): 


[ax + by, z] =alx, z] + bly, z], (9.1) 
[x, x] = 0, (9.2) 
[z, [y, 2] + [z, [z, yl] + [y, [z, 2] = 0. (9.3) 
Equations (9.1) and (9.2) imply |, y] = —[y, 2], while eq (9.3) is the Jacobi identity. 


The purpose of the following analysis is to show that the Lie algebra of the general linear group 
of dimension n over the real numbers (the group of n x n invertible matrices), GL (n, ®), can be 


'This chapter follows [5] 
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represented by a rotation group in an appropriate vector space, M (let (p,q) be the signature of 
the new vector space). Furthermore, since the rotation group, SO (p,q) can be represented by 
the spin group, Spin (p,q), in the same vector space. Then by use of geometric algebra we can 
construct any rotor, R € Spin(p,q), as R= e2 where B is a bivector in the geometric algebra 
of M and the bivectors form a Lie algebra under the commutator product.” 


The main trick in doing this is to construct the appropriate vector space, M, with subspaces 
isomorphic to #” so that a rotation in M is equivalent to a general linear transformation in a 
subspace of M isomorphic to R”. We might suspect that M cannot be a Euclidean space since 
a general linear transformation can cause the input vector to grow as well as shrink. 


9.2 Simple Examples 


9.2.1 SO(2) - Special Orthogonal Group of Order 2 


The group is represented by all 2 x 2 real matrices* R where* RR’ =I and det (R) =1. The 
group product is matrix multiplication and it is a group since if R,R’ = J and R,R} =I then 
(R, Ry) (Ri Ry)” = Rion Re = Ria = RyRy =f (9.4) 

det (R,R.) = det (R,) det (R,) = 1. (9.5) 

SO (2) is also a Lie group since all R € SO (2 
the coordinate 0: 


can be represented as a continuous function of 


) 
cos (#) —sin (6) 
(0 


sin(™) cos (@) 


In this case the spin representation of SO (2) is trivial, namely? R (6) = e29 = cos (5) + Isin (3), 
R(0)' = e-2® = cos (S) — Isin($) where J is the pseudo-scalar for G(R”). Then we have 
R(0)a=R(0)aR(0)'. 


?Since we could be dealing with vector spaces of arbitrary signature (p,q) we use the following nomenclature: 
SO (p,q) Special orthogonal group in vector space with signature (p, q) 
SO (n) Special orthogonal group in vector space with signature (n, 0) 
Spin (p,q) Spin group in vector space with signature (p, q) 
Spin (n) Spin group in vector space with signature (n, 0) 
GL (p,q,k) General linear group in real vector space with signature (p, q) 
GL(n,") General linear group in real vector space with signature (n, 0) 


R(@)= (9.6) 


3We denote linear transformations with an underbar. 
4TIn this case I is the identity matrix and not the pseudo scalar. 
°For multivectors such as the rotor R there is no underbar. We have R (a) = Ra = RaR! where a is a vector. 


9.3. PROPERTIES OF THE SPIN GROUP 149 


9.2.2 GL(2,8) - General Real Linear Group of Order 2 


The group is represented by all 2 x 2 real matrices A where det (A) 4 0. Again the the group 
product is matrix multiplication and it is a group because if A and B are 2 x 2 real matrices 
then det (AB) = det (A) det (B) and if det (A) 4 0 and der ) # 0 then det (AB) 4 0. Any 
member of the group is represented by the matrix (a = (a', a”, a°,a*) € R*): 


i 
A@=|% % |. (9.7) 
Thus any element A € GL (2,8) is a continuous function of a = (a',a?,a°,a*). Thus GL (2, ®) 
is a four dimensional Lie group while SO (2) is a one dimensional Lie group. 


Another difference is that S'O (2) is compact while GL (2, ®) is not. SO (2) is compact since for 
any convergent sequence (4;), 

lim R(6;) € SO (2). 

Ioo 
GL (2,8) is not compact since there is at least one convergent sequence (a;) such that 

lim det (A (a;)) = 0. 

1 0Co 


After we develop the required theory we will calculate the spin representation of GL (2, ®). 


9.3. Properties of the Spin Group 


The spin group, Spin (p,q), of signature (p,q) is the group of even multivectors, R, such that 
RR' =1. Then if R € Spin (p,q) the rotation operator, R, corresponding to R is defined by 


Rr = RxRt Vx € RP, (9.8) 


9.3.1 Every Rotor is the Exponential of a Bivector 


Now we wish to show that any rotor, R, can be written as R = e? where B is a bivector. 


To begin we must establish some properties of the rotor R, = e?« where B, is a 2-blade®, which 
is to say that B, defines a plane a by the relation x \ B, = 0, x € R?4. Since B, is a blade we 


°Remember that an n-blade is the wedge (outer) product of n vectors and defines and n-dimensional subspace. 
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have B? € ®. We evaluate R, as follows: 


Be 0% 
By 
i= e08 (v —B?) + sin (v —B?) = (9.9) 
Bes 03 


R, = cosh (/B2) + sinh (/B2) Ba (9.10) 


Likewise, if we have an expression of the form: 
Rag=S,+ Ba, (9.11) 


where S, is a scalar and B, a bivector blade we can put it in the form of eq (9.9) or eq (9.10) 
if S? — B? = 1. This implies that R = e+, where B? = +1 is a normalized 2-blade and 


6, = Sin (/—B2) or 6, = Sinh™ (./B2) depending on the sign of B?. 


Now consider the compound rotor R = R,R, = e?« ce? where both B, and B, are 2-blades. If 
B, and By, commute (the planes defined by B, and B, do not intersect) then 


R= R,Ry = e®* > = eFat Bs — 9 eFa — RR,, (9.12) 


and the rotors commute. If B, and B, do not commute the planes defined by B, and B, must 
intersect on a line as shown in figure (9.1). 


Define the basis vectors for intersecting planes as shown in figure (9.1) by 
e, A unit vector in the plane B, and is normal to ey. 


e,, A unit vector along the intersection of a and b. 
e, A unit vector in the By and is normal to ey. 


where unit vector means that the square of the basis vector can be +1. The metric tensor for 
vector space spanned by these basis vectors is 


e 0 a 
g=( 0 e 0 |, (9.13) 
a 0 e 
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e=+1 
ev =+1 
e? =+1 
€,° eG = a 


Figure 9.1: Basis for Intersecting Planes 


where a = e, : e;’. With this notation we have 


R, = e8°2em = C, + Sx€a€m (9.14) 

Ry = eFoemes =C+ SphemEd; (9.15) 

where Cy, Sa, Cy, and 5S, are the functions of 6, or 4 appropriate for the sign of the square of 
the bivectors €,€,, and e,,e,. Now evaluate R,Rb: 
R=R,R=S+B ( ) 

S = ae? S484 + CrCy (9.17) 

B= S,Cheaem + €? SaSr€a€b + Cy S,emep (9.18) 

(9.19) 

(9.20) 


B= S75? owe: 5.550.026 0 8 Hee OS Sees Ss 
ab m b-m~a'’b av~m~b\a a~bYa’b 
S? — B® = C20? + ee2,C2S? + e262, S20? + €2€292S?. 


a~m~ a 


Evaluating S* — B? for all combinations of e?, ef, and e?, and using trigometric identities to 
simplify S? — B? gives: 


7Some of the nomenclature associated with orthogonality is confusing. If a = 0 we say the planes defined by 
B, and By are orthogonal. However, if B, and By define subspaces with only the zero vector in common we say 
the subspaces are orthogonal. This is equivalent to B, - By = 0. 
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e ef e2, S? — B? 

-l -1 -1 sin? (0, ) sin” (O)) + sin? (0,) cos? (6) + sin * (2) co co s° (6a ) + cos? (8,) cos? (8)) = 1 

-1 -1 1 sinh? (6,) sinh? (6,) — sinh? (0,) cosh? (6,) — sinh? (4, \ cos? (0,) + cosh* (0,) cosh* (6,) 

-l 1 -1 — sin? (,) sinh? (0) + sin? (0,) cosh? (@,) — cos? (9) sinh? (6,) + cos? (8) cosh? (0)) = 

-l 1 1 — sin? (6)) sinh? (0,) + sin? (0) cosh? (0,) — cos? (0p) sinh? (0,) + cos? (04) cosh? (0,) = 

1 -1 -l — sin? (6p) sinh? (6,) + sin? (4) cosh? (8) — cos? (04) sinh? (@,) + cos? (0p) cosh? (6,) = 1 
1 -1 1 —sin? (0, ) sinh’ (9)) + sin? (0,) cosh? (@,) — cos? (0) sinh? (@,) + cos? (@,) cosh? (6,) = 1 
1 1. -1 sinh? (6,) sinh? (8 ) — sinh? (6,) cosh? (0) — sinh’ (Op) cosh” (0,) + cosh? (0,) cosh? (6,) = 1 
Ey ’at: a sin? (8,) sin? (,) + sin? (8) cos? (9) + sin? (@,) cos? (0,) + cos? (8,) cos? (O)) = 1 


This proves that if the blades B, and By, do not commute we can construct a 2-blade Bay such 
that 


B 
B* <0, Bap = Sin“ /—B?) — 
ePa,b =_ ePa ePo — S +- B, Bap — ca e (9.21) 


2 ee Arey 2 
B?>0, Bay =Sinh (vB?) a 


Now consider the general rotor R that is composed of products of the form 
Rae ase (9.22) 
where the B;’s are all 2-blades. 


Equation (9.22) is reduced to R = e? where B is a bivector (not necessarily a blade) via the 
following operations on adjacent terms in the product. 


1. If B; Biz. is grade 4, B; and B;,, define orthogonal subspaces and e”' and e?'+! commute. 
If needed reverse the order of e?' and e?+1 to get one or the other term closer to a term it 
does not commute with. 


2. If B; Bix is a scalar (grade 0), B; and Bj, define the same subspace and commute. Sub- 


stitute e?:+Fi+1 for e?: eFi+1 and remove term e? +, 


3. If B;Bj.1 is grade 2, B; and B;,, define intersecting subspaces and do not commute. Use 
q (9.21) to substitue e?4'+1 for e?! e?'+1 and remove term e?'+!, 


Repeat operations 1, 2, and 3 until R is in the form R = e*...e%* where s < r and all of the 
Bi commute with each other. Then we have 


R= cPitet Bs (9.23) 


and the assertion that ever rotation can be expressed as the exponential of a bivector is proved. 
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9.3.2 Every Exponential of a Bivector is a Rotor 


Now we need to prove the converse, that the exponential of a general bivector is the spin repre- 


sentation of a rotation. 


Let B be a general bivector and note that B! = —B. Then let R = e? so that 


phe OBO! a Bl) ey 
(8)! = @) ~ 9B) -y CB) 


i=0 , i=0 , i=0 
Likewise since B commutes with itself we have 


ae ae (9.24) 


RR — 2 (e?)! pO Se gh — AO 
Now consider the following function of \ where ag is any vector, 
AB 


a(Ay= e 2 age? 


and develop the power series expansion for a (A) 


d 1 AB AB 

a = Ao (—Bao + apB)e2 
=e 2 (ao - B)e® 

da _ AB AB 

i 2 ((ag9: B)- Bye? 

Ma _ AB dB 

De 2 (((a9: B)- B)- B)e2 


r 


1. (9.25) 


(9.26) 


(9.27) 


(9.28) 


(9.29) 


a 
Thus every derivative —— is a vector since the dot product of a vector and a bivector is always 


a vector and the Taylor series expansion of eq (9.26) about A = 0 is 


dM 
a(A) = a9 + a0A-B+ = (ao B)- B+. (9.30) 
Thus a (A) is a vector for any \ (take \ = 1) and R=e~® and Ray = e~7 age? so that we have 


(Rao) - (Rao) = (e ? age? | : (c ? age? | 


(9.31) 
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B 
2 


Thus e~2 preserves the length of ag and thus must be a rotation operator. 


9.4 The Grassmann Algebra 


Let VY” be an n-dimensional real vector space with basis {w;} : 1 <i <n and A is the outer 
(wedge) product for the geometric algebra define on V”. As before the geometric object 


DING Nisce INU (9.32) 


is a k-blade in the Grassmann or the geometric algebra where the v,’s are k independent vectors 
in Y, and a linear combination of k-blades is called a k-vector. The space of all k-vectors is just 
all the k-grade multivectors in the geometric algebra G (V,,). Denote the space of all k-vectors in 
v" by AF = A*(V") with 


dim (A‘) = (;,) - a (9.33) 


Letting A} = V” and A®° = & the entire Grassmann algebra is a 2”-dimensional space.® 


ve are (9.34) 
k=0 
The Grassmann algebra of a vector space Y” is denoted by A (V"). 
9.5 The Dual Space to VY, 
The dual space (V")* to V" is defined as follows. Let {w;} be a basis for V” and define the basis, 


{w7} for (V")" by 
bij. (9.35) 


8Consider the geometric algebra G (V”) of V” with a null basis {w;}: w? =0,1<i<n. Then 


(w;+ wy)? = w; + ww; + ww; + wi 
0= Wi Wj + Wj Wi 
Wj Wj = —WjWi 
0 = 2u; - w;. 


If the basis set is null the metric tensor is null and w;w; = w; A w; even if 7 = 7 and the geometric algebra is a 
Grassmann algebra. 
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Again the dual space Y"* has its own Grassmann algebra A*(V”) given by 
MOP =ae= yo Ar (9.36) 
k=0 


A* (Y") can be represented as a geometric algebra by imposing the null metric condition (w; \ =) 
as in the case of A(V"). 


9.6 The Mother Algebra 


From the base vector space and it’s dual space one can construct a 2n dimensional vector space 
from the direct sum of the two vector spaces as defined in Appendix F 


Ran =" @ (Yr) (9.37) 


with basis { wi, w't. An orthogonal basis for R”” can be constructed as follows (using equa- 
tion 9.35): 


e;-e; = (w; + w;) - (w; + w;) 


— * e P x 
= Wi Wj + Wi W; TW; W; TW; W,; 


= WW; — W; Ww; —W;- w+; Ww; 

=O; (9.41) 
e;- €; = (w; + w7) - (w; — vj) 

= Wi Wi — Wi W; + W; Wi WwW; W; 

= 0. (9.42) 


Thus R”” can also be represented by the direct sum of an n-dimensional Euclidian vector space, 
&", and a n-dimensional anti-Euclidian vector space (square of all basis vectors is —1) E” 


rn — BM @ BM (9.43) 
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In this case E” and E” are orthogonal? since e;-é€; = 0. The geometric algebra of ”” is defined 
and denoted by 
Rn = G(R") (9.44) 


has dimension 2?” with k-vector subspaces ee = G* (R™") and is called the mother algebra.'® 


From the basis {e;, €;} we can construct (p + q)-blades 


Ej, gah, b=, 8), (9.45) 

where!" 
Ep = Ci, €iz Ci = “p,0 DS as oe Se. (9.46) 
Eig = Cj, Cj «+ Ej, = Koa 1 <4, < Jo Se fats 74 <n. (9.47) 


Each blade determines a projection of E,,, of #" into a (p+ q)-dimensional subspace jv? 


defined by (see section 5.1.2 and equation 1.42)!2 
be 1 S 
Ei, g@ = (@- Epq) Bog = 5 (@- (—1)?*" EB, gay) - (9.48) 
A vector, a, is in 4 if and only if 
GN Epa = 0= ab, a+ (—1)"! Ey ga: (9.49) 


°Every vector in E” is orthogonal to every vector in By 
10 As an example of the equivalence of £"@ £” and V" SY" consider the metric tensors of each representation 
of R?:2. The metric tensor of E? 6 E? is 


10 0 0 
0 1 0 0 
0 0 -1 O 
0 0 0 -1 


with eigenvalues [1,1,—1,—1]. Thus the signature of the metric is (2,2). The metric tensor of V? 6 V?* is 


oOro°e 
Fe oOo7oc oe 


1 0 
0 1 
00 |” 
0 0 


with eigenvalues [1/2,1/2,—1/2,—1/2]. Thus the signature of the metric is (2,2). As expected both representa- 
tions have the same signature. 

‘1Since the {e;, €;} form an orthogonal set and specifying that no factors are repeated we can use the geometric 
product in equations 9.46 and 9.47 instead of the outer (wedge) product. 

!2The underbar notation as in E,, , (a) allows one to distinguish linear operators from elements in the algebra 
such as Ep q- 
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For p+ q =n, the blade E,,, determines a split of R”” into orthogonal subspaces with comple- 
mentary signature!’, as expressed by 


RE" = NPT | YW, (9.50) 


For the case of g = 0 equation 9.48 can be written as 


1 
Ea = 5 (a+a*), (9.51) 
where a* is defined by 
C=C Bae: 4, (9.52) 


It follows immediately that e* = e; and (€;)* = —é;.'* The split of R”” given by equation 9.37 
cannot be constructed as the split in equation 9.43 since the vectors expanded in the { wi, w* 
basis cannot be normalized since they are null vectors. Instead consider a bivector!® K in R”” 


K=)S°K,, (9.53) 
i=1 
where the K; are distinct commuting blades, kK; x AK; = 0 (commutator product), normalized to 
Kk? =1. The bivector K defines the automorphism K : R”" > #™” 
a= Kaz=axkK=a-K. (9.54) 


This maps every vector a into the vector @ which is called the complement of a with respect to 
Le 


Each K; is a bivector blade that defines a two dimensional Minkowski (since K? = 1) subspace 
of R™”. Since the blades commute, K; x AK; = 0, they define disjoint subspaces of R”” and since 


3The signature of R?-4 is (p,q) as opposed to R?*4 with signature (q, p). 
M4This is obvious by 


= (=1)"*? Be,Eo! = (—1)" 1)" 6 BB = C1)" a; =e, 
+1 = Caner ae (é;) E,E = Sl eee (é;) a (é;). 


15Tn general the basis blades for bivectors in Rp,» do not commute since the dimension of the bivector subspace, 


Re, is 
2n\ _ (2n)! i ee 
(3) = age re 1) = 2n* — n. 


If one has more than n bivector blades the planes defined by at least two of the blades will intersect. 


158 CHAPTER 9. LIE GROUPS AS SPIN GROUPS 


there are n of them they span R”". Since K? = 1 there exists an orthonormal Minkowski basis 
e; and é€; such that 


ee 0, (9.55) 


Then if a € R”” and using equations 9.55, 9.56, and from appendix B equations B.2 and B.5 we 
have 


a: K; =a: (e;€;) = (a : é;) e;+ (a * e;) e;, (9.57) 
(a- K;)-K; = (a- (e;€;)) - (e:€:) 
= (a- e;) e; — (a- &;) E; 
2 Gueljes Ge) eau (9.58) 
(a- Ki) - Kj = ((a- &) (ei: €;) — (@- e:) (Gi - €;)) e; (9.59) 
+ ((a- ej) (E&- ej) — (a &) (ej e;)) Ej; =0, Vi FZ 
Then Z 7 
Ka=)/a-Ki= ((a- e;) €; — (a: &) e;) (9.60) 
and!® 


=a. (9.61) 


16Remember that since K; defines a Minkowski subspace we have for the reciprocal basis e’ = e; and é? = —é;. 
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Also 


3 


= 0, (9.62) 


a =" ((a-e;)’ - (a-&)’), (9.63) 


a= 5 ((a-&)° —(a-ei)’), (9.64) 
i=1 

a’ +a’ =0. (9.65) 

Thus defining 
a4, =atG@=atKa=axta:K, (9.66) 

we have 
(axz)” = 0, (9.67) 
a,-a.=0-@ =25 (a-e)’. (9.68) 
i=1 


Thus the sets {a;} and {a_} of all such vectors are in dual n-dimensional vector spaces (a; € V” 
and a_ € Y¥"*), so K determines the desired null space decomposition of the form in equation 9.37 
without referring to a vector basis.!” 


The K for a given #”” is constructed from the basis given in equations 9.38 and 9.39. Then we 
have from equation 9.60 


‘The properties K; x K; =0 and K? = 1 allows us to construct K from the basis {e;, é€;}, but do not specify 
any particular basis. 
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Also from equation 9.60 


Kw; = (w;- e;) €; — (w;- €;) e;, 
= 5 (((e:+6)-€:)&—((e +€) Ae), 
lene 
= 5 (€; + e;), 
= Wi, (9.71) 
Kw; = (w; €;) e; — (w; é;) ej, 
1 
= 3 “ = 6;) + €))€;— ey =e, ) +6;) 2:4 
er (e; _ é;) ; 
= —w7}. (9.72) 


The basis {w;, w;} is called a Witt basis in the theory of quadratic forms. 


9.7 The General Linear Group as a Spin Group 


We will now use the results in section 1.15 and the extension of a linear vector function to 
blades (equation 1.77)!® and the geometric algebra definition of the determinant of a linear 
transformation (equation 1.79) 1% 


We are concerned here with linear transformations on #”” and its subspaces, especially orthog- 
onal transformations. An orthogonal transformation R is defined by the property 


(Ra) - (Rb) =a. (9.73) 
R is called a rotation if det (R) = 1, that is, if 
HPs as Lesa: (9.74) 


where En.» = EB is the pseudoscalar for R,, (equation 9.45). These rotations form a group 
called the special orthogonal group SO(n). 


From section 1.10 we know that every rotation can be expressed by the canonical form 


Ra = RaR’, (9.75) 


ee ) 


=f@AL@A.. 
det (f) Iw 


here I is the pseudesoalie: 
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where FR is an even multivector (rotor) satisfying 
RR =1. (9.76) 


The rotors form a multiplicative group called the spin group or spin representation of SO(n), and 
it is denoted by Spin(n). Spin(n) is said to be a double covering of SO(n), since equation 9.75 
shows that both +R correspond to the same R. 


From equation 9.76 it follows that R~' = R' and that the inverse of the rotation is 
Ria = RraR. (9.77) 
This implies that from the definition of the adjoint (using RR5 in appendix C) 
a- (Rb) = (aRbR') = (bR'aR) = b- (Ria). (9.78) 
The adjoint of a rotation is equal to its inverse. 
“It can be shown that every rotor can be expressed in exponential form 
R=+e2, with Rt =+e72?, (9.79) 


where B is a bivector (section 1.10.2) called the generator of R or R, and the minus 
sign can usually be eliminated by a change in the definition of B. Thus every bivector 
determines a unique rotation. The bivector generators of a spin or rotation group 
form a Lie algebra under the commutator product. This reduces the description 
of Lie groups to Lie algebras. The Lie algebra of SO(n) and Spin(n) is designated 
by so(n). It consists of the entire bivector space R7,,,. Our task will be to prove that 
and develop a systematic way to find them. 

Lie groups are classified according to their invariants. For the classical groups the 
invariants are nondegenerate bilinear (quadratic) forms. Geometric algebra supplies 
us with a simpler alternative of invariants, namely, the multivectors which determine 
the bilinear forms. As emphasized in reference [4], every bilinear form can be written 
as a: (Q ()) where Q is a linear operator, and the form is nondegenerate if Q is 
nonsingular (i.e., det (Q) 4 0).” 7° 


To prove that e2 is a rotation if B isa general bivector consider the function?! 


a(A)=e2 ae 2 (9.80) 


20From reference [5]. 
21 We let a = a (0). 


162 CHAPTER 9. LIE GROUPS AS SPIN GROUPS 


where a is a vector. Then differentiate a (,) 


da B AB AB AB 1s B 
= 9g 


ay. 3 
= 5a (A) —a(A) 5 
= B-a(d) (9.81) 
ne = 8: (GG) = 8 (Ba) (9.82) 
d'a 
ca [eas ee cee al cae a (9.83) 


r copies of B 
so that the Taylor expansion of a (1) is 


1 
-? =a+B-a+—B-(B-a)+ 


ae oy 


gb (B (Bea) +e. (9.84) 
Since every term in eq (9.84) is a vector then e2 ae? 
then (using RR5 in appendix C) 


B B B B B B B B 
ezae 2)-(e2be 2)=(fe2zae 2e2 be 2 


is a vector. Finally let a and b be vectors 


ep eh (9.85) 


B 
2 


Ba : Bo SRS . B 
Thus e2 is a rotor since e2 ae 2 is a vector and the transformation generated by e2 preserves 


the inner product of two vectors. 


For an n-dimensional vector space the number of linearly independent bivectors is (5) = Tom = 


, which is also the number of generators for the spin group Spin (n) which are ea = 
cos (5) + sin 4) e;\e; V 1<i<j<nwhere 6;; is the rotation in the plane defined by e; and 
e,;. Thus there is a one to one correspondence between the bivectors, B, and the rotations in an 
n-dimensional vector space. 


n2—n 


Q is invariant under a rotation R if 


(Ra) - (QRb) = a- (Qb). (9.86) 
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Using equation 9.78 transforms equation 9.86 to 


or since the a and b in equation 9.87 are arbitrary 


RIQR = Q = RQR', 


RR'QR = RQ, 
QR = RY. 
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(9.87) 


(9.88) 
(9.89) 
(9.90) 


Thus the invariance group of consists of those rotations which commute with Q. 


This is obviously a group since if 


RQ a QR, 
RQ =; QRo, 


then 


R, RQ = RiQRo, 
RRQ re QR, Ro. 


e Qb = b*. Then 


Eon. 
T B ROR ES", 


T READER, 


As a specific case consider the quadratic form wher 
Qb = (-1)" 
QR = (-1)" 
RQ = (-1)" 
E, RoR! E,' = RE,bE,'R', 
Eh RE ys 
fy = REE. 


r 


(9.95) 


(9.96) 


Thus, the invariance of the bilinear form a- b* is equivalent to the invariance of the n-blade E,,. 
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We can determine a representation for R by noting the generators of any rotation in #”” can be 
written in the form 


R=e™, 
= { conto J * sabi J 


e sooh (6/3) bad Gn, pS (u'e; + w'e;) (v’e; + 0'e;), (9.97) 


i<j 


where u-v = 0 and | (wv)? | = 1. Equation (9.97) is what makes corresponding the Lie algebras 
with the bivector commuatator algebra possible and also allows one to calculate the generators 
of the corresponding Lie group. The generators of the most general rotation in R”” are e;e,, 
é;e;, and €;é;. The question is which of these generators commute with E,,. The answer is 


eje;E, = Eneie; (9.98) 
exe; i, = —,e,€; (9.99) 


Equation 9.100 is obvious since e; and €; give the same number of sign flips, n, in passing through 
E,, totalling 2n an even number. Likewise, in equation 9.98 as e; and e; give the same number 
of sign flips, n — 1, in passing through E,, totalling 2(n — 1) an even number. In equation 9.99 
e; produces n — | sign flips in traversing E,, and e; produces n sign flips for the same traverse 
so that the total number of sign flips is 2n — 1 and odd number. 


Thus the reqired rotation generators for R are 
ei; = eje;, for l< j = 1, aotany TU, (9.101) 
e;; = E;€;, fori < j = 1, se nabs (9.102) 


Also note that since e;; x €4; = 0, the barred and unbarred generators commute. Any generator 
in the algebra can be written in the form 


B=a:e+B:6é, (9.103) 
where - 
a:e= Ss" aes, (9.104) 
i<j 


and the a” are scalar coefficients. So that @ : e is the generator for any rotation in R” and 
(3: € is the generator for any rotation in R”. The corresponding group rotor is (we can factor 


9.7. THE GENERAL LINEAR GROUP AS A SPIN GROUP 165 


the exponent since e and € generators commute) 
R= erler Be) — ppaie Qi He (9.105) 
This is the spin representation of the product group SO(n) @ SO(n). 


In most cases the generators of the invariance group are not as obvious and in the case of the a-b* 
form. Thus we need some general methods for such determinations. Consider a skew-symmetric 
bilinear form GO 


a- (Qb) = —b- (Qa). (9.106) 
The form Q can be written 
a+ (Qb) =a-(b-Q) = (a@Ab)-Q, (9.107) 


where Q is a bivector.?? We say that the bivector Q in involutory if Q is nonsingular and 


eae al (9.108) 


Note that the operator equation 9.108 only applies to vectors. 


Note that 


(Ra \ Rb)-Q=R(aNb)-Q 


= (Rt R(a/Nb) ) R'). Q 
= (R(a/Nb) R'Q) 
= ((aA b) R'QR) 
=(aAb)- (RQ). (9.109) 


22This selection is done with malice aforethought. To generate the memembers of GL (n, R) we do not need the 
most general linear transformation on R” since the dimension of that group is 4n? and not the n? of GL (n, R). 

?3This is obvious from the properties of the bilinear form that if a-(Qb) = (aA b)-Q then a- (Qb) is linear in 
a and b and is skew-symmetric 


b- (Qa) = (bAa)-Q 
—(aAb)-Q 
=a: (Qi 


Finally the maximum number of free parameters (coefficients) for the bivector, Q, in an n-dimensional space 


is ee) a mint) This is also the number of independent coefficients in the n x n antisymmetric matrix that 


represents the skew-symmetric bilinear form. 


166 CHAPTER 9. LIE GROUPS AS SPIN GROUPS 


Then equation 9.107 gives for a stability condition 


(Ra) - (QRb) = a: (Qb) , 
((Ra) - (Rb))-Q = (ab) -Q, 
((Ra) A (Rb) -Q = 


(aA b)- (R'QR) =(aAb)-Q, 


ROR=@. 
QR= RQ (9.110) 


From equation 9.110 we have that generators of the stability group G(Q) for Q must commute 
with Q. To learn more about this requirement, we study the commutator of Q with an arbitrary 
bivector blade a A b. Since aA b=a x band ax Q=a-Q the Jacobi identity gives 


(a\b)xQ=(axb)xQ, 

x Q) x b+ax (bx Q), 

x Q)Ab+aA(bxQ), 

=(a-Q)Ab+aA(b-Q), 

= (Qa) Ab+aA (Qb), (9.111) 


=a 
=(a 


and then 
((a Ab) x Q) x Q = ((Qa) Ab+ aA (Qd)) x Q. (9.112) 


Now using the Jacobi identity and the extension of linear functions to blades we have 


((Qa) Ab) x Q = ((Qa) x b) x Q, 

= (Qa) x (bx Q) —b x ((Qa) x Q), 
= (@a) A (b-Q) —bA ((Qa) -Q), 
= (Qa) A 
=(@ 


(Qb) —bA (Q?a), 
a) b+ Q(aAb). (9.113) 


a 


Similarly 
(a A (Qb)) x Q=Q (ab) +aA (Q’b) , (9.114) 
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so that, since Q is involutary (Q? — +1) 


((a Ab) x Q) x Q = (Q?a) Ab+2Q (ab) +aA (Q’b), 
=2(Q(aAb)+aQNb). (9.115) 


By linearity and superposition since equation 9.115 holds for any blade aA 6 it also holds for any 
bivector (superposition of 2-blades) B so that 


(Bx Q)xQ=2(QB+B). (9.116) 
If B commutes with Q then B x Q = 0 and 


QB+B=0, (9.117) 
QB==B. (9.118) 


Thus the generators of G'(Q) are the eigenbivectors of Q with eigenvalues +1. 


Now define the bivectors E* (a,b) and F (a,b) by 


E* (a,b) =aAb+ (Qa) (Qd), (9.119) 
F (a,b) = (Qa) Ab—aA (Qb). (9.120) 


Then using equation 9.111 we get 


B* (a,b) x Q = (Qa) Ab +.aA (Qb) ¥ (Qa) A (Qb) = (Qa) A (QB), 
= (Qa) Ab+an (Qb) —aN (Qb) — (Qa) \ 0, 
= 0, (9.121) 
F (a,b) x Q = (Qa) A (Qb) + aA (Q°b) — (Qa) Ab— (Qa) A (QS), 
Ska Abe. Vb; 
=0. (9.122) 


Thus E* (a,b) and F'(a,b) are the generators of the stability qroup for Q. A basis for the Lie 
algebra is obtained by inserting basis vectors for a and b. The commutation relations for the 
generators E* (a,b) and F' (a,b) can be found from equations 9.119, 9.120, and B.13. Evaluation 
of the commutation relations is simplified by using the eigenvectors of Q for a basis, so it is best 
to defer the task until Q is completely specified. = 
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As a concrete example let Q = K where K is from equation 9.53 and use equations 9.69 and 
9.70 to get?* 


Ei; = Et (ex, €;) =e; Ae; — (Kke;) A (Ke), 


ep = F (ei, e;) =e; /\ (Kve;) = (Ke;) A ej, 
1 


At this point we should note that if a bivector B is a linear combination of Ej, Fi;, and F;, B 
will commute with Q. Then R = e? also commutes with Q since ez only can contain powers of 
B. Also Ej;, Fi;, and F; are called the generators of the Lie algebra associated with the bilinear 


form defined by K and the number of generators are n= ne nen tn=n?. 


The structure equations for the Lie algebra (non-zero commutators of the Lie algebra generators) 
of the bilinear form K are 


Ei X Fi =2(Ki — K;) (9.126) 
Ei; x K, = —Fy (9.127) 
F, x K; = -E; (9.128) 
Ei; Xx Ey =—Ej (9.129) 
Fy x Fa = Ej (9.130) 
Fy; x Eq = Fi (9.131) 


The structure equations close the algebra with respect to the commutator product (see ap- 
pendix G for how to calculate the structure equations). 
Thus (using the notation of eq (9.104)?°) the rotors for the stablity group of K can be written 


er = er (CEB Fty:K) (9.132) 


where the n? coefficients are aj;;, 3ij;, and 7;. 


The stability group of Kk can be identified with the general linear group GL(n, #). First we must 
show that A does not mix the subspaces V” and Y™ of R””. Using equations 9.71 and 9.72 we 


4Since K? = 1 we only need consider E+ (e;,€;). 
2°The a: E notation simply says that indices of the scalar coefficients, a;;, are balanced by the indices of the 
bivectors, £;;. So that a: EB = en OH Ory eK = Sy: 
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have (where W,, and W; are the pseudoscalars for V" and V"*). 


W,, = W 1... Wn, (9.133) 
W, = wy... Ww, (9.134) 
K (Wy) = (-1)" Wt. (9.136) 
Thus when restricted to V” or VY", K is non-singular since the det (’) 4 0. 
Since each group element leaves Y” invariant 
RK (W,,) = RW, = RW,,R' =W,, (9.137) 
we can write a 
Rw; = S— wrpry- (9.138) 
k=1 
Where using eq (9.35) we get 
pij = 2w} - (Rw;) =2(w7Rw;R'). (9.139) 


The rotations can be described on V” without reference to #”. For any member FR of GL (n, R) 
we have 


where detyn is the determinant of the linear transformation restricted to the V” subspace of #”” 
and not on the entire vector space. First note that 


Ws -Wi = (wi... wt) - (wp... w1) = (wt... wkw,...wy) = 2. (9.141) 
Therefore, 


RW; = Widetyn (R™) (9.142) 


9.8 Endomorphisms of ft” 


For a vector space an endomorphism is a linear mapping of the vector space onto itself. If the 
endomorphism has an inverse it is an automorphism. By studying the endomorphisms of R,,, the 
geometric algebra of #R” we can also show in an alternative way that the mother algebra, Rn, 
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is the appropriate arena for the study of linear transformations and Lie groups. First note that 
(~ is the symbol for “is isomorphic to” ) 

End (Rn) ~ RK”, (9.143) 
since End (R,,) is isomorphic to the algebra of all 2” x 2” matrices. 


For an arbitrary multivector A € ®,, left and right multiplications by orthonormal basis vectors 
e; determine endomorphisms of ,, defined by 


e,: A> e, (A) =e,A, (9.144) 
é,: A é, (A) = Ae; (9.145) 

and the involution operator (A = A) is defined by® 
(AB) = AB (9.146) 


e,e; (A) = e; (e;A) = exe; A, (9.148) 
e,€; (A) =e, (Ae;) = e;Ae;, (9.149) 
€,e; (A) = €; (eA) = e,Ae; = —e;Ae;, (9.150) 
é,€, (A) = &; (Ae;) = Aeje; = —Aeje; (9.151) 
Thus 
(e,e; ats €,€;) (A) = (e,e; + €;€;) A= 20;;A, 
ee; + ee; = 20;4;, (9.152) 
(e,é; é,€;) (A) = e,;Ae; — e,Ae; = 0, 
C6) or Ee = 0, (9.153) 
(€,€; + €;€;) (A) = —A (eye; + exe;) A = —26,A, 
6Since any grade r basis blade is of the form e;, ...e;, with the j,’s in normal order then 
@;,..-€j, = €;, ...€;, = (-1)" ej, ...e;, 
and the involute of any multivector, A, can be determined from equations (9.146) and (9.147). Likewise 


Cj, ++ €j, = €j,.--€j,- 
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Thus equations (9.152), (9.153), and (9.154) are isomorphic to equations (9.40), (9.42), and (9.41) 
which define the orthogonal basis {e;, é;} of R””. This establishes the isomorphism?” 


Ran ~ End (Rn). (9.155) 


The differences between F,,,, and End (R,,) are that Ry» is the geometric algebra of a vector space 
with signature (n,n), basis {e;, €;}, and dimension 27”, while End(,,) are the linear mappings 
of the geometric algebra #,, on to itself with linear basis functions {é;, é,} that are isomorphic 
to {e;,é;} and have the same anti-commutation relations (dot products) as {e;, €;}. 


Note that in defining é,(A) = Ae;, the involution, A, is required to get the proper anti- 
commutation (dot product) relations that insure e; and €; are orthogonal and that the signature 
of End(R,,) is (n,n). 


Additionally, the composite operators 


e,€;: A > €,€; (A) = eAe; (9.156) 
e,€; (AB) = e,ABe; 
= e;Ae,e;Be; 
= €,€; (A) €,é; (B) (9.157) 


preserve the geometric product and generate Aut (R,,), a subgroup of End (R,). 


edit | 2% 


Chapter 10 


Classical Electromagnetic Theory 


To formulate classical electromagnetic theory in terms of geometric algebra we will use the space- 
time algebra with orthogonal basis vectors Yo, 71, Y2, and 3 and indexing convention that latin 
indices take on values 1 through 3 and greek indices values 0 through 3. The signature of the 
vector space is given by y6 = —y? = 1 and the reciprocal basis is given by y° = yo and 7‘ = —7;. 
The relative basis vectors are given by ¢; = 7;7o. These bivectors are called relative vectors 
because their multiplication table is the same as for the geometric algebra of Euclidean 3-space 
(section 1.8.1) since 6? = 1 and G4; = —Gj4;. Vector accents are used for relative vectors 
since it is usually used to denote vectors in ordinary (Euclidean) 3-space. We also have that 


F10203 = Yo 273 = I. 


We now define relative electric and magnetic field vectors 


EB = Ey = E*yi = E'G; (10.1) 
We also have 
V0 = (7°00 +. ¥'8:) Yo (10.3) 
=-V (10.5) 
where we have defined Z 
and also have 7 
YoV = Oo + YWY'O: = A — WV: = Oo + ViVi = 00 + V (10.7) 
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Finally we have the 4-current J = Jy, where J° = p then 
Jy = pt Jv = pt Jia; = pt J. (10.8) 


10.1 Maxwell Equations 


The four vacuum Maxwell equations are then given by (for this section x is the 3-D vector 


product and remember in three dimensions a x b = —I(a A b)) 
V-E=p (10.9) 
V-B=0 (10.10) 
-V x E=IVAE=&B (10.11) 
Vx B=-IVAB=HE+S (10.12) 
Now multiply equations 10.10 and 10.11 by J and —I respectively and add 10.11 to 10.12 and 


10.9 to 10.10 to get 


V+ (B+ 1B) =p (10.13) 
Va (£418) =-a (2 +18) -F (10.14) 
VA (E+1B) +0 (B+ 1B) =-F (10.15) 

Let 7 7 
f= hah (10.16) 

so that 

V-F=p (10.17) 
VAF+0F =—J. (10.18) 


Now remember (7 = ¥°) 
VF=VAF+V-F 
VF =-F +p—J 
VF+F =p-J 
wVF =p-J 
VF =p%0-— J'y0¥iV0 
VF =pyt Jini (10.19) 
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or 


VF=J. (10.20) 


Equation 10.20 contains all four vector Maxwell equations. Evaluating F' gives 
F = E'yoy + E072 + E*y0713 + BY yas — Beis + Beye (10.21) 
with the conserved quatities E? — B? and E'- B since 
FP? = B?_ B49 (é B) faarR (10.22) 


and both scalars and pseudoscalars are unchange by spacetime rotations (Lorentz transforma- 
tions). 


The components of F’ as a tensor are given by F¥” = (y” Ay") - F 


0 —-E! —E? —E 
E' 0 -Be B 
ee 0 —B! 
He =p? Bp 0 


Fu = (10.23) 


10.2. Relativity and Particles 


Let x (A) = x” (A) y, be a parameterization of a particles world line (4-vector). The the proper 
time, T, of the world line is the parameterization A (7) such that 


(2) =A. (10.24) 


Since we are only considering physically accessible world lines where a particle velocity cannot 
equal or exceed the speed of light we have 


(2) > 0. (10.25) 
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The function \ (7) can always be constructed as follows - 
dx _ dxdx 


dx\? dx\? (dd\? 
(=) =(&) (B) = oe 
dX 1 
= ; (10.28) 
( 
Me 
(=) 
ae 10.29 
[ (=) ( ) 


From the form of equation 10.29 and equation 10.25 we know that 7 (A) is a monotonic function 
(7 (Az) > T (Ai) V Ag > Az) so that the inverse function (7) exists x (7) = x (A(7)). 


Now define the 4-velocity v of a particle with world line x (rT) as 


dx dx” 


= FS y=, 10.30 
"dr dr ee? aes 


the relative 3-velocity (bivector) B as (remember that x° is the time coordinate in the local 
coordinate frame and 7 is the time coordinate in the rest frame of the particle defined by 8 = 0) 


- _ dx’ a _ dz _ 
and the relativistic y factor as 
dx® 
ee (10.32) 
dt 
Thus we can write 
dx® dx* 
= j 10.33 
Cee (vo a) (10.33) 
u=Sy (v0 +. Bro) a (1 a 3) Yo (10.34) 
v= (1-8) =1 (10.35) 
1 
Dr z 10.36 
ate: (10.36) 
1 
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We also define the relativistic acceleration 0 by 


_ adv 
v= —. 
dt 
Note that since v? = 1 we have 
d d 
=, U0) =v = Ww-H=0. 


10.3. Lorentz Force Law 


From section 10.2 we have for the relativistic 4-velocity 


= =7(1+8) % 


Then the covariant Lorentz force law is given by (q is the charge of the particle) 


dp _ dx 


gq Faq (EA + E'B + EB") 90 


dr ‘dr 


= (-B’e° ae B88? 1 E") a 
i (B'B° —B*6} + E?) Yo 
L (-B’? nS oy + FE) 73) 


= (EA) w+ ((B+5%8)-8)) 


ERC 


(10.38) 


(10.39) 


(10.40) 


(10.41) 
(10.42) 


where the component gy (E- 8) is the time derivative of the work done on the particle by the 


electric field in the specified coordinate frame. 


10.4 Relativistic Field Transformations 


A Lorentz transformation is described by a rotor L = BR which is composed of a pure spatial 


rotation R and a relativistic boost rotation B. The general form of B is 


B=cosh (=) — sinh (5) B 


(10.43) 
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where ( is the magnitude of the relative 3-velocity (bivector) and 8 is a unit relative vector 
(bivector) in the direction of the 3-velocity. Then 


6 = tanh (a) (10.44) 

1 
cosh (a) = y = ——== 45 
sinh (a) = ¥8. (10.46) 


If a boost is performed the relationship between the old basis vectors, y,, and the new basis 
vectors, ¥), is 


\, = By B. (10.47) 


If a boost is performed on the basis vectors we must have for the electomagnetic field bivector 
(remember BBt = B'B = 1) 


F= Fey A, = Fy ny, (10.48) 
F = F*” By, Bt By, Bt = FY’ By, Bt = FY yyy (10.49) 
BYFB = Fy, = FY Bly,yB. (10.50) 


As an example consider the case of B= Yio and 


B= cosh (5) — sinh (=) 710 (10.51) 


where the velocity boost is along the 1-axis (x-axis). Then (after applying hyperbolic trig iden- 
tities and double angle formulas) 


BIFB= E'yon 

+ (—B? sinh (a) + E” cosh (a) 
+ (B* cosh (a) — E* sinh (a)) 
+ (B’ sinh (a) + E* cosh (a)) yoy 

+ (—B? cosh (a) — E* sinh (a)) 1173 

+ Byo73. (10.52) 


YoY2 
172 


an 


2S. 3 


nS” 
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Now using equations 10.44, 10.45, and 10.46 to get 
BIFB = BX 

+7 (8B + B°) vor 

+ y (—BE? + B°) 7% 

+ (6B? mig E°) YoY3 

ian (-BE° = B?) W173 


+ Biyays (10.53) 
Equating components of B'F'B gives 
Ei = Fl Bi pl 
EB? = 7(E?— 6B’) B? =7(B? + BE) (10.54) 


ES =7(E5+ 8B?) B’ =+7(B — BE?) 


for the transformed electomagnetic field components. 


10.5 The Vector Potential 


Starting with 


VFHaV-F+VAF=J (10:55) 
and noting that V- F' and J are vectors and V A F is a trivector we have 
Weed, (10.56) 
VAF= 0. (10.57) 
By equation 3.13 we can then write 
F=VAA (10.58) 
VRESVNY KA = 0. (10.59) 


The equation V - F' = J gives (use the fact that VG = VAG+V-G where G is a multivector 
field) 


J=V-(VAA) 


Aa OY AY, (10.60) 
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We also note that if A is a scalar field then 
F=VA(A+VA)=VAA (10.61) 


since V A (VA) = 0. We have the freedom to add VA to A without changing F’. This is gauge 
invariance. 


The Lorentz gauge is choosing V - A = 0, which is equivalent to V - (VA) = —V- A, so that the 
equation for the vector potential becomes 


VWA= J. (10.62) 


pdf 


10.6 Radiation from a Charged Particle 


We now will calculate the electromagnetic fields radiated by a moving point charge. Let 2, (7) 
be the 4-vector (space-time) trajectory of the point charge as a function of the charge’s proper 
time, T, and let x be the 4-vector observation point at which the fields are to be calculated. The 
4-vector separation between the observer and the radiator is 


X =x -%p(T) = (a? - es (7)) Yo + (F = 7%) Y0- (10.63) 


The critical relationship involving X is that it must lie on the light cone of x, (7) since we are 
dealing with electromagnetic radiation. That is X must be a null vector, X? = 0. X (x,r)* =0 
implies there is a functional relationship between 7 and x so that we may write 7, (x). 7, (x) is 
the retarded time function. We use retarded time to mean the smaller solution (earlier time) 
of X (2, ty? = 0 since the equation must always have two solutions. The larger solution would 
be the advanced time and would correspond to signal travelling faster than light. In all these 
calculations c = 1 and the relative 3-velocity (bivector) of the moving charge is By (section 10.2). 


The 4-velocity and acceleration of the point charge are given by 


dx = 
i= ae = Fy, (1 + B,) vo (10.64) 


du 
US aE 


(10.65) 
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The vector potential of a moving point charge of charge q is given by the Liénard-Wiechart 
potential! 


1 
yp (r,t) = 7 = (10.66) 
™ (1-a- 8) |F- 7 r 
A(7,t) = 8 (% (2) 9 (Ft) (10.67) 
where S 
me C= Tp 
i= 2 (10.68) 
7 = i 
or in terms of 4-vectors 
4 rae q Vv 
A=( ie A(t) = 10. 
y (7,t) + A(F,t)} Yo rer (10.69) 


and we must calculate V A A. In the following derivation whenever we write X we mean 
X (x,7,(x)). X is always evaluated at the retarded time. 


To start consider that if we know 7, (x) implicitly defined by X (x,7,)” = 0 then the function 
X (x, 7; (x))* = 0 (identically equal to zero). Thus V (X (a7 (x))”) = 0. First note that (by 
the symmetry of gy) 
V (Xx?) = "Or, (AEX gi) 
4 ((OyX#) XY + XOX”) Gyn 

SOL Xe Gig 

= 24" (O,X) -X (10.70) 
Thus y" (0,X)-X = 0 is equivalent to V (X?) = 0. But 


(OX) AX Sy" (Op) * XS" (Open) eX 


-1 (Sn) un (Ss) 


p OTr 
= "00 X" Guy — Y" ( i vm) x 


OT, Ox" ut 
Or, {Ox 
= gw" X" — 9" ann Ez .) 7.4 
= X — V7, (up: X) =0 (10.71) 


‘en.wikipedia.org /wiki/Liénard-Wiechert_potential 
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or 


VS : 10.72 
a orn (10.72) 
Another simple relation we need is 
Vy, = nue 
Up —< ey. Ort Vv 
7 yOUp OTp 
OT, Oxt '” 
= (V7) Up 
Xb, 
= : 10.73 
ra (10.73) 
The last quantity we need is V (X - v) 
V(X -v) = "0, ((2” — a) A" Gin) 
On" 08,07 Ov?! Or, 
= vb EE 7 > Gee eas 
: (= OT, a] ena OT, Oxt! ae ,) 
. Or, Oxe OTe sty AOU 
= FInd — 7" Ox ar, Sunitp 1 Voce HH "Or, 
= Up — VT (Up * Up — X + Up) 
a aa a a a 
X (X + %, — 1) 
= 10.74 
Up xX P Up ( ) 


Now we can calculate VA 
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An (X + Up 
—— v _ XA Up — (X + Up) X AU 
Am (X + Up)” (xa ce X + Up ) 
Pe. ) Uv -u,) — av ‘i 
= Fe ay A Ate + (K Ntp) (Xty) — (Kp) (KAp)) (10.75) 


Note that equation 10.75 is a pure bivector so that V-A =0 and VA A=VA so that F = VA. 
Now expand (X A Up) (X - Up) — (X - Up) (X A vp) using the definitons of A and - for vectors to 


get 
q 1 . 

Fo= — , (X Av,+ =X (v1, Av,)X ]. 10.76 
ee p+ 5X (yA op) X)] (10.76) 


1 
The first term in equation 10.76 falls of as — and is the static field term. The second term falls 
r 


1 
off as — and is the radiation field term so that 
r 


X (i x 
pe Ag (10.77) 
8 (X - Up) 


Appendix A 


Further Properties of the Geometric 


Product 


For two vectors a and b we can revise the order of the geometric product via the formula 


ba = 2a-b— ab. 


By repeated applications of equation A.1 we can relate the geometric product aa; 
geometric product a, ...a,a - 


AQ, ...A,p = 2(A-+ a1) AQ... A, — G1 AAg... Ay 


= 2(a-a1)a2...a, — 2(a- ag) aa3... a, + 41A2043... A, 


r 


250 (-1)"" (@- ag) a... Op + (—1)" ay. . apa. 


k=1 
Thus 
= (aa Gg 201) ee yen re Gi cate 
: k=1 
- k+1 
a-(@1...4-) = (—1)""" (a: ag) 1... Gg... Oy 
k=1 


Now consider the expression a- (a; A... / a,) and how it can be reduced. Since a; A 
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(A.1) 


...@, to the 


(A.2) 


(A.3) 


... Aa, only 
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contains grades r, r — 2, ..., we have 


= 5 (aay. -dy ~ (1)! ay. ga) 
9 Gre pW? (Oy ALO phn g Os (A.4) 


The term we need is the r — 1 grade part so that 


1 r 
a-(ai1A...AQ@,) = = (aa,...a, —(-1) a1... aya), 4 


2 
= aay cere een fe 
= So (-1)" (a-an) (iA... Nag A... A Gr) 
= (1? (aay) (avn Aden 2 Nae). (A.5) 


If A,, B,, and C; are pure grade multivectors of grade r, s, and t we have 


A,-(Bs-C:) =(A-ABs)-G:;Vr+s<tandr,s>0 (A.6) 
A,-(B,-C,) =(A,- Bs)-C;Vrt+t<s. (A.7) 


To prove equation A.6 


Ae (B, . aS Ay: ns ie 


gos 


(A, B.C) ,_ ee higher grades) 


S {A BiOe) secs (A.8) 


and 


(A, \ B;) S Ch= a. . Ct 


Ee Ct) .-( (r+s) 


(Ar 
((A 
(is (ApBsUie. Gia higher grades) 


t-(r-+s) 
=(ABsCt)s (rts) (A.9) 


187 


To prove equation A.7 


A, + (B,- Cy) = Ap: (BsC;) 
= (A, (B sCt) Nes a (r+t) 
= A,BsCt) . (p44) + higher grades) 
= (A, BsCz) ._(r42) (A.10) 
and 
(A, : Bs) -C, : = (A, a : ri Ct 
=((A Ct), —(r+t) 
= (is A,BsCt) . (p44) + higher grades) 
= (A, B Ct). (r+t) * (A.11) 


If A, is the pseudoscalar of an n-dimensional subspace of the vector space and B; is a blade 
(Hestenes puts a line over the subscript or superscript to indicate that the pure grade multivector 
is a blade or as he says is simple) then 


First remember that Ps, (B) = (B-A,)A;!. Since both A, and B; are pede they both 
can be written as the geometric product of orthogonal vectors. Let A, = u;'...u4 and Bs = 
uf...uAuP_,...u?. The expression for B; indicate that it could contain r of the same basis 


vectors of A,. Also remember that 


A A 
7 os u 
= = 7 (A.13) 
(un) (uf) 
and evaluate 
BAA tite ge Uy te (A.14) 


Since r <n the number of orthogonal vectors in equation A.14 and hence the grade of B;A,, is 
n+s—2r>0. If s <n and subspace defined by B; is contained by the subspace defined by A, 
then the grade of 


BA He «tees on US (A.15) 
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isn —s and B;A, = B;-A,. Thus 


a ee. 


Now prove that for the integers 1 <7, <---<i,<nand1 <j, <--- <j, <n we have 


CU Ase US Wie Nac GG) SO Od (A.17) 


Start by letting A, = uj, A uj, A... A Uj, be a pseudoscalar for the vector space, then from 
equation 1.61 we have 


ul = (-1)" 1 A... A ty A... AtmAz?. (A.18) 
Since u” is in the space defined by A,, we have uA, = u" - A,. Now consider the progression 
utA, =u- An =(-1)2 uA... Ati, A... AUn (A.19) 
Cu? A u') A, = (a? A u) -A, 
=u". (u"- A,) from equation A.6 (A.20) 
= (-1)"* u®- (up A... Ad A... A tn) (A.21) 


= (-1)9) SO (1) hu ug (ur A. A tha Ae At A At) (A.22) 
k#iy 

Sy Se iy he Na nae Oi, (A.23) 
Equation A.23 is obtained by applying equation A.5 to equation A.22. The reason that in 
equation A.23 the second power of —1 is 72 — 2 and not ig — 1 is that in equation A.21 u;, has 
been removed from the wedge product. Since 

(a Kawa") A,= Cu" Axes. NU") A, 
— (u’* ix (k= dN u")) A, 


= (Cts JNU) A! (A.24) 

so by induction 
(ul? A... Au?) Ap = (DEF a AL At A. Ata, A A tn (A.25) 
ui A... Au = (-1)E =O (up AL. A tig Ao. A ty, Aw. AU) Az! (A.26) 
ur AL. Aut = (<1) (12218 (ay AL A tig Ao At Ae Attn) ATL (A297) 
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Now let the indices {7,41,...,%n} be the complement of the indices {7),...,7,} in the set {1,...,n} 
with the condition that 0 < i,4, < +++ <i, <n. Then we may write equation A.27 as 


r(r+l1) 


ut A... Aut =(-1)? (-1)2" (uy, A... Au, Ag? (A.28) 


n 


—1) (it Anca da OP Ne AW) 

—1)7" (-1)%= Hug A... Ati) « (tig A. A Uy,) Az?) 
1) FD) (as A. Adi)» (thingy A. A Uj,) Az?) 
) 


—1)E FD (4, AL. Ady) A (tig Ae AUj,)) Ant 
4(-1 . 

12 GD a A A ti A Ugg Noe AUy,47) — (A.29) 
Ji sd jr 
PO ac0s (A.30) 
The step from equation A.29 to equation A.30 requires a detailed explanation. First the r.h.s. 
of equation A.29 is zero unless there are no repetitions in the index list [7,,..., 71, jr4i,--+5 In]; 
but there are no repetitions if aes ... 6)" #0. Because both the %;’s and the j;’s are ordered 
we must match 7; = 7; V1 <1 <r so that in the Kronecker delta’s the index subscripts match. 
Also note that in equation A.29 that the subscript list |7,,...,71,Jp41,---;Jn] is not in ascending 
order so that if not zero (u;, A... A Ui, AU;,4, A... A Uj,) Ay? = £1. The factor (—1)2i=1 GY) 
is required to order the subscript list. For the nonzero products we have 

IN crc NU) UP Nae ae = (—1)" 7 uy A. A(—1)2 ay, A Ui, A---AUj, AZ). (A.31) 
Now remember that 7; 1 < 1] <r is the absolute position index of a basis vector in the pseudoscalar 
product. To move uj, in equation A.31 to its correct position in the pseudoscalar product requires 
jr — 1 transpositions and hence the factor (ye sy Repeating the proceedure for basis vectors 
Uu;,_, to u;, reduces the r.h.s. of equation A.31 to 1. 


Appendix B 


BAC-CAB Formulas 


Using the python geometric algebra module GA! several formulas containing the dot and wedge 
products can be reduced. Let a, b, c, d, and e be vectors, then we have 


a: (bc) = (b-c)a—(a-c)b+(a-b)c (B.1) 

a:-(bAc) =(a-b)c—(a-c)b (B.2) 
a-(bAcAd) = (a-d)(bAc) —(a-c) (bAad) 

+(a-b)(cAd) (B:3) 

a:-(bAcAdAe)=-—(a-e)(bAcAd)+(a-d)(bAcAe) 

—(a-c)(bAdAe)+(a-b)(cAdAe) (B.4) 
(a-(bAc))- (de) =((a-¢) (6-e) —(a-b) (c-e))d 

+ ((a-b) (cd) (a0) (b-d))e. (B5) 


If in equation B.2 the vector b is replace by a vector differential operator such as V, 0, or D (we 
will use D as an example) it can be rewritten as 


(a-D)c=a-(DAc)+(a-é)D 
=a-(DAc)+D(a-é) 
=a:(DAc)+D(é-a) (B.6) 


Cyclic reduction formulas are 
a:-(bAc)+c:-(aAb)+b-(cAa)=0 (B.7) 
‘Alan Macdonald website http: //faculty. luther .edu/~macdonal/vagc/index.html 
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a(bAc)—b(aAc)+c(aAb) =3aAbAc (B.8) 
a(bAcAd)—b(aAcAd)+c(aAbAd)—d(adbdAc) = 

4aNbAcAd (B.9) 


Basis blade reduction formula 


(aA b)-(cAd) = ((aAb)-c)-d 
= (a-d)(b-c) —(a-c)(b-d) (B.10) 


But we also have that 
((aA b)-c)-d=(a-d)(b-c) —(a-c) (b-d) (B.11) 
which gives the same results as equation B.10. Since for any bivector blade B = a A c we have 
(B-c)-d=B-(cAd). (B.12) 


By linearity equation B.12 is also true for any bivector B since B is a superposition of bivector 
blades. 


Finally one formula for reducing the commutator product of two bivectors 


(aA b) x (cA d) =(a-d)bAc—(a-c)bAd 
+(b-c)a\d—(b-djaNc (B.13) 


Appendix C 


Reduction Rules for Scalar Projections 
of Multivectors 


If A is a general multivector define the even, At, and odd, A~, components of A by 


A® = (A)y + (A)o +> (C.1) 
AS =A) eA te (C.2) 


Then the following rules are used for the reduction of the scalar projections of multivector 
expressions. 


Reduction Rule 0 (RRO) (A) =0. 
Reduction Rule 1 (RR1) Let A =a,---a,, then ifr is even A = A* and if r is odd A = A~ 


If r = 2 we know that A = At. Assume RR1 is true for r even and multiply the result by 
a vector a then 


Aa = ((A}g CAs tear) @ 
=(A)at ((A)pa), + ((A)ga)g to (C.3) 


Now assume that A = A™~ and multiply by a to get 


Aa = ((A), + (A); feos) a 
= ((A)1 a) + (A), @)9 + ((A)3 @)o + (Ag a)g toe: - (C.4) 
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Reduction Rule 2 (RR2) For any two general multivectors A and B 


AtBt = (At Bt) (C.5) 
AB” =(AB’)" (C.6) 
A’ Bt = (A Bt)” (C.7) 
AT BS {ATB-) = (C.8) 


AT = 97, (A)o; and AT = 97, (A)o,_1 and likewise for BY and B~. Thus the general term 
in At B® if of the form 


Cs) ee (Ao (Bas) = ((Aa (Bas), pear 


+ (Ala (Bbai) (C9) 


which has only even grades. The general term in A*B™ if of the form 


(Aa (Baja = (Abas (Baia) se pga) + (Adar (Bex) wra-oycg t 
+ ( (Ao; (Bay) (C.10) 


I2-+5)-1| 


which has only odd grades. The general term in A~ B™ if of the form 


(A)»; (B)os_4 -_ (Aas Baa)» a (Aoi (Bait) yeas 


+ ((A)oia (Baj-1) (C.11) 


|2i-+3)—-2| 
which has only even grades. 
Reduction Rule 3 (RR3) If B, is an s-grade multivector ((A), Bs) = 6,5 (ABs) 
The lowest grade in (A), B, is |r — s| so if r # s there is no zero grade. 
Reduction Rule 4 (RR4) If A, is an r-grade multivector (A,B) = A, * (B),, 
We have from RR3 that (A,B) = (A, (B),.) = A, * (B),. 
Reduction Rule 5 (RR5) (A;...A,) = (Aj Aji... ApAr.-- Aj—1) 
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Using equation 6.17 we have 
U Aa Bigs ces A oA aae ole aN At) 
= (Ay (Ar Ag... Ag-1)) 
(A (ArAd--. An-a))') 
= (A; (Aj Ag... Ag—1)) - (C.12) 
No just keep applying equation C.12. 
Reduction Rule 6 (RR6) (AB) = (AtBt+A™B-) 
A*B~ and A~B* are odd and have no zero grade. 
Reduction Rule 7 (RR7) ((B),-a)-b=(B),-(aAb) 


Since (B), is a grade 2 mulitvector it can be written as the linear combination of bivectors 
c; \ d;. Now apply equation B.10 in Appendix B directly 


((c; A dj)-a)-b = (ce; Adj): (a AB). (C.13) 

In applying the reduction rules one must be carefull of the multivector nature of the multivector 
derivative. Note that is ~ is a vector then 

1. Vy is an even multivector (Vw = (Vv), + (VY),). 

2. Oy is a vector (Oy = (Oy),)- 

3. Ovy is an even multivector (Ovy = (Ovu)y + (Ovu)s): 
If w is a spinor (even multivector w = ~t) then 

1. Vw is an odd multivector (Vi) = (Vw) ). 

2. Oy is an even multivector (Oy = (Oy)*). 


3. Ovy is an odd multivector (Ovy = (Ovy) ). 


Appendix D 


Curvilinear Coordinates via Matrices 


To use matrix algebra to transform from linear to curvilinear coordinates start by defining the 
matrices 


g = [ua Uy] = [gig] (D.1) 
g = [uv -v’] = [97] De) 


where the u;’s are a set of fixed basis vectors for the vector space. Then define the coordinate 
vector x by the coordinate functions * (0) 


x (0) = 2' (0) u; (3) 


where @ is the curvilinear coordinate tuple 0 = (0;,...,4,). In the case of 3-D spherical coordi- 
nates 0 = (r,0,@) where @ is the azimuthal angel measured in the xy plane and ¢ is the elevation 
angle measured from the z axis. Then the position vector, x, in terms of the coordinate tuple, 


(r,0, 0), is 
x =r (sin (A) (cos (¢) uz + sin (d) uy) + cos (8) uz) (D.4) 


The unnormalized basis vectors are then given by 


— Ox 
‘ 6e" 


(D.5) 


If we represent x as the row matrix 


x = [x (6,...,0")] (D.6) 
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then a3 
a 
and the matrix J = [Jj;] is 
J = Jacobian (2). (D.8) 


Again for spherical coordinates (J is the transpose of the Jacobian matrix) 


sin (?)cos(¢) — sin(@) sin (@) cos (0) 


J’ = | rcos(@)cos(¢) rcos(@)sin(¢) —rsin (6) (D.9) 
—r sin (¢) r cos (@) 0 
and 


The square of the magnitudes of the basis vectors are 
ee, = (ite) : (Ji tes) 


= Jie One 


[ei6ij] = diag (J" 9J) 
[lei| 5:j] = \/ diag (J*gJ) 
. Define the normalized matrix a = |a;;] by 


where |e;| = sgn (e?) \/|e?| so that |e;| could be positive or negative depending on the signature 
of g and the actual position vector. We now can write 


é; = AijU;. (D.16) 


For the reciprocal curvilinear basis vectors é’’s we would have 


é' (0) =a" (0) w’. (D.17) 

To calculate a”? note 
e* €;, = a oj “Un = 6 (D.18) 
= Gap. 6" = 0) (D.19) 


251g Meee O) (D.20) 
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or in matrix terms 


aa’ =1 (D.21) 
a= (aT) (D.22) 
a= FaraTyiteate (a’) é. (D.23) 


We also have 


Gi; = 6; + 6; (D.24) 
= Os 40 lia te (D.25) 

= Oi Dann (D.26) 

G = aga" (D.27) 
GI=é.é (D.28) 
= a'™oJ™y™ - u™ (D.29) 
oaks aaa a (D.30) 
G=aga (D.31) 


A 


0é; 
The final items required for the curvilinear coordinate transformation are the ——’s which are 
k 


A 


0é; 
needed for the calculation of geometric derivatives in curvilinear coordinates. The —’s are 


OO; 
simply calculated by noting 
Oe; . 
Am 0é; Am A m 
59, — b Hse »€; = Pig; = Tpim (D.33) 
Oa; Oa; 
kim = 0? uP ug = PP D.34 
RE Oper dO, 4 De 
Oa; 
=o ts D.35 
On ag, (D.35) 


i een _ 001 (D.36) 
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Again for spherical coordinates we have for the I’’s 


and 
OG, O€g _ Oey -. 
ar ~° or or” 
Dep _ i OCp:. 2% Oég _ 
ao 06° a9 
has = cos (6) é4 Oe = —sin (6) é4 065 = —cos (0) é, + sin (8) é9. 


Og 


(D.37) 


(D.38) 


(D.39) 


(D.40) 


Appendix E 


Practical Geometric Calculus on 
Manifolds 


The purpose of this appendix is to derive a formulation of geometric algebra/calculus that is suit- 
able for use in a computer algebra system. For this purpose we must reduce or delay the divisions 
required for the calculation of quatities in the algebra and calculus and to also precalculate those 
quantities required for the acts of mulitplication and differentiation. 


Start with a vector manifold defined by x (@) where @ is a n-tuple, @ = (0',...,0”), and w is a 
vector function of @ in a vector space of dimension > n. Then a basis for the tangent space to 
the manifold is defined by 


Ou 
ei = api (E.1) 
and the metric tensor by 


gi; (9) is all that is needed to define the geometric algebra on the tangent space of the manifold. 


To define the intrinsic and covariant derivatives we need to calculate the reciprocal basis to the 
tangent space basis. To do so we use the relations 


“—— 
esl 
w 

er 


e; = ge? 


e'=g'e; (E.4) 


ZU 
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where g‘! is the inverse of g;;. Additionally, we have from equation 5.167 


n(n — 1) 
E2 =(e,A...Ae€n)°=(-1) 2 det(g). (E.5) 
where det (g) is the determinant of the metric tensor and E,, = e; \.../A €,. Now define 
é = (-1)"(e,A...A8;A...A€n) En. (E.6) 
Then 
yee a (E.7) 


ER (=1)° Y?? det (9) 
So that calculating é’ requires no division. 
To define the geometric calculus of the manifold the projections of the derivatives of the tangent 


vectors onto the tangent space of the manifold is required to define the covariant derivative of a 
multivector field on the tangent space of the manifold. 


The projection of the derivative of the tangent vector into the tangent space is the covariant 
derivative of the tangent vector, D;e; so that 


e; Oe; 
Die; = ef (« : at) = eg!" (« ; at) (E.8) 
and ; F 3 
D;e; = De; = e* (« . oa) E.9 
J yee det @ j mse eae det (g) x OO" ( ) 


so that D;e; can be calculated without any divide operations (other that those involving caclu- 
lating the derivatives of the coefficients of the tangent vector derivatives). We have identified a 
common divisor for all the covariant derivatives of the tangent vectors. 


de; 
The next step is to compute D;e; we must compute e;, - aa First differentiate g,; with respect 
nee Agi; 8 de; 0 a 
Gij = ej ej _ (An ej 
0k = 56n ©: T Ci Dok a0! “Ej Ci” Bk (E.10) 
Oe; Oe; 
Where the final form of equation E.10 results from the vector manifold relation aa = oR . Now 
cyclicly permute the indices to get 
Ogn Oe; Oe; 
Hoes adh goat E.l1 
aor ape ER BB oe) 
Ogjx OE; Jex 
= + 7 =, E.12 
api aps SRT I” HB oy 
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Now add equations E.10 to E.11 and subtract E.12 and divide by 2 to get 


i (54 Ohi _ “Bi ee OE es One (E.13) 


2\ 00k ' a6) ae) = Abk 


Thus 
Die; = ek (E.14) 


which can then be pre-calculated for all combinations of basis vectors and coordinates and the 


results placed in a table or dictionary. Note that while 0, = is a scalar ooperator, D; and 


: ok 
D,, are linear operators, but not scalar operators as can be seen explicitly in equation E.14. 


Thus the manifold can be defined either by specifying the embedding function, #(@), or the 
metric tensor, g;;(@). For computational purposes if a (@) is specified g;; (@) will be computed 
from it and Q;;, computed from g;; (8). 


Now consider an r-grade multivector functon F; on the manifold given by (using the summation 
convention) 


BS, KLINE (B.15) 


where both F“*" and the e;, can be functions of 6, the coordinates and that i; < --- <i) < 
-++ <a, V1 > %, >n. The the intrinsic derivative of F, is defined by (summation convention 
again for repeated indices) 


OF, = e30;F, (E.16) 
= OF" ei (e;, N...NG;,) + Fe) Ss (Cis Neve Kh OCx Non JVs): (E.17) 
i=1 


The covariant dervative, D is defined as the projection of the intrinsic derivative onto the tangent 
space of the manifold (the tangent space can be defined by the pseudo scalar E,,) 


DF, = Pz, (OF). (E.18) 


Thus 


DF, = 0; Fo @) (E;, awed €;,.) + Ft ed Se (E;, Dextre D7; ae'silS €i,.) ; (E.19) 
l=1 
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For computational purposes 


il 


ee O,F re) (e,, A... A i, 
are dae (g) J ( ) 
4 Fretr @ (ei, [Ntctecd Dje; Neitacl\ e;,) (B.20) 
det (g)° s | 
1 ae 
a O,F"-% @ (E;, \ ane, \ €;,) 
(—1)"°-Y? det (g) 
1 os ty 
5 Feral S~ (6, A... A (Que) A... ei,) (E.21) 
5) 41 J MP 
det (9) l=1 


Again for computational purposes D;e; will be calculated in terms of a linear combination of e;’s 
and then the connection multivectors computed from ordered combinations of 7,...i, 1<r<n 


Dis O,F-"@ (e,, A...N e;,) 


(1)? det (g) ” 
POC Gs cote te (E.23) 


det (g)” 


It is probably best when calculating DF, (or DF = D(F))+-:-+ D(F),,) to have the option 
of returning (—1)"-)? O;F8-e) (e;,, A... Ae;,) and F"""C {i,...i,} separately since the 
det (g) scaling for each term is different. 


To calculate the outer and inner covariant dervatives use the formulas 


DAF =(D(F)o), +--+ (D(F)y4),, (E.24) 
DSi = DE 3) ch ee) pa (E.25) 


Appendix F 


Direct Sum of Vector Spaces 


Let U and V be vector spaces with dim (U) =n and dim(V) = m. Additionally, let the sets of 
vectors {u;} and {v;} be basis sets for U and V respectively. The direct sum, U 6 V, is defined 
to be 

USV={(u,v)|ucU,veV}, (F.1) 


where vector addition and scalar multiplication are defined by for all (a, bi), (a2,b2) Ee U @V 
andae Rr 


(a1, b,) + (a2, bz) = (a1 + ae, by + ba) (F.2) 
a (a1, b,) = (aa, ab;). (F.3) 


Now define the maps iy : U ~ U @V andiy :V ~U®V by 


Thus (u,v) = iy (uw) + iy (v) and the set {iy (uj) ,iy (vj) |0<i<n,0 <j < m} form a basis 
for U ®V so that if a € U @V then (a’, 8’ € RK) and using the Einstein summation convention 


w = aliy (us) + Piv (0). (F.6) 


If from context we known that 2 € U @ V and we are expanding x in terms of a basis of U @ V 
we will write as a notational convenience wu; for iy (w;) and v; for iy (v;) so that we may write 
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Also for notational convenience denote for any x € U @V, ry = iy (a) and zy = iy (a). 


Likewise we define the mappings py: U@V > U and py: U@V > V by 


pu ((u,v)) =u (F.8) 
v. (F.9) 


For notational convenience if py ((u, v)) = 0 we write u = py ((u, v)) and if py ((u,v)) = 0 we 
write v = py ((u,v)). We always make the identifications 


u + (u, 0) 
v + (0,v). 


Which one to use will be clear from context or will be explicitly identified. 
One final notational convenience is that for (u,v) € U @ V we make the equivalence 


(u,v) Out. 


Appendix G 


sympy/galgebra evaluation of GL (n, #) 
structure constants 


The structure constants of GL (n, R) can be calculated using the sympy python computer algebra 


system wiht the latest galgebra modules (https://github.com/brombo/sympy). The python 


program used for the calculation is shown below: 


#Lie 
from 
from 
from 
from 


Algebras 

sympy import symbols 

sympy.galgebra.ga import Ga 

sympy. galgebra.mv import Com 
sympy.galgebra.printer import Format, xpdf 


#General Linear Group E generators 


def E(i,j): 
global e,eb 
B = e[i]*e[j] — eb[i]*eb[j] 


i 
Bstr = ’E{’+i¢+j 
print ’%’ + Bstr + ’ ‘ 
return Bstr, B 


#General Linear Group F generators 


def F(i,j): 
global e,eb 
B = e[il*eb[j] — eb[i]xe[j] 
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v] 


Bstr = ’F{’+i¢+j 
print %’ + Bstr + ’ 
return Bstr, B 


| + 


} 
* B 


#General Linear Group K generators 
def K(i): 
global e,eb 
B = e[i]*eb[i] 
Bstr = “K{’ +i+ 7}? 
print ’*%’ + Bstr + ’ =’,B 
return Bstr, B 


#Print Commutator 
def ComP(A,B): 
AxB = Com(A[1] ,B/1]) 
AxBstr = ’*%’ + A[O|] + r’ \times ’ + B[O] + 7? =’ 
print AxBstr, AxB 
return 


Format () 
(glg ,ei ,ej ,em,en,ebi,ebj ,ebm,ebn) = Ga. build(r’e_i e_j e_k e_] \barf{e}_i \bar 


6 Sf 1 er fre) y. ko een, lent 
eb = {’i’:ebi,’j’:ebj,’k’:ebm,’1’:ebn} 


print r’#\centerline{General Linear Group of Order $n$\newline}’ 
print r’#Lie Algebra Generators: $1\le i < j \le n$ and $1 \le i <1 \le n$’ 


#Calculate Lie Algebra Generators 
Bij = B01 3?) 
By Se) 


Ki = K(’i’) 
B= Fa Sa) 
Piss PC er) 


print r’#Non Zero Commutators’ 
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#Calculate and Print Non—-Zero Generator Commutators 
ComP(Eij , Fij) 

ComP( Eij , Ki) 
ComP(Fij , Ki) 
ComP( Eij , Eil) 
ComP( Fij , Fil) 
ComP( Fij , Eil) 


v] 


xpdf(paper=’ letter’ ,pt=’12pt’ ,debug=True, prog=True) 


Only those commutators that share one or two indices are calculated (all others are zero). The 
TREX output of the program follows: 


General Linear Group of Order n 
Lie Algebra Generators: 1<i<j<nandl<i<l<n 
Fig = 4 N€j — & AG; 
B= Cee eed, 
K;,=6e; A é; 
Fy =e, \e,—-€; A & 
Fy =e,N& +e &; 


Non Zero Commutators 
Ey Pp = Lea \ 6; = 267K 6; 


Hig Xk he Hh =e, Ne; 
Fig xX Ky = —e; Ne€3 + & AG; 
fy x Ly = eg Ne +e AE; 
Eek Ey = GING 8 Ne 
Le Lag CPN. Cy CIN; 


The program does not completely determine the structure constants since the galgebra module 
cannot currently solve a bivector equation. One must inspect the calculated commutator to see 
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what the linear expansion of the commuatator is in terms of the Lie algebra generators. For this 
case the answer is: 


Ei; x Fy =2(K; — K;) (G.1) 
By, Xx Ky=—Fy (G.2) 
Fy x K;, = —Ey (G.3) 
Eo hy k (G.4) 
Fy; x Fy = Ej (G.5) 
Fy; x Ea = Fy. (G.6) 


Appendix H 


Blade Orientation Theorem 


A blade only depends on the relative orientation of the vectors in the plane defined by the blade. 
Since any blade can be defined by the geometric product of two orthogonal vectors let them be 


e, and e,. Then any two vectors in the plane can be define by: 


A= Azz + Ayey 
b= Dpex yey 


and any rotor in the plane by 


R= ab = (a-b) + (azby — aybz) Cxey 


as long as 
RRi=1 
but 
Rezey = €ze,R 
and 


Re,R'Re,R!' = Re,e,R' = €,€,RR' = €7€y 


and absolute orientations of e, and e, does not matter for e,e,. 


Ld 


Appendix I 


Case Study of a Manifold for a Model 
Universe 


We wish to construct a spatially curved closed isotropic Minkowski space with 1, 2, and 3 spatial 
dimensions. 


To do this consider a Minkowski space with one time dimension (e? = 1) and 2, 3, or 4 spatial 
dimensions (unit vector squares to —1) as indicated below 


2 2 2 2 2 
€g9 = —e] = -—eg = -e3 = —Eg = 1 (1.1) 


The vector manifolds will designated by 


Vector Function Spatial Dimensions Coordinates Components 


x® 1 TP €0, €1, €2 
xX) 2 T, p,9 €0, €1, €2, €3 
xe) 3 T,p,9,0 €0, €1, €2, €3, €4 


The condition that 7 does parametrize the time coordinate of the manifold is that 


Ox ax® ‘ 


OT OT (1.2) 


eee = 
since then the time coordinate is given by (all spatial coordinates are fixed) 
[vee dr =F (1.3) 
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Then the manifolds in 1, 2, and 3 spatial dimensions are defined by 


XY = t(r)eg +r (rT) (cos (E)er + sin (E)es) (1.4) 
X® =t(r)e9 +r (r) (cos (E)er + sin (£) (cos Bey + sin des)) (1.5) 
X®) = t(r)egt 

r (rT) (cos (Eyer + sin (£) (cos Oe2 + sin 6 (cos ge3 + sin es))) (1.6) 


Where r (7) is the spatial radius of the manifolds. In order for equation 1.2 to be satisfied for all 
the three manifolds we must have 


ele a 


1 
For the specific case of r (7) = oe the derivatives and t (7) are given by 


Figure I.1: Expansion Radius of Model Universe, r (7), where 7 is the time coordinate. 


dr 
ema 

dt 

— a af ty? 
adr 


ti : (rv 1+7?2+sinh' r) 


215 


Spatial coordinates for the model universes are given by p a spatial radius and the angular 
coordinates @ and ¢ if required. A visualization of the 1-D spatial dimensional manifold is shown 


in figure I.2. e, and e, are 
ax) 


ax] 
= an Cp = ae 


OT 


and metric tensors for the three cases are 


e; 


w _ 9X9 ax 


ar, nr 


where p,v = {7T, p,0, d}. If we define h (7, p) as 


Then the differential arclength given by the metric tensor, gj, is 


(ds)? = gy da"dx” 


and the metric tensors for the three cases are (using sympy to do the algebra) 


1-h? h 
m=" A) 


1-h? A 0 
9?) = h —1 0 ; 
0 0 —- (rsin (£)) 
1-h? A 0 0 
h —1 0 0 
9) = 0 0 — (r sin (2)) 0 


0 0 0 _ (rsin ( 


) sind) 


(1.8) 


(1.10) 


(1.11) 


(1.12) 


(1.13) 


(1.14) 


216 APPENDIX I. CASE STUDY OF A MANIFOLD FOR A MODEL UNIVERSE 


Tg 


(ae) 


aSIOAIUY) GT ‘ZT olnsty 


If we renormalize eg and eg to be unit vectors 


The metric tensors gi) and gs) become 


1-h2 bh O 

oS h -1 0 

0 Ot 
PR he 0G 
18) — ke <i 70-. °0 
uv 0 -1 0 


and 
det (gf:)) = det (90) = det (gf?) =-1= 1 


For the 1-Dimensional space the differential arc length is 
(ds)? = (1 —h?) (dr)* + 2hdrdp + (dp)’ 
so that for the light cone, ds = 0, we have the differential equation 


dp dp . 
=> 2 | | — 
walls (@) oo 


dé 
Note that this equation also applies to the 2 and 3 dimensional case if we set — = 


Solving for ae gives 
dt 
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(1.15) 


(1.16) 


(1.17) 


(1.18) 


(1.19) 


(1.20) 


(1.21) 


dp _ 
ae 


(1.22) 


Using the integration factor for linear first order differential equations the solution to equation 


1.22 is (p (7) = 0) is 
- 


p(r)= arr) f So 


TO 


(1.23) 
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Note that r (7) and ar (rT) have the same solution p(7). Now consider the case that r (7) = 7", 
then 


=, +rln (=) 


p(t) = ent | (r')"dr' = de 1 ey (1.24) 


TO 


Typical p(r)’s for various —1.5 < 7 < 1.5 are shown in the light cone plot. If 7 > 0 the speed of 
light is greater than c in flat space. If 7 < 0 the speed of light is less than c in flat space. Note 
that if the universe is curved, but not expanding or contracting the speed of light is the same as 
c in flat space. 


Photon Path (" (T)= (<)’) 


1 f 1 
1.0 1.5 2.0 2.5 3.0 
Te 


TO 


Figure 1.3: Lightcone in Curved Space 


2 
For a time periodic universe r (7) = sin (=) where T' is twice the period of the universe and 
p(T) =0, then 


CSC 


p(7) =T sin (=) In 3 uae A (1.25) 


r) +e (=) 
CSC T co T 


The plot of equation 1.25 is shown in figure I.4 
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Photon Path (r (r) = sin (4)) 


Figure I.4: Periodic Light Cone 
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I.1 The Edge of Known Space 


Another kinematic question to answer is under what conditions light cannot access parts of the 
universe. The critical quantity is 


\(7) = ran (1.26) 


where J (the angular distance around the universe restricted to 0 < J < 7m) is the measure of 
how far from the observer you are. Since the universe is spatially periodic the maximum value 
of J is x. If A(7) > 1 for some finite 7 you can access the distance defined by ¥. Substituting 
equation [.23 into equation 1.26 gives 


Te if? sar 
Alt )= aE rr) (1.27) 


T dr’ 
: eo 0. (1.28) 


First consider a linear expansion model of the form 


romm(vee(5-9) a 


To 
1 «o(>) 
e \To/ —1 (1.30) 


so that the connection condition is 


where r (7) =7o. Then 


In a linearly expanding universe the photon time of flight increases exponentially with distance. 
Now consider super-linear expansion of the form (7 > 1) 


r(r) =o (1 +a (= = i) ) (1.31) 


1 

ae) 

ans | ; BSG (1.32) 
0 0 b+ yw 

but the integral in equation I.32 does not have a closed form solution unless we let T — oo. In 


that case! we can write : ; 
- a nese (=) 2 (<) (2) (1.33) 
1) 1) T TO 


A contour plot of the left side of equation I.33 is shown below 


Then 


‘fo 7Zx is 3.241-2 in ”Gradshteyn and Ryzhik” 


1. THE EDGE OF KNOWN SPACE 221 


1.5 

1.0 
1.44 

0.8 
1.3 

0.6 

> = 

0.4 

1.2.4 
fF 0.2 

1.14 

0.0 
1.0 T T T T T T T 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 


Figure I.5: Photon Exclusion Zone 
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0 
The right side of equation 1.33, (<) (2), is interpreted as follows - 
as To 


0 0 
1. — is the fractional distance around the closed spatially periodic universe. — = 1 is as far 
1 


as one can go before the distance from the observer starts to decrease. 


ih Tr 
2. — is a measure of inflation. Immediately after an inflationary epoch ae 
TO i) 


0 
Thus equation I.33 determined the maximum distance — that a photon can propagate in a finite 
‘3 


amount of time. 


: é : ie TO ay. i 
Another question to consider is under what conditions — will increase as 7) increases or when 


TO 
will the following be true 
r (7) > r (To) ; (1.34) 
cm To 
Equation 1.34 is equivalent to 
n-1 
a (= = i) Si (1.35) 
TO 
. a ee . 
so that if 7 > 1 then —— will eventually grow as 79 increases. 


TO 
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